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X XIX. BAND ZWEITES HEFT 1960 


Druckentnahme durch Lingsschlitze 
Von F. Krause 


I, Einleitung.' Die vorliegende Arbeit behandelt die Messung des statischen Druckes in einer 
strémenden Fliissigkeit. Seit den grundlegenden Untersuchungen von Prandél und seinen Mitarbeitern? 
benutzt man zur Druckentnahme schlanke Stromlinienkérper, deren Oberflache eine oder mehrere 
Druckentnahme-Offnungen besitzt. 

Zwischen dem so gemessenen Druck und dem Druck, der am Ort der Wandéffnung vorhanden war, 
bevor der Stromlinienkérper eingefiihrt wurde, bestehen nun verschiedene Abweichungen, die sich auf : 
drei Ursachen zuriickfiihren lassen. nee 

: Die erste besteht in der Verdrangung der Strémung durch den eingebrachten Rotationskirper. Sie 
kann mit Hilfe der Singularitatenmethode berechnet und durch geeignete Formgebung weitgehend 
vermieden werden* 4, ‘ 

Die zweite besteht in der Rauhigkeit des Bohrlochrandes® und Welligkeit® der festen Wand in der i 
E ueepong der Druckentnahmeéffnung und abt sich durch sorgfaltige Herstellung weitgehend ver- a 
meiden. 

Die dritte Ursache schlieBlich besteht in der Sekundarstrémung, die sich iiber und in der Wand- 
6ffnung ausbildet. Fiir rande Bohrungen hat man den Einflu8 untersucht, den der Durchmesser der 
Bohrung” *9, die Tiefe der Bohrung®, die Neigung der Bohrlochachse®, Abrundung des Bohrlochran- 
des?°, die Parameter der auBeren Strémung in unmittelbarer Umgebung des Bohrloches" und ein Druck- 

gradient der AuBenstrémung” auf die Druckentnahme ausiiben. Die beobachteten Abweichungen be- 
tragen gréBenordnungsmaBig 1% des Staudrucks. 

Die vorliegende Arbeit behandelt ebenfalls den Einflu8 der Sekundarstrémung auf die Druckent- 
nahme und ist aus einer Fortfihrung der Arbeit von A. K. Ray" entstanden. Im Gegensatz zu den 

-bisherigen Untersuchungen wurde die Wandéffnung nicht als Bohrung, sondern als langer Schlitz 
_ausgebildet. — 
Theoretische Uberlegungen lassen vermuten, daB ein langer Schlitz die Druckverteilung, die sich 
iiber ihm einstellt, mittelt. Ist der Druck der urspriinglichen Strémung in Strémungsrichtung kon- 
stant, so wird der Druck im mittleren Teil des Schlitzes voraussichtlich ziemlich stérungsfrei tiber- 
tragen. Die Sekundarstrémung in der Schlitzéffnung erzeugt nur an der Vorder- und Hinterkante 
_ Storungen der urspriinglichen Druckverteilung. Es ist daher zu erwarten, daf der Druckentnahmefehler 
um so kleiner wird, je langer der Schlitz ist. : 

Andererseits tritt bei einem langen Schlitz der EinfluB eines Druckgradienten in Strémungsrichtung 
-starker hervor. Zwischen den beiden Schlitzenden besteht dann ein Druckunterschied, der bei gege- 
benen Druckgradienten um so grofer ist, je langer der Schlitz ist. Es fragt sich nun, an welcher Stelle 
der Schlitzlange der veranderliche AuBendruck mit dem Druck im Schlitz tibereinstimmt. Da diese 


1 Gekiirzte Fassung der gleichlautenden Diplomarbeit der Math.-nat. Fakultat der Universitat Gottingen 
_vom Dezember 1956. 
2-H. Kumbruch, Messung strémender Luft mittels Staugeraten, VDI-Forschungsheft 240 (1921). 
3 EF. Owen u. F.C. Johansen, The design of pitot-static-tubes, ARC, R. u. M. 891 (1925). 
: 4 W. M. Schulze, G. C. Ashley, jr. u. J. R. Erwin. Several combination probes for surveying static and total 
pressure and flow direction, NACA TN 2830. 

8 J. Polzin, Ing.-Arch. 10 (1939), S. 326. ees 

; 6 H. Peters, Druckmessung. Handbuch der experimentellen Physik, Bd. 4, Teil 1, Leipzig 1931. 
-- % G. Fuhrmann, Theoretische und experimentelle Untersuchungen an Ballonmodellen, Dissertation, Got- 
_tingen 1912, S. 46. 2 

8 J, Fage, Proc. Rey. Soc. London 155, (1936), S. 576. 
9 A, Myadzu, Ing.-Arch. 7 (1936), S. 35. 
10 R, Herrmann, Experimentelle Untersuchung zum Widerstandsgesetz des Kreisrohres bei hohen Reynolds- 

chen Zahlen und groGen Anlauflangen, Dissertation, Leipzig (1930). : 
11 4, K. Ray, Einflu8 der Bohrlochgro®e auf die Anzeige des statischen Druckes fiir verschiedene Reynolds- 
ahlen, Dissertation, Gottingen (1955). 

12 W. Lincke, Physik. Z. 32 (1931), S. 905. 
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Stelle zundchst nicht bekannt ist, wurde bei den Experimenten der AuBendruck in Schlitzmitte mit 
dem angezeigten Druck verglichen. Aus der hier beobachteten Druckdifferenz kann dann bei bekann- 
tem Druckgradient der AuBenstrémung die Stelle des Druckgleichgewichtes berechnet werden. 
Im Folgenden wird nun zunichst der theoretisch zu erwartende KinfluB der Sekundarstrémung im 
Schlitz auf die Druckentnahme fiir eine Str6mung ohne und mit Druckgradient qualitativ geschildert. 
AnschlieBend wird eine Versuchseinrichtung beschrieben, mit der die Stelle des Druckgleichgewichtes 
in Abhangigkeit von der Reynoldschen Zahl bei verschiedenen Druckgradienten ermittelt werden kann, 
und im letzten Abschnitt wird das experimentelle Ergebnis mit den theoretischen Vorstellungen ver- 


glichen. 


2. Verhalten bei konstantem Druck. Im Innern des Schlitzes wird sich infolge der Reibung eime | 
stationare Strémung ausbilden, wie sie schematisch in Abb. 1 gezeigt ist. Der Antrieb durch die Rei- | 
bungskrafte der auBeren Strémung wird ausgeglichen durch die Haftreibung an den Seitenwanden des 
Schlitzes. In den tieferen Teilen des Schlitzes muB sich dann aus Kontinuitatsgriinden eine Rickstré- 
mung ausbilden. Im mittleren Teil des Schlitzes erhalten wir also etwa ein Geschwindigkeitsprofil wie 
es in Abb. 2 gestrichelt eingezeichnet ist. 


Vorderkante, Totwasser Minterkante, Stougebiet 


—-— Trennungstiiche 


— — — Ceschwinaigkeitsverteilung 


Geschwindigkeitsverteilung 
bei Annahme einer 
Trennungstlaohe 


es <.rw ruhende Flissigheit unter 
“252 der Irennungstiache %, 


Abb. 1. Druckstérungen an den Schlitzenden. 


Abb. 2 (rechts). Geschwindigkeitsverteilung im mittleren Teil 
des Schlitzes. 


Wir denken uns nun die Schlitzvorkammer sehr lang und breit. Mit wachsender Tiefe klingt die 
Rickstrémung infolge von Reibungswiderstanden immer mehr ab, so daB schlieBlich die Flissigkeit 
ruht. Dabei wird die kinetische Energie, die die bewegten Fliissigkeitsteilchen aus der AuBenstrémung | 
erhalten haben, in Warme verwandelt und kann daher nicht zu einer Anderung des hydrostatischen — 
Drucks weit unten in der Manometer-Leitung beitragen. Der durch die Wandéffnung entnommene 
Druck mu8 dem in der Wandéffnung gemittelten entsprechen. | 

Im mittleren Teil des Schlitzes stimmt der Druck mit dem konstanten Druck der ungestérten |} 
AuBenstrémung iiberein, da alle Stromlinien parallel sind. An der Schlitzvorderkante entsteht ein Tot- | 
wasser mit Unterdruck (Abb. 1), an der Schlitzhinterkante ein Staugebiet mit Uberdruck. Bei der 
Mittelung des Drucks heben sich beide Stérungen teilweise auf. Die verbleibende Differenz ist meist es 
positiv und erzeugt bei kurzen Schlitzen einen entsprechenden Druckunterschied. Die Hafthedingung | 
an den Schlitzseitenwanden 148t nun nur einen Stérbereich in GréBenordnung der Schlitzbreite zu. 
Die Stérungen an den Schlitzenden klingen daher in einem Bereich ab, der der Schlitzbreite b propor- 
tionalist. Je langer nun der Schlitz im Verhaltnis zu seiner Breite ist, um so mehr iiberwiegt die fehler- 
freie Druckiibertragung im mittleren Teil gegentiber den stérenden Einfliissen an den Enden. Es ist 
daher zu erwarten, da die Differenz zwischen dem AuBendruck in Schlitzmitte und dem im Schlitz 
angezeigtem Druck, die von der Strémung im Schlitzinneren herriihrt, mit wachsendem Verhaltnis I/b 
(Schlitzlange zu Schlitzbreite) immer kleiner wird. 


ret Pinte Siege 
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3. Verhalten bei Druckgefiille. Zur Vereinfachung der folgenden Uberlegungen soll die Strémung 
im Schlitz durch die Annahme einer Trennungsflache idealisiert werden. Wie in Abb. 2 angedeutet, 
ersetzen wir die tatsachliche Geschwindigkeitsverteilung durch eine, die einen Knick hat. Die Tren- 
nungstlache ist dann der geometrische Ort dieses Knicks und trennt die ruhenden und bewegten 
Flissigkeitsteilchen innerhalb und auBerhalb des Schlitzes. 

Unter der Wirkung der verschiedenen Aufendriicke wird sich diese Trennungsflache so verwolben, 
daB sich durch die Zentrifugalkrafte in den gekriimmten Stromlinien der AuBenstrémung ein stetiger 
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Ubergang vom variablen Druck der AuBenstrémung auf den konstanten Druck unter der Trennungs- 
flache einstellt (Abb. 3). Uberall, wo der AuBendruck p(x) héher ist als der konstante Innendruck p,, 
wird die Trennungsflache nach innen gedriickt und erscheint, von auBen gesehen, konvex gekriimmt. 
Dort, wo der AuBendruck kleiner als der Innendruck ist, wird die Trennungsflache konkav gekriimmt. 
. Sa bezeichnet die Stelle x =x,, an der Innen- und AuBendruck iibereinstimmen 
Pi = P\*w))- 

Zunachst werde eine Strémung mit Druckabfall in Strémungsrichtung betrachtet. An der Vorder- 
kante des Schlitzes beginnt die Trennungsflache und kriimmt sich konvex in den Schlitz hinein. Ent: 
sprache der Innendruck dem AuBendruck in Schlitzmitte, so lage der Wendepunkt in Schlitzmitte, und 
die Trennungsflache wiirde unterhalb der Wandebene enden. Die auBen vorbeistrémende F liissigkeit 
staut sich an der Hinterkante. Eine entsprechende Druckerhéhung des Innendrucks Pp, der ja dem 
Druckmittel entspricht, verschiebt den Wendepunkt in das Gebiet héheren Druckes, also in Richtung 
auf die Vorderkante. Dadurch nahert sich das hintere Ende 
der Trennungsflache der Wandebene und der Aufstau am 
Schlitzende wird Kleiner. Es wird sich also ein solcher Innen- 
druck einstellen, der einem AuBendruck stromaufwarts der 
Schlitzmitte entspricht. Zwischen der Lage x,, dieses Gleich- 
gewichtspunktes und der Druckdifferenz Ap zwischen dem 
urspriinglichen Druck py in Schlitzmitte und dem angezeigten 
Druck p; besteht die Bezeichnung (Abb. 3) 


(7) 
Ap = Py — Pi = —%w 5 - (1) 


Bei einer Strémung mit Druckanstieg in Stroémungsrichtung 
(dp/dx > 0) wiirde sich auf Grund der gleichen Uberlegungen 
eine Trennungsflache ergeben, die etwa spiegelbildlich zur iisi.d. Veewuiliang oinier: Trasamioetliouh 
ersten verlauft. Bei einem Innendruck, der dem AuBendruck bei Druckabfall. 

in Schlitzmitte entspricht, endet die Trennungsflache oberhalb 

der Wandebene. Ks ergibt sich eine stérende Saugwirkung. Der Wendepunkt wandert in das Gebiet 
kleineren Druckes, erscheint also wieder in Richtung auf die Vorderkante verschoben. Der Faktor x,, 
in Gleichung (1) ist stets negativ. Andererseits kann man sich nur sehr schwer einen Mechanismus vor- 
stellen, bei dem der Gleichgewichtspunkt vor der Schlitzvorderkante liegt. Fiir die Lage des Druck- 
gleichgewichtes ergibt sich demnach folgende Abschatzung: ; 


1 Lap 


4, Versuchseinrichtung. Die Versuche beziehen sich auf die Druckentnahme aus einer laminaren 
Grenzschicht. Als AuBenstrémung wurde die Strémung in einem langen, flachen, rechteckigen Kanal 
(Querschnitt 18 x 170mm), Abb. 4, gewahlt. Das Geschwindigkeitsprofil der ungestérten Aufen- 
strémung entspricht der voll ausgebildeten laminaren Strémung zwischen zwei unendlich breiten, 
parallelen Wanden. Die Dicke einer vergleichbaren Grenzschicht entsprache gréSenordnungsmabig 
einer halben Kanalhéhe, also etwa 9 mm. Als Versuchsflissigkeit wurden zwei verschiedene Schmieréle 
verwandt, deren kinematische Zahigkeit in Abhangigkeit von der Temperatur in Abb. 5 dargestellt 
wurde. Mit Hilfe der Erwarmung des Ols beim Betrieb und durch Drosselung der DurchfluBmenge war 
es modglich, die auf den gleichwertigen Durchmesser und auf die mittlere DurchfluBgeschwindigkeit u 
bezogene Reynoldsche Zahl Re = ud,/y der Kanal-Strémung in den Grenzen 


50 < Re < 2200 
zu variieren. - : se 
Der Druckgradient der Kanalstrémung ist immer negativ. Um auch Strémungen mit positivem 
Druckgradienten verwirklichen zu kénnen, wurden noch zwei schlanke Diffusoreneinsatze mit den 
Erweiterungswinkeln 3°28’ und 5°14’ verwandt (Abb. 4). é 
Die untersuchten Schlitze waren geometrisch ahnlich. Ihre Form ist aus Abb. 6 und ihre Anordnung 
im Boden des Kanals aus Abb. 7 zu ersehen. iS ' 
Als MaB der Schlitze mige die Schlitzbreite b dienen. Sie variiert von 2 bis 10 mm und ist damit 
Kleiner oder gerade so grof wie die halbe Kanalhéhe h/2 = 9 mm, die die Grenzschichtdicke vertritt. 
Die in Abb. 7 ebenfalls eingezeichneten drei Vergleichslécher dienen zur Bestimmung des statischen 
‘Druckes im Kanalquerschnitt. Samtliche Druckanbohrungen im Kanalboden sind so weit von den 
6* 
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Abb. 6. Schlitzabmessungen. Abb. 7. Anordnung der Mefischlitze. 
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Kanalseitenwanden entfernt, daB die Anderung der Wandschubspannung im Kanalquerschnitt ver- 
nachlassigbar ist. 

Die im ersten Abschnitt geforderte weite und lange Schlitzvorkammer wurde durch Beruhigungs- 
kammern verwirklicht, wie sie im Schnitt in Abb. 4 gezeigt sind. Die Kammern sind tiefer als der gréBte 
Schlitz lang ist und enthalten zusatzlich noch Bleche, um die Riickstrémung in der Kammer mit 


Sicherheit durch Reibung abbremsen zu kénnen. Erst durch diese MaBnahme wurden die Messungen 
reproduzierbar. 


5. Messungen. Gemessen wurde: 


1. die Druckdifferenz zwischen einem Schlitz von der Breite b und dem seitlich daneben in Schlitz- 
mitte liegenden Vergleichsloch 


Ap* = PVergleichsloch a Pi > (3) 
2. der Druckgradient dp/dx im glatten Kanal und im Falle der Diffusoreneinsatze die dem Druck- 


gradienten entsprechende Druckdifferenz Ap; zwischen den beiden Wandbohrungen an der Oberseite 
des Diffusors, die L = 5 cm voneinander entfernt. waren; 


3. die Temperatur der MeBfliissigkeit T mit einem Thermometer stromabwarts der Schlitze; 


4. die Dichte @ [kpsec?/m*] der MeBfliissigkeit mit einem Araometer im groBen Sammelbehalter. 


Bei der Bestimmung der gesuchten Druckdifferenz Ap = py — p; besteht eine Schwierigkeit zu- 
nachst darin, da der Druckentnahmefehler des Vergleichsloches nicht bekannt ist. Um diesen Fehler 
zu bestimmen, wurde die Druckdifferenz Ap” fiir alle 5 Schlitze méglichst gleichzeitig gemessen und 
uber der Schlitzbreite aufgetragen (Abb. 8a). Der Fehler des Vergleichsloches wurde dann durch 
lineare Extrapolation auf die SchlitzgréBe b = 0 ermittelt. Um den so ermittelten Fehler zu priifen, 
wurde die Messung dann bei ungefahr der gleichen Reynoldschen Zahl, aber anderer Zahigkeit und 
anderem Staudruck wiederholt, und der zu dem neuen Staudruck gehérende Vergleichslochfehler wieder 
durch die gleiche lineare Extrapolation ermittelt. Mit Hilfe der beiden so gewonnenen Vergleichsloch- 
fehler wurde aus der gemessenen Druckdifferenz Ap* dann die gesuchte Druckdifferenz Ap berechnet, 
und das-Ergebnis beider Eliminationsverfahren dimensionslos (Re Ap/q) aufgetragen (Abb. 8b). Bei 
dieser Darstellung wichen die MeBpunkte der beiden Verfahren nicht allzu viel voneinander ab. Dabei 

__ ist der Unterschied, wie sich spater zeigen wird, zum Teil auf die doch etwas verschiedenen Re- Werte 
_ guriickzufihren. Damit erscheint die lineare Extrapolation in erster Naherung berechtigt. 


L— 1.Messung 
4 2.Messung 


Abb. 8. Lineare Extrapolation. a) Direkte Messung; b) Dimensionslose Darstellung. 


Die Dichten der beiden Olsorten unterscheiden sich kaum voneinander und waren fast temperatur- 
unabhangig. Es geniigte, mit der mittleren Dichte 9 = 89,0 kp sec?/m4 zu rechnen. Die kinematische 
Zahigkeit y ist, der konstanten Dichte wegen, eine Funktion der Temperatur T allein, so daB nach einer 
einmaligen Eichung (Abb. 5) die Messung der Oltemperatur eine Zahigkeitsmessung ersetzt. 
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Der Druckgradient dp/dx wurde im glatten Kanal mit Hilfe von Wandanbohrungen direkt ge-_ || 


messen. Die mittlere DurchfluBgeschwindigkeit im Kanal u ergibt sich dann aus dem Widerstands- 


gesetz der laminaren Rohrstrémung’: 


y 
Dabei ist x das Kanalseitenverhaltnis aka 
X= Kecatbroite 7 (7106 » fe 
d, bezeichnet den ,,gleichwertigen™ Durchmesser 
dg 4 Rell = 335 .10-2m. : 
Mit u kann man dann eine mittlere Reynoldsche Zahl 
Re = "oe! (7) 
und einen mittleren Staudruck 
as () 


ausrechnen. : 
Im Falle der Diffusorstrémung konnten die entsprechenden Mittelwerte leider nur abgeschatzt 
werden. : 
An Stelle der wirklichen Druckgradienten im Diffusor wurde die Druckdifferenz A Pr, die an der 
Oberseite des Diffusoreinsatzes gemessen wurde, mit dem Abstand L dieser MeBbohrung diffidiert. Es 
ergibt sich ein Naherungswert 


op Ae Apr d 
Sas (°) 


OxXp. 


fir den Druckgradienten, der wahrscheinlich kleiner ist als der wirkliche Druckgradient. 


Die mittlere Geschwindigkeit Zp, bezieht sich auf die Diffusorhéhe hp; des Diffusorquerschnitts, 
der in der Mittelebene der fiinf MeBschlitze liegt. 
Ebene Strémung vorausgesetzt ergibt sich 


u 


By. (" ae mittl. Diffasorhbhe fe fir Diffusorerweiterung 3°28’, 
= ji — = 


Kanalhohe ~~ |1,21 fiir Diffusorerweiterung 5°14’. 
Mit der so definierten mittleren Geschwindigkeit ist auch der entsprechende Staudruck gp; = + 


gegeben. 
Die Reynoldschen Zahlen Rep; im mittleren Diffusorquerschnitt und Re im stromaufwarts ge- 
legenen Kanalquerschnitt stimmen in ebener Strémung aus Kontinuitatsgriinden iiberein: 


Re = Rep; (11) 


Damit sind die gesuchten Mittelwerte der Diffusorstroémung auf die bekannten Gréfen eines stromauf- 
warts liegenden Kanalquerschnitts zuriickgefiihrt. 

Anzumerken ist noch, daf die Messung der Druckdifferenz Ap* durch Schwankungen des statischen 
Druckes stark gestért werden kann. Die gleiche Schwankungsamplitude verandert zeitlich die Druck- 
anzeige des Vergleichsloches in anderer Weise als die der viel gréBeren Schlitze. 

Als Ursache dieser Druckschwankungen stellten sich mitgerissene Luftblascheu und der DurchfluB 
sich ausbreitender Turbulenzgebiete heraus. 2 

Mitgerissene Luftblaschen anderten die Férdermenge der benutzten Schleuderpumpe erheblich. Zu 
ihrer Beseitigung wurde (Abb. 4) der Sammelbehalter B unstetig verengt, um geniigend Druck im Ein- 
laufrohr G zur Pumpe A zu haben. Weiterhin wurde der Uberlauf C so unterteilt, daG ein weiteres 
Beruhigungsbecken fiir das zulaufende Ol entstand und ein drittes Beruhigungsgefa8 H wurde vor dem 
Einlauf der Mefstrecke eingeschaltet. AuSerdem wurde im Uberlaufbecken C das Einlaufrohr zur 


1 A. Betz, Mechanik unelastischer Fliissigkeiten. Hiitte Bd. 1, 28. Aufl., S. 782. 
® J. Rotta, Ing.-Arch. 24 (1956), S. 285, 


| 
| 
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MeBstrecke mit groben Kieselsteinen iiberdeckt, damit hier keine Luftblaschen angesaugt werden 
kénnen. Infolge dieser MaBnahmen entwichen die beim Anstellen des Kanals mitgerissenen Luft- 
blaschen in etwa 4/, Stunde. 

Die sich ausbreitenden Turbulenzpropfen kiénnen nicht vermieden werden und lassen die bei 
eréBeren Re,-Werten gemessenen Druckunterschiede etwas fragwiirdig erscheinen. 


6. Darstellung und Diskussion der Messungen. Um die MeBergebnisse gut tibersehen zu kénnen, ist 
es verniinftig, sie dimensionslos darzustellen. Dabei verwendet man zweckmaBig solche Faktoren, 
welche die wichtigsten EinfluBgréfen beriicksichtigen. Die GréBe Ap, fiir die wir uns interessieren, 
hangt offenbar ab von dem gemessenen Druckabfall J dp/dx langs eines Schlitzes und der Reynoldschen 
Zahl Re = wd,j/y. Méglichwerweise hat auch das Verhaltnis J/d,, der Schlitzlange zur Kanalhéhe einen 
Einflu8, besonders bei den groBen Schlitzen. Es wird also 


So ee | 
Ap ae ; Re; z,] (12) 
sein. 
Die GréBe Ap wird am besten mit der Druckanderung iiber dem Schlitz (dp/dx) 1 dimensionslos 


gemacht. Diese Darstellung hat den Vorteil, da8 man den Druckgleichgewichtspunkt x,, unmittelbar 
erkennen kann. Nach (1) und Abb. 3 gilt ja fiir die Gleichgewichtslage x,, 


Ap Loe 
op ie F 
ox! 


Der Druckgradient dp/dx soll so dimensionslos gemacht werden, daB der entsprechende Parameter 
nicht mehr von der Reynoldschen Zahl Re =ud,/y abhangt. Dazu wird das Widerstandsgesetz 
4 =A,/Re fiir laminare Strémungen benutzt. Entsprechend Gleichung (4) ist demnach 


ép ho q 


éx Re da ° 
Macht man dp/dx mit 1/Re q/d,, dimensionslos, so erhalt man 
r7) 
we ae as Aes (13) 


In einer glatten Kanalstrémung ist A) eine nur vom Kanalquerschnitt abhangige Konstante. In einer 

laminaren anliegenden Diffusorstrémung wird J, auBerdem auch noch durch den Diffusorerweiterungs- 

winkel « und durch den Quotienten aus DurchfluBmenge Q und kinematischer Zahigkeit y mitbestimmt. 
Fiir die dimensionslose Darstellung der MeBergebnisse empfiehlt sich also die Form 


ay = — 5 =f{tes Res 7) | (14) 


Da im glatten Kanal J, konstant ist, bleibt nur eine Abhangigkeit von Re und 1/d,, iibrig. Dementspre- 


- chend ist in Abb. 9 Ap/(0l/dx) als einparametrige Kurvenschar abhangig von Re mit dem Parameter 


Id, aufgetragen. Um den Einflu8 dieses Parameters tibersehen zu kénnen sind jeweils die Geraden 
Ap|l (dp/dx) = 0,3 in Ordinatenrichtung um einen Betrag versetzt, der der Schlitzlange / direkt pro- 
portional ist. - 

Man erkennt unmittelbar, daB die Gleichgewichtslage x,,/1 von der Reynoldschen Zahl Re unab- 
hangig ist, solange die mit der Schlitzbreite gebildete Reynoldsche Zahl den Ubergangswert Re, = 
bu/d,,; = 60 nicht iiberschreitet. Innerhalb dieses Re-Bereiches ist — Xyoll auch von der dimensions- 
losen Schlitzlange l/d,, unabhangig, wenn die SchlitzgréBe I/d,, < 1,1 bleibt. In diesem Bereich ergibt 
sich 


Ap _ we 15 
p= — f= 0,39 . (15) 
oa 


Der Gleichgewichtspunkt zwischen angezeigtem Druck und Druck der AuBenstrémung ist bei allen 
Schlitzen um 39% der Schlitzlange von der Mitte in Richtung auf die Vorderkante verschoben. Fiir 
Schlitze, bei denen Re, = ub/y > 60, wandert der Gleichgewichtspunkt langsam zur Schlitzmitte. 
Die Messungen in den Diffusoreinsatzen kénnen in der dimensionslosen Form (14) nicht recht dar- 
gestellt werden. Da sie namlich vor allem im Hinblick auf Strémungen mit vermindertem Druckgra- 
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dienten durchgefiihrt wurden, ist die BezugsgréBe (dp/dx) I, die beim glatten Rohr zur dimensionslosen 
Darstellung der gemessenen Ap-Werte benutzt wurde, oft so klein, daB sie infolge der MeBungenauig- 
keit erheblich verfalscht werden kann. Als Grundlage fiir die Darstellung der Diffusormessungen 
wurde daher die dimensionsbehaftete Form (12) zugrunde gelegt (Abb. 10). 


5 ees Ta\ Fak SS 
Bae Ufliyy=185 
» eS ES es 
Verbindungsiinie aller Abszissenpunkte 
fiir die Re,=Re - 2 =69 
dy 
OM es aie 
03 - \ ENE = —< 
02 \ Yadgy = 411 
ony x 
H-O4 = iB fs 
2 Ud, =0,740 
03 ne Oo, fe ge 7 
gutoze 
dips : 9 Udy, =0555 
93 << 2 ——— - 
02}-O4 = — — a Yatyy = 0,369 
03 Re 
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Abb. 9. MeBergebnisse beim glatten Kanal. 
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Abb. 10. MeBergebnisse beim Kanal und beim Diffusor. 
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Bei allen MeBpunkten war Re, = u b/y < 60, und I/d ie Pai 

3ei al jar =< 60, 1 < 1, so daB die Parameter Re, I/d - 
aussichtlich wieder keinen EinfluB haben. Als Gnabhensice Veranderliche hatte man ee 
die Druckanderung iiberm Schlitz (dp/dx) 1 wahlen kénnen, als abhangige Veranderliche Ap. Die 
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Gleichgewichtslage — x,,/l ergabe sich dann als Tangens des Richtungswinkels, der zu dem entsprechen- 
den Me8punkt gehért. Um die Darstellung iibersichtlicher zu machen, und um wieder dimensionslose 
Werte zu erhalten, wurden Ordinate und Abzisse noch jeweils mit dem zugehérigen Re/q-Wert multi- 
pliziert. Die Tangensbeziehung wird dadurch nicht gestért. AuBer den MeBpunkten der Diffusorstré- 
mung wurden auch noch einige MeBpunkte der Kanalstrémung eingetragen. Sie sind durch ,,0° 
Diffusorerweiterung“ gekennzeichnet. 


Man kann dieser Darstellung folgende Ergebnisse entnehmen: 


a) Die bei verschiedenen Re-Zahlen und an verschiedenen groBen Schlitzen gewonnenen MeBpunkte 
fallen unmittelbar zusammen. Auch fiir die Druckentnahme aus einer Diffusorstrémung gilt, daB das 
Ergebnis von der Reynoldschen Zahl Re und der dimensionslosen Schlitzlange I/d,, unabhangig ist, so- 
lange Reb kleiner als etwa 60 bleibt, und Uda < 1-1 ist. i 


b) Die MeBpunkte lassen sich durch Kurven annahern, die innerhalb der MeBgenauigkeit durch den 


Nullpunkt gehen. Die Druckdifferenz Ap ist also, wie erwartet, in einer Strémung ohne Druckgradient 
sehr klein. 


c) Die positive Steigung der Naherungskurven bedeutet, daB der Punkt x,, des Druckgleichgewichts 
stromaufwarts der Schlitzmitte liegt, im Einklang mit den theoretischen Uberlegungen. 


d) Die Naherungskurven sind Gerade durch den Nullpunkt. Das bedeutet, daB die Lage des 
Druckgleichgewichts im Bereich der gemessenen Druckgradienten von diesen unabhangig ist. Die 
Linearitat bleibt auch dann noch erhalten, wenn der Druckanstieg bis fast zur Ablésung und Riick- 
strémung gesteigert wurde. 


_e) Die Steigung der Geraden, die die in den beiden Diffusoreinsatzen bestimmten MeBpunkte an- 
nahert, betragt rd. 0,75. Das wiirde bedeuten, da der Druckgleichgewichtspunkt 25° der Schlitz- 
lange stromaufwarts vor der Vorderkante liegt. Das ist sehr unwahrscheinlich. Der gemessene Druck- 
gradient Ap,/L ist vermutlich zu klein. Wahrscheinlich liegt das daran, daB die vordere Druckanboh- 
rung S, (Abb. 4) aus Raumgriinden zu weit nach vorn gelegt wurde. Die Druckverteilung lings des 
Diffusoreinsatzes besitzt ahnlich wie ein endlich dickes symmetrisches Profil in Nahe der Vorderkante 
eine ausgesprochene Saugspitze. Die Tatsache, daB der gemessene Druckgradient wesentlich zu klein 
ist, deutet darauf hin, daB die vordere Bohrung S, den Druck stromaufwarts der Saugspitze gemessen 
hat. Bei korrekter Durchfiihrung der Diffusormessungen wiirde sich wahrscheinlich dieselbe Lage des 
Druckgleichgewichts ergeben; wie sie im glatten Kanal gemessen wurde. 


7. Zusammenfassung. Ziel der vorliegenden Arbeit war, die Fehler zu untersuchen, die bei der 
Druckentnahme durch einen langen Schlitz durch die Sekundarstrémung im und in Nahe des Schlitzes 
hervorgerufen werden. Insbesondere sollte die theoretisch gewonnene Vermutung experimentell ge- 
priift werden, daB der angezeigte Druck dem Druck der urspriinglichen Strémung an einer Stelle ent- 
spricht, die zwischen Schlitzvorderkante und Schlitzmitte liegt. 

Die Druckdifferenz Ap zwischen AuSendruck in Schlitzmitte und angezeigtem Druck wurde fiir 5 
geometrisch ahnliche Schlitze (//b = 6, 1/b = 1) in einer laminaren Kanal- und Diffusorstrémung 
gemessen. Die qualitativen theoretischen Uberlegungen werden bestatigt, der Druckunterschied hangt 
in erster Linie von der Druckanderung iiber dem Schlitz ab und verschwindet in einer Str6mung ohne 
Druckgradient. Das Druckgleichgewicht befindet sich an einer Stelle, die in gewissen Grenzen unab- 
hangig vom Druckgradienten 11% der Schlitzlange hinter der Vorderkante liegt. Voraussetzung ist 
allerdings, daB die Riickstrémung im Schlitzinnern durch Reibung gentigend gebremst wird und eine 


gewisse Reynoldsche Zahl der Schlitzbreite 
b 


=! 


~ 60 


=| 


nicht tiberschritten wird. 


Hiermit danke ich Herrn Professor Dr. Betz herzlich fiir die Anregung und stete Forderung 
der Arbeit. Weiterhin danke ich den Mitarbeitern des Max-Planck-Instituts fiir Strémungsforschung 
und der Aerodynamischen Versuchsanstalt fiir ihren hilfreichen Rat bei der Durchfihrung. AuBer- 
dem michte ich der BV-ARAL Aktiengesellschaft fiir die Spende von 550 kg Hochleistungsél 
danken, das fiir die Versuche besonders geeignet war. 


(Eingegangen am 9. Mai 1959.) 
Anschrift des Verfassers: Dipl.-Phys. Fritz Krause, Gottingen, Bunsenstr. 10. 
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Erginzungen zum Grenzschichtquadraturverfahren von E. Truckenbrodt * 


Von N. Scholz 


1. Einleitung. In der Entwicklung der Berechnungsverfahren fiir die Grenzschicht kann man ins- 
besondere im Hinblick auf turbulente Grenzschichten zwei Richtungen unterscheiden: Einerseits 
wird durch eine Verfeinerung des als Grundlage dienenden physikalischen Inhaltes der Ausgangs- 
gleichungen eine bessere Ubereinstimmung der Rechenergebnisse mit der Wirklichkeit angestrebt. 
Andererseits soll eine méglichst weitgehende Vereinfachung und Vereinheitlichung des Rechen- 
ganges erreicht werden, so da diese Verfahren nicht nur von Spezialisten gehandhabt, sondern in 
der ingenieurmaBigen Praxis von Hilfskraften oder Rechenanlagen mit lohnendem Aufwand ausge- 
fiihrt werden kénnen. Wegen der einheitlichen Behandlung laminarer und turbulenter Grenzschichten 
im ebenen und rotationssymmetrischen Fall erscheint das von E. Truckenbrodt' entwickelte Quadratur- 
verfahren im letzteren Sinne besonders geeignet und wurde deshalb auch bereits fir umfangreiche 
systematische Untersuchungen herangezogen” °. 

In der vorliegenden Arbeit werden fiir dieses Verfahren einige Erganzungen mitgeteilt, die zum 
einen den Anwendungsbereich erweitern, zum anderen die Treffsicherheit der Ergebnisse verbessern 
und die numerische Durchfiihrung der Rechnung vereinfachen sollen. Die Kenntnis der Trucken- 
brodtschen Arbeit wird daher vorausgesetzt. Im folgenden werden wir den ebenen und rotations- 
symmetrischen Fall gemeinsam behandeln. Fiir den ebenen Fall ist dabei in allen Gleichungen der 
Wandradius in der Ebene senkrecht zur Strémung R bzw. R/l = 1 zu setzen. 


2. Bezeichnungen. Es werden folgende Bezeichnungen verwendet: 


x und y ' Koordinaten in Richtung der Wandbegrenzung und senkrecht dazu, 

1 Bezugslange, im allgemeinen Lange der betrachteten Grenzschicht, 

R(x Wandradius bei rotationssymmetrischen Begrenzungen, 

U(x), o(x), T(x) Geschwindigkeit, Dichte und Temperatur auBerhalb der Grenzschicht, 

u(y), e(y), T(y) Geschwindigkeit, Dichte und Temperatur innerhalb eines Grenzschichtprofiles, 

T), Ty, T,, Tz Temperaturen des Strémungsmittels bei isentropem Aufstau, an der Wand, an der Wand 

Es infolge Reibungswarme und Bezugstemperatur nach Gleichung (5), 

0, O*, &, 6 Grenzschichtdicke, Verdrangungsdicke, Impulsverlustdicke und Energieverlustdicke des 
= Stromungsgrenzschichtprofiles, 

Ops Op Dicke und Warmegewinndicke des Temperaturgrenzschichtprofiles, 

k, KorngréBe der der Wandrauhigkeit aquivalenten Sandrauhigkeit, 

of Dimensionsloser Reibungswiderstandsbeiwert der Gleichdruckgrenzschicht (ebene Platte), 

cf _ _ Dimensionsloser Schubspannungsbeiwert, 

= 07/05 6/@ Grenzschichtdickenverhaltnisse, 

Hyp= 67/6 VerhaltnisgréBe zwischen Temperatur- und Strémungsgrenzschichtprofil, 

L Formparameter des Grenzschichtprofiles nach Gleichung (18), 

n, tt Exponenten des Dissipations- und des Widerstandsgesetzes, 

w Exponent des Zahigkeitsgesetzes als Funktion der absoluten Temperatur, 

x Isentropenexponent, 

a, a,b, d,e, A, B Konstanten, 

C,k, p,q, m,c - Exponenten, 

@, A,é = Oc Dimensionslose GréBen, definiert durch Gleichung (4) bzw. (23), 

Reg = U By értliche Reynoldszahl des Grenzschichtprofiles, 

Re, = Ug |/v Bezugs-Reynoldszahl, 

Pr Prandtlzahl, 

M Machzahl, - 

r Erhitzungsfaktor, definiert durch Gleichung (6). 
Als Indizes stehen: 

20 fur GréBen der Bezugsstrémung oder bei einer Grenzschicht mit Gleichdruck, 

0 fiir Groen des Anfangszustandes der Grenzschicht, 


* Uberarbeitete Fassung eines Berichtes des Institutes fiir Strémungsmechanik der Technischen Hoch- 
schule Braunschweig (Prof. Dr. H. Schlichting). 

1 E. Truckenbrodt, Ing.-Arch. 20 (1952), S. 211. 

* L. Speidel u. N. Scholz, Untersuchungen iiber die Strémungsverluste in eb S lgi -For- 
AaeeGane steno) g gsverluste in ebenen Schaufelgittern, VDI-For 


* J. Weissinger, Mitt. Nr.9 d. Inst. f. angew. Math. d. T. H. Karlsruhe (1958). 
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/ eee fiir GréBen bei konstanten Stoffeigenschaften (zur Unterschei g liber 
chenden Gré8en bei er andexlichion. Se hen} pie ities entepre 

fiir mit der Integrationsvariablen x’ veranderlichen Gréfen, 

A fiir Werte am Ablésungspunkt der Grenzschicht, 

ry fiir Werte am v-ten Rechenpunkt. 


3. Berechnung der Impulsverlustdicke. a) Kinflu8 des értlichen Formparameters: Aus der 
Integration des Energiesatzes der Grenzschicht erhalt E. Truckenbrodt eine geschlossene Lisung fiir die 


Impulsverlustdicke, indem er fiir den ortlichen Formparameter H(x) den konstanten Wert der Gleich- 
druckgrenzschicht einsetzt. Ohne diese einschrankende Annahme kann die Gleichung fiir die Impuls- 
verlustdicke #(x) geschrieben werden!: 


Foy ee ey 
Co + aa Flee | U(x yt" REx) +” de’ 
Ud\ Bains 
ep 
Hierbei ist C, eine Integrationskonstante fiir die Impulsverlustdicke bei « = X. Die zu dieser Glei- 
chung fiihrenden Voraussetzungen, da8 namlich die dimensionslose Dissipation nur eine Funktion 


der drtlichen Reynoldszahl Re,? und die Gré8e E unabhangig von H und damit von der Integrations- 
variablen x’* sind, kénnen als weitgehend erfiillt angesehen werden*. Schreibt man Gleichung (1) 


fiir den Gleichdruckfall auf, so ergibt sich sofort die Bedeutung der GroBe E/H1+™, namlich wegen 


Buf = efoof2 . 
= See Oy wie 


Ure R(x)!” : (1) 


BLA NEP A Ne : ] 
Hiermit erhalt man eine fiir praktische Rechnungen 
bequeme Gleichung fiir die Impulsverlustdicke in 
der Form® | 
1 
(x) ti ag 7 UNO. ay Of os (ge g 
eee le UG RE) 
(ae) 7 
mit 
On — [2% Hy (U,\2 Rt +" if 
Scare 0.) 0 Of 02 Q3 g4 05 06 


x Abb. 1. Impulsverlustdicke der turbulenten Reibungsschicht am 
U(x’) 3+2n R(x’) l+tn , Profil NACA 65(216)—222. « =10,1°, Rel = 2,6 102, Messung 
+ ( d ag (4) nach NACA-Rep. 772, — — — Rechnung nach Truckenbrodt, 
co 
Xo 


l 


U l — .— Rechnung nach Gleichung (3). 


Gegeniiber Truckenbrodt enthalt Gleichung (3) das Glied H,,/H(x), das den Einflu8 des értlichen 


Formparameters auf die Impulsverlustdicke darstellt und im Druckanstieggebiet wegen H, Ts ec | 
den Impulsverlust erhéht. Abb. 1 zeigt fiir das Rechenbeispiel von Truckenbrodt (dort Abb. 15) die 
Verbesserung, die durch Gleichung (3) erzielt wird’. Bei der Berechnung von # ist zunachst die 
QuadraturgréBe O(x) zu ermitteln. Mit dieser kann, wie im folgenden noch gezeigt wird, die Qua- 


dratur des Formparameters erfolgen, und anschlieBend mit H(x) die Impulsverlustdicke erhalten 
werden. Damit ist keine Koppelung der beiden zu lésenden Gleichungen fiir die Impulsverlustdicke 
und den Formparameter vorhanden, so daf ein iterativer Lésungsgang nicht erforderlich ist. 


1 Siehe FuBnote 1 von Seite 82; dort Gleichung (34). 

2 Siehe FuGBnote 1 von Seite 82, dort Abb. 2 im laminaren und Abb. 4 im turbulenten Fall. 

3 Siehe FuBnote 1 von Seite 82, dort Abb. 6. 

4 Nur im laminaren Fallist E im Druckabfallbereich etwas von H abhangig. om ef 

5 Hierbei ist zur Vereinheitlichung des Rechenganges in der Integrationskonstanten H(x)~ Hy gesetzt 
worden. Diese Vereinfachung ist einerseits fiir Werte von x nahe bei x, gerechtfertigt, andererseits spielt fur 
grdBere Lauflangen x — x, die Integrationskonstante gegeniiber dem Integralwert eine immer untergeordnetere 
Rolle. 

6 Die von Truckenbrodt (Abb. 14) erzielte gute Ubereinstimmung fiir die laminare Grenzschicht der Howarth 


Strémung wird auch mit Gleichung (3) nicht schlechter, da hierbei H/H,, zwischen 1,00 und 1,03 liegt. 
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b) Einflu® veranderlicher Stoffeigenschaften!: Bei einem Abweichen der értlichen 
Temperatur innerhalb der Reibungsschicht von der Temperatur T der AuBenstrémung liegen ent- 
weder ein Warmeiibergang zwischen Wand und Strémungsmittel (Heizen, Kithlen) oder eine durch 
die Reibungswarme entstandene Aufheizung der Grenzschicht oder auch beides gleichzeitig vor. 
Die Riickwirkungen dieser Temperaturunterschiede auf die Reibungsschicht sind durch die Ab- 
hangigkeit der Stoffeigenschaften (Dichte und Zahigkeit) des Str6mungsmittels von der Temperatur 
verursacht. Es liegt nahe und ist von verschiedenen Forschern vorgeschlagen worden, den EinfluB 
der unterschiedlichen Temperaturen innerhalb der Reibungsschicht dadurch zu beriicksichtigen, daB 
die als konstant angenommenen Stoffeigenschaften des Strémungsmittels bei einer bestimmten, von 
der Temperatur der AuBenstrémung abweichenden Temperatur Tz, die als Bezugstemperatur be- 
zeichnet wird, angenommen werden. Fiir diese Bezugstemperatur wird die Beziehung 


Tg = T + 4(T,— T)—6(T,— Ty) (5) 


nahegelegt, wonach diese sich aus der AuBentemperatur T, einem infolge Reibungswarme von der 
Temperaturdifferenz T,— T und einem infolge Warmeiibergang an der Wand von der Temperatur- 
differenz T,— Ty, abhangigen Glied additiv zusammensetzt. Die Reibungstemperatur T, hangt in 
bekannter Weise von der Machzahl ab, namlich 


Ps Poe Te 


M2, se 


Der Erhitzungsfaktor r ist eine Funktion der Prandtlzahl. 

Wie der Verfasser! gezeigt hat, erweisen sich fir gasférmige Medien mit von der Temperatur 
praktisch unabhangiger Prandtlzahl die in Tabelle 1 zusammengestellten Konstanten fiir die Bezugs- 
temperatur TJ’, als brauchbar. 


Tabelle 1. 
Werte der Konstanten in Gleichung (5) und (6). 


Reibungsschicht 


laminar 
turbulent 


Eine Kinfiihrung der Bezugstemperatur in den von E. Truckenbrodt benutzten Ansatz fir die 
Dissipation des Grenzschichtprofiles erméglicht die Integration des Energiesatzes der Grenzschicht 
auch fiir den Fall veranderlicher Stoffeigenschaften!: 

1 


d(x ) FF ene) ak ae (Cfo) po O(x)i +” 
ap SS H —_ = ERS Se Ze 
I (x) I 7) o(x) ooh +2Hp R(x) (7) 
Qo \ UX ads 
mit 
Lf 2 Oh Ho 00 Uo 342Hp R lin 
O(x) (opsdiglt g SIDE 05 (e) “ 


I es U 7 d=. (8) 


ae ee 7( Ts. here | ae (Ben A 
Hierbei bringt der Formparameter Hy = 67/6 eine Koppelung zwischen der Strémungs- und Tem- 


peraturgrenzschicht. Sein genauer Wert ist jedoch fiir das Ergebnis (7) von untergeordneter Be- 
deutung. Es kann daher die Naherungshbeziehung! 


a Ae! Tyx) — T,(x) 
H,(x) == M%x) SS a (9) 
die iibrigens fiir warmeundurchlassige Wand (Ty = T,) exakt gilt, zur Auswertung von (8) herange- 
zogen werden. Dariiber hinaus ist es in vielen Fallen im Rahmen der Naherung ausreichend, fiir 


Hy einen konstanten Mittelwert langs der Grenzschichtlange einzusetzen. 


if Ei ee . 2 . . 
see 71059), 8. 38 Darstellung der Uberlegungen dieses Abschnittes hat der Verfasser gegeben: NV. Scholz, 
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Fir U = U, wird aus (7) der Reibungsbeiwert der Gleichdruckgrenzschicht (ebene Platte) mit 
veradnderlichen Stoffeigenschaften erhalten. Das Verhaltnis zum Reibungsbeiwert der Gleichdruck- 
grenzschicht mit konstanten Stoffeigenschaften ist daher 

1—nw 
Cf ou eeleg Ceren (10) 
(Cro) p oa, = (75) 
Abb. 2 gibt das zahlenmafige Ergebnis dieser Beziehung in Abhangigkeit von der langs der Grenz- 


schicht konstanten Wandtemperatur und der Machzahl der Anstrémung als Parameter bei turbu- 
lenter Grenzschicht wieder. 


ce) Einflu8 der Reynoldszahl und Wandrauhigkeit: Wie J. Rotta! fiir die turbulent 
Grenzschicht gezeigt hat, beschrankt sich dieser Einflu8 auf eine diinne wandnahe Schicht, deren 
Ausbildung im wesentlichen unabhangig vom iibrigen Ce. 
Grenzschichtverlauf ist. Man kann daher, wie es auch mes 
vom Verfasser? bereits durchgefiihrt worden ist, die zur 
Integration der Impuls- und Energiegleichung der Grenz- 
schicht erforderlichen Ansatze ebenso fiir hydraulisch 
rauhe Wandbeschaffenheit anschreiben, wobei lediglich 
das értliche RauhigkeitsmaB #/k, anstelle der drtlichen 
Reynoldszahl Re, bei hydraulisch glatter Wand und 
selbstverstandlich andere Zahlenwerte der Konstanten 
auftreten. Letztere sind schlieBlich in (3) bis auf den 
Exponenten n und den Formparameter H, in dem 
Reibungsbeiwert der Gleichdruckgrenzschicht c;,, zu- 
sammengefaBt, dessen Abhangigkeit von der Reynoldszahl 
bzw. Wandrauhigkeit wohlbekannt ist. 


Fir den Exponenten n 148t sich ein Zusammenhang 
mit dem Exponenten n des Reibungsgesetzes der Gleich- 
druckgrenzschicht 


(od al 0 OY O8 12 16 0 
foo a n ( ) Abb. 2. Verhaltnis der Widerstandsbeiwerte der ebenen 
Re, Plattengrenzschicht mit und ohne TemperatureinfluB nach 


Gl. (10) in Abhangigkeit von der Wandtemperatur. Stro- 


herstellen®, namlich mungsmittel: Luft mit Pr = 0,72, m = 0,76 und n = 0,19 


n (turbulent, Re; = 107), —-— Tp = TR (kein Warme- 
n= —— (12) ibergang), — — — Ty = Too (Kihlung auf Aufen- 
1+ n° temperatur). 


Bekanntlich nimmt n mit zunehmender Reynoldszahl ab. Es kann aus dem universellen log- 
arithmischen Reibungsgesetz nach Prandtl-Schlichting 


Cf. = 0,455 (log Re,)—?8 (13) 
durch Differentiation zu ae 
na log Re, (14) 


gewonnen werden (Abb. 3). Fiir laminare Grenzschicht ist bekanntlich n = 1/2. 

Der Formparameter H der turbulenten Gleichdruckgrenzschicht ist in Untersuchungen von 
J. Rotta4 als Funktion des Schubspannungsbeiwertes cr bzw. Res enthalten. Unter der vereinfachenden 
Annahme, da8 H,, langs einer Gleichdruckgrenzschicht im wesentlichen gleich dem Wert am Ende 
der Grenzschicht ist, kann jedem Wert von Re, mit Hilfe des Reibungsgesetzes (13) ein Wert von 
Re, zugeordnet werden, so das damit H..(Re;) gewonnen wird. In Abb. 3 ist diese Beziehung eben- 
falls eingezeichnet und mit Messungen von F. Schultz-Grunow® und W. Tillmann® verglichen, wobei 
sich eine recht brauchbare Ubereinstimmung ergibt. Im laminaren Fall ist H, wiederum konstant 
gleich einem Wert von etwa 2,6. 


1 J. Rotta, Ing.-Arch. 19 (1951), S. 31. 

e Ni Scholz, Jeiee. d. Scheib. Techn. Ges. 45 (1951), S. 244. 

3 Siehe z. B. FuGnote 1 von Seite 82, dort Gleichung (37a). - 

4 J. Rotta, Mitt. Nr. 8 d. MPI f. Strémungsforschung (1953), S. 1. 

5 F. Schultz-Grunow, Luftfahrtforschung 17 (1940), S. 239. ; 

6 W. Tillmann, Untersuchungen iiber Besonderheiten bei turbulenten Reibungsschichten an Platten, KWI f. 


Strémungsforschung, UM 6627 (1945). 
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Die Werte von 7 und H,, fir die turbulente Grenzschicht bei hydraulisch rauher Wand lassen 
sich nunmehr durch Zuordnung einer der relativen Wandglatte I/ks aquivalenten Reynoldszahl Re, 
bestimmen. Vernachlassigt man den fiir jede Rauhigkeitsstruktur speziellen Verlauf der Rauhig- 
keitsfunktion im Ubergangsgebiet ,,rauh—glatt“!, so setzt sich die Rauhigkeitsfunktion aus den 


sm Ais ie 
ort qt 028 : + | ak Baca 
| [ 
008+ 16+ O24 | + 
Goel 45+ O20 os | } t 1 
ok 4k 06 ee 
904+ 13+ ane re 


mec: 
PL = ected Re, 


7 8 I 1 
10 0 70 a4 


-0,16 10 0 


70° 70° 10" 


Abb. 3. KenngréBen der Gleichdruckgrenzschicht n, Hoo und Loo als Funktion der Reynoldszahl Re; bzw. der relativen Glatte I/k, fiir Hoo: 
© Messung ,,glatt‘t nach Rotta (FuBnote 4 von Seite 85), + Messung ,,rauh‘* nach Rotta (FuBnote 4 von Seite 85), fir Lao: nach Verfasser, 
—-+-— nach E. Truckenbrodt, — — — nach Tabelle 3. 


beiden den ,,glatten“ und ,,rauhen“ Bereich wiedergebenden Geraden zusammen, deren Schnittpunkt 
die Zuordnung der dquivalenten Reynoldszahl zu einer relativen Wandglatte bestimmt. Danach 
befindet man sich mit einer Grenzschicht dann im hydraulisch rauhen Bereich, wenn die zugehérige 
Wandschubspannung gréfer als die entsprechende fiir hydraulisch glatte Wand ist. Auf Grund der 
Berechnungen des Schubspannungsbeiwertes c¢ als Funktion der Gértlichen Reynoldszahl Reg bzw. 
der értlichen Wandglatte 3/k,; und des Formparameters H von J. Rotta laBt sich fiir die Grenze 
,rauh—glatt“ die Interpolationsformel 


==. rlatt:: 
a (15) 


log Us = 1,58 + 0,30 a} te 


angeben, die fiir H = 1,4 den als Faustformel bekannten Wert U k,/y = 100 liefert. Mit zuneh- 
mendem H steigt, wie auch anschaulich einleuchtet, die zulassige Rauhigkeit bei gleicher AuBen- 
geschwindigkeit U an. Setzt man fiir die Gleichdruckgrenzschicht in (15) H = H,, nach Abb. 3, so 
ergibt sich die in Abb. 3 als zweite Abszissenskala eingetragene Zuordnung von l/k, und Re,, womit 
auch n und H,, als Funktion von I/kg vorliegt. 

Mit der Einfithrung einer Aquivalenten Reynoldszahl fiir die hydraulisch rauh wirkende Wand 
ist auch das Reibungsgesetz der rauhen ebenen Platte festgelegt. Ersetzt man in dem Reibungs- 
gesetz (13) die Reynoldszahl mit (15) und H = H,, nach Abb. 3 durch die relative Glatte I/ks, so 
ergibt sich ein Reibungsgesetz, daB sich durch folgende Gleichung interpolieren laibt: 


Cf. = 0,418 [2 + log (I/kg)]->**. (16) 
Vergleicht man dieses Gesetz mit dem von L. Prandtl und H. Schlichting? (Abb. 4), so findet man 


besonders bei groBenRauhigkeiten niedrigere c;-Werte, was mit den Beobachtungen von W. Tillmann 
ibereinstimmt. Auf die Ursache dieses Unterschiedes ist bereits von J. Rotta? hingewiesen worden. 


* Diese Annahme wird auch von Rotta und vom Verfasser gemacht. Uber den Begriff der Rauhigkeitsfunktion 
und deren Verlauf vgl. z. B. H. Schlichting, Grenzschichttheorie S. 485ff., 3. Aufl. Karlsruhe 1958. 

® H. Schlichting, a.a.O. Gleichung (21, 33). 

3 Siehe FuBnote 4 von Seite 85. 
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Aus (16) ergibt sich durch Differentiation und Verelei i i 
< é gleich mit (11) eine 14 los : x 
fiir den Exponenten n bei hydraulisch rauher Wand, namlich ie au (14) analoge Bezichung 


me 1,09 
log (2 + Iks) ° (17) 


Wahrend der EinfluB der Reynoldszahl bzw. Wandrauhigkeit auf den Formparameter H,, und damit 
auf die Impulsverlustdicke nach (3) und (4) durchaus zu beriicksichtigen ist, kommt der Exponent n 
hierbei in einer solchen Kombination vor, daB seine Anderung mit der Reynoldszahl oder Wand- 
rauhigkeit auf das Ergebnis einen nur untergeordneten Einflu8 hat. 


Cr 

070 + z “ = 
0,009 5 
0,008 
9007, il 2 | z 
0,006 = 

> 
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— 
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Abb. 4. Reibungswiderstandsbeiwert der hydraulisch rauhen ebenen Platte nach Prandtl-Schlichting ( 
und nach (Gleichung 16) (— — —). 


) 


Tabelle 2 gibt eine Zusammenstellung dieser Werte, wobei bei der turbulenten Reibungsschicht 
die gréBeren Werte dem Blasiusschen Reibungsgesetz (n = 1/5) und die kleineren dem Falknerschen 
(n = 1/7) entsprechen. 


Tabelle 2. Werte der Exponenten in Gleichung (3) und (4). 


Reibungsschicht n | l+n | 34 2n 


1 | 2 5 
0:17-—0,25 | 1,7F=1,25 © |) 3.33350 


4, Quadraturformel fiir den Formparameter. a) Umformung der Quadraturformel. E. Truk- 
kenbrodt fiihrt fiir die laminare und turbulente Grenzschicht den neuen Formparameter 


de 3 OL (18) 
(eae 


ein und gewinnt aus einer Kombination von Impuls- und Energiesatz der Grenzschicht eine geschlossen 
integrierbare Differentialgleichung fiir L, deren Lésung schlieBlich lautet* 


0,9 
0,14— 0,20 


laminar 
turbulent 


é 
L() = #1 +ng + = [4 In Rep — B—In 7, | db". (19) 
Ey 


Hierbei sind A und B empirische Konstanten und = 9© mit @ nach (4) und C einer weiteren em- 
pirischen Konstanten?. 


1 Siehe FuGnote 1 von Seite 84; dort Gleichung (57). 

2 Nach Truckenbrodtist z. B.im turbulenten Fall A = 0,07 log e, B = 0,23 und C = 4,0. Ferner unterscheidet 
sich € um einen Faktor von der entsprechenden Gréfe bei Truckenbrodt, der sich jedoch in (19) wieder heraus- 
kiirzt. Hieraus ergeben sich gewisse Vorteile bei der numerischen Rechnung. 


sae 
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Zunachst kann in (19) die in Re, enthaltene Impulsverlustdicke # durch die GréBe @ nach (3) 


ersetzt werden. Dann schreibt sich Gleichung (19): 


E é a 
cpa aa tee — 
& 0 
E é 
+ pAy [ moa + [In gt + Aln(Ra%)—B je (20) 
& s 


Hierbei ist das dritte und vierte Integral unabhangig von der speziellen Geschwindigkeitsverteilung 
und kann integriert werden!. Bei dem ersten Integral ist es weiterhin zweckmabig, die Tntegrations- 
variable £’ in die bei der numerischen Rechnung zugrundeliegende Variable x‘/] tiberzufihren. Mit 
diesen Umformungen wird schlieBlich 


"2 I U A n@ A 21 

eet AOS In Us 1 ae 51 é 9 ( ) 
Hierbei bedeutet 

AA Cfo 22 

L =i + Aln(Re, a (22) 


den Formparameter der Gleichdruckgrenzschicht, da fiir U = konst und J) = 0 alle weiteren Glieder 
der Gleichung (21) verschwinden. In diesen tritt aber auch die Reynoldszahl nicht mehr auf, so dal 
L,, allein den Einflu8 der Reynoldszahl auf den Formparameter bestimmt. Damit gilt (21) auch 
fiir hydraulisch rauhe Wand, da lediglich fiir L,. ein entsprechender Wert einzufiihren ist, der sich 
ohne weiteres aus der zu der Wandrauhigkeit aquivalenten Reynoldszahl ergibt. Die GréBe A ist 
aus einer Quadratur langs der Wandkontur zu gewinnen: 


Upper. 
UU. Fen 


; H dO x” : 
—6[\a+24m ge +t (F mle (is) | (23) 
; T 


x 


A = f5(Lp— Lo —In Ino) 


Die numerische Auswertung des in (23) auftretenden Integrales gestaltet sich dadurch besonders 
einfach, da die GréBen 0 und dO@/d(x/l) bereits aus der Berechnung der Impulsverlustdicke (4) 
zahlenmaBig bekannt sind. Dariiber hinaus kann im Rahmen der Naherung der zweite Summand 
im Integranden der Gleichung (23) im ebenen Fall (R/l = 1) gegen den ersten vernachlassigt werden, 
da das Argument des In nahezu eins und A sehr klein gegen eins ist?. 

Als Konstanten zur Bestimmung von L treten in (21) A, C und statt B der Formparameter L,, 
auf. Nach Truckenbrodt ergeben sich A, B und C aus der Linearisierung einer Funktion K (Res; H), 
durch die eine geschlossene Integration der Differentialgleichung fiir L erst erméglicht wird’. Eine 
gewisse Unsicherheit in der Bestimmung dieser Konstanten geht jedoch, wie man aus dem Aufbau 
der Gleichungen (21) bis (23) und durch Beispielrechnungen feststellen kann, bei den Konstanten A 
und C praktisch in das Ergebnis nicht ein. Die Konstante B, die in die GréSenordnung von L fallt, 
und damit mafgeblich eingeht, kann durch die Einfiihrung des Formparameters L,, der Gleichdruck- 
grenzschicht unbestimmt bleiben. Wegen des eindeutigen Zusammenhanges L(H) kann L,, aus dem 
Formparameter H,, der Gleichdruckgrenzschicht bei glatter und rauher Wand ermittelt werden 
(Abb. 3). 

Gleichung (21) fiir den Formparameter L gilt gleichermaBen sowohl fiir laminare wie fur turbu- 
lente Grenzschicht, wobei lediglich die Konstanten verschieden sind. Die Werte dieser Konstanten 
sind in Tabelle 3 zusammengestellt. 


1 Dieser Schritt wurde von L. Speidel in einem unveroffentlichten Bericht des Institutes fiir Stromungs- 
mechanik der T. H. Braunschweig bereits durchgefihrt. 
* Auch im rotationssymmetrischen Fall kann dieses Glied dadurch vernachlassigbar klein gemacht werden, 


da man die Bezugs-Reynoldszahl anstatt mit | mit einem mittleren Radius der Korperkontur bildet, wie man 
anhand von Gleichung (20) erkennt. 


3 Siehe FufSnote 1 von Seite 82, dort Abb. 10 und 11. 
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Tabelle 3. Werte der Konstanten in Gleichung (23). 


Reibungsschicht 


laminar 


| glatt | 0,0685 Re, — 0,486 | 
turbulent = | rauh | 0,0665 (I/k,) — 0,340 0,03 | 


b) EinfluB der veranderlichen Stoffeigenschaften. Wird das vorstehende Gleichungs- 
system formal auf die Grenzschicht mit veranderlichen Stoffeigenschaften iibertragen, indem fiir die 
Impulsverlustdicke nunmehr @ nach (7) eingefiihrt wird, so bleiben die Gleichungen (21) bis (23) 
unverandert erhalten. Dies bedeutet wiederum, daB der Einflu® der Netmidellichen Stoffeigen- 
schaften der Grenzschicht im wesentlichen in einer Anderung des Formparameters L,, der Cleich- 
druckgrenzschicht zum Ausdruck kommt. An Hand von (22) ergabe sich fiir diese Anderung 


c fo.) 
L.. —(Lx)r=T. = Aln ane a (24) 
Im Jaminaren Fall (A = 0) ware danach L,, unabhangig von einem TemperatureinfluB auf die Grenz- 
schicht, wahrend im turbulenten Fall (A = 0,03) eine geringe Erhéhung der Abliseempfindlichkeit 
bei Erhéhung der Temperatur innerhalb der Grenzschicht vorhanden ware. Dieses Ergebnis steht 
zumindestens fiir warmeundurchlassige Wand in qualitativer Ubereinstimmung mit bisherigen Er- 
kenntnissen. Im Rahmen des Naherungscharakters der gesamten Rechnung kann dieser geringe 
Hinflu8 jedoch wohl vernachlassigt werden. Ubrig bleibt demnach lediglich der Temperatureinflu8 
auf die Impulsverlustdicke, die wiederum in die Berechnung von L eingeht. : 

Der Zusammenhang LH ; H), wie er durch (18) gegeben ist, wurde mit Hilfe der Definitions- 
gleichungen fiir 6*, # und 6 gewonnen, die fiir konstante Stoffeigenschaften gelten. Der Zusammen- 
hang mit den entsprechenden Gréfen bei veranderlichen Stoffeigenschaften laBt sich naherungsweise 
mit Hilfe der Bezugstemperatur Tp angeben. Fiir die Verdrangungsdicke gilt zunachst streng 

5 3 5 


a POU Fe. __ u(y) _ fey) uy) /e) * 5 
=f Me) = [ (“re + [ Seve (ayy 1) = rare + 3r- 05) 
0 0 0 


Ferner gilt naherungsweise fiir die _Impulsverlustdicke 


é é 
_ fey) uy) u(y) a rTP uy) u(y)’ T 
o=| 0(6) U ee oe mr | “a — “7 )y = 7, Or-t= ; (26) 
6 6 
und fiir die Energieverlustdicke 
é é 
= u? 2 oe u? LB ie ; 
5 fre harm E [Ba Eder 
Daraus folgt ; ; 
H & (A)r-1e (28) 
und mit (25) : 
Hw 2 (H)r-1. + Hy Hl. : (29) 


Mit zunehmender Temperatur in der Grenzschicht wachst also die Verdrangungsdicke der Grenz- 
schicht bei gleicher Impulsverlustdicke stark an, wahrend die Energieverlustdicke unverandert bleibt. 
Um aus dem Formparameter L bei einer Grenzschicht mit veranderlichen Stoffeigenschaften das 
Grenzschichtdickenverhaltnis H zu gewinnen, hat man somit zunachst tiber den universellen Zu- 
sammenhang L(H) ein aquivalentes (H)7_7,, zu ermitteln und damit nach (29) H zu berechnen. 


c) Bestimmung der Funktion H(L). In die Beziehung H(L), die durch (18) gegeben ist, 
geht die Abhangigkeit H(H) ein, fiir die im turbulenten Fall einige Unsicherheit vorliegt?. Die von 
J. Rotta angegebenen zweiparametrigen Beziehungen A(H; cy) und c(H; Res) baw. c;(H; 3/ks) lassen 
sich auf Funktionen H(H) mit dem Parameter Rey baw. /kzg fiir hydraulisch rauhe Wand umrechnen. 


1 Zur Vereinheitlichung des Rechenganges kann K(L) im laminaren Fall abweichend von Truckenbrodt durch 
eine einzige Gerade ersetzt werden, so daf sich nur ein Wert von C, dafir aber eine Nullpunktsverschiebung fiir 
L ergibt, die bei der Beziehung H(L) zu beachten ist. 

2 Siehe FuBnote 1 von Seite 82, dort Abb. 5. 
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Dabei zeigt sich, daB diese Funktionen im Bereich der praktisch interessierenden Parameterwerte 
nur so wenig voneinander abweichen, da sie durch eine mittlere Kurve ersetzbar sind, die in Abb. 5 
enthalten ist. Dieser Zusammenhang, der nunmehr sowohl fiir hydraulisch glatte als auch fir rauhe 
Wand giiltig ist, wurde fiir die 
Funktion L(H) zugrundegelegt. 
Die Integration von (18) ergibt 
dann die ebenfalls in Abb. 5 
eingetragene Beziehung L(H), 
bei der abweichend von Trucken- 
brodt die Integrationskonstante 
so bestimmt wurde, daB L = 0 
bei H = 1,3 zu liegen kommtt. 
Gegentiber Truckenbrodt ergibt 
sich in der Nahe der turbulenten 
Ablésung eine merkliche Ab- 
weichung. Das maximal erreich- 
bare L des Ablésungspunktes 
liegt entsprechend Hy, = 2,4 bei 
L,—=— 0,264. In Tabelle 4 sind 
die Werte der Funktionen H(H) 
und L(H) fir laminare und 
turbulente Grenzschicht zusam- 
mengestellt. 


Damit kann die GroBe L, fur 
die turbulente Grenzschicht bei 
Gleichdruck aus dem in Abb. 3 
enthaltenen Grenzschicht- 
dickenverhaltnis H,, ermittelt 
werden. Ein Vergleich dieser 
Werte mit den sich aus (22) 

“4 ergebenden Werten unter Ver- 


Ah 18 


O12\- 46 


008- 745 


~g0ul- 42 


-G08\- 47 


-072— 40 
40 42 14 1 j 
. ie : Bros 108 eal a ah 56 wendung der von Truckenbrodt 
bb. 5. Die zugrundegelegte Funktion H(H) und die daraus berechnete Funktion L(H). + 43 
(— — — L(H) nach Truckenbrodt). angegebenen Konstanten zeigt®, 
; a da lediglich beihohen Reynolds- 
Tabelle 4. Die Funktionen H(H) und L(A) fiir laminare vaion oe Cleidhcetnk one 


und turbulente Grenzschicht ‘ 
schicht nach Truckenbrodt etwas 


laminar? | Sepeeet kleinere L-Werte liefert. Da der 
L Verlauf von L,, tiber Re, un- 
mittelbar den EinfluB der Rey- 
ED) 1,59 0,034.7 1,20 1,855 0,127 noldszahl auf die Lage des Ab- 
ee Py rae nee oe ie lésungspunktes wiedergibt, wird 
25 157 00042 135 774 | 0,042 auch der turbulente Ablésungs- 
2,6 1,56 —0,0008 1,40 1,751 —0.079 punkt nach Truckenbrodt bei 
8 Ses —0,0083 1,45 | UPA N= DG hohen Reynoldszahlen etwas 
, nee Seyi ee . ea | Ba oe erreicht als bei Zugrunde- 
344 1,52 | —0,0173 1,70 1,045 | - —0,197 ogung Ser Lie" neve Oe 
3,6 1,52 —0,0182 1,80 1,622 | —0,217 fassers. Weiterhin erkennt man, 
4,038 | 1,52 | —0,0190 1,90 1,602 | —0,231 da bei sanftem Druckanstieg 
(Ablésung) ae es —0,242 die Reynoldszahl bzw. Wand- 
220 1558 ee rauhigkeit einen erheblichen 
2,30 1,548 0,260 EinfluB auf die Lage des tur- 
2,40 1,540 —0,264 bulenten Ablésungspunktes ha- 

(Ablésung) ben kann. 


+ Damit entspricht L = 0 etwa bei Re = 10? dem Gleichdruckfall (Abb. 4), wahrend nach Truckenbrod 
L=0 bei H=1,4 lag. Diese Nullpunktsverschiebung ist bei dem Vergleich i “Abb 5 beriicksi pig wenden elk 
; ee ies unter Verwendung der Hartree-Profile. : . ee ee 

ach (22) ist hierbei zu setzen: A = 0,07 log e, B = 0,23, C = 4, Cfo nach (14), fe ist di 
des Nullpunktes fiir L gegeniiber Truckenbrodt zu beriicksichtigen. pe PSS He: 


XXIX. Band 1960 


N. Scholz: Erganzungen zum Grenzschichtquadraturverfahren 91 


Tabelle 5. Rechenschema zur Berechnung laminarer und turbulenter Grenzschichten bei hydraulisch glatter und rauher 
Wand fiir ebene und rotationssymmetrische Begrenzungen. 


ee 


| 


1 Ifd. Nr. des Rechenpunktes ty 23 

2 relative Lauflange | x/l 24 

3 relativer Wandradius Ril ton 

4 relative AuBengeschwindigkeit WiiGics 

5 relative Wandtemperatur Ty/T x ae 

x—l] Daa 
- 6 Cire M2, - (4)? 
x—1 / 

i et 3 MS / (1+ (6)) TEs 28 

8 (7)U@—1) 0/0 x 29 

9 r (6) 
10 | (7) (9) +0) T,/Ta | 30 
11 (5)—(10) 31 
12 (11)/(7) 32 
13 (6) + (12) Hr 33 
14 3 + 2-(13) c 34 

n 
15 255 (1 + (13)) 35 
16 (14) + (15) q 36 
17 a: (7): (9) 37 
18 b- (12) 38 
19 (7) + (17) + (18) Tp/T x 39 
20 (18)* 40 
21 (3)? 41 
(4)4 42 


(7) - (8) - (20) - (21) - (22) dO/d = 
Oy + J (23) d(x/l) | ) 
(24)1/P | 

(aye | 

(25)/(3) + (8) (26) ©) | Fok et|2) 
(24) | 

d- (23) - (28) + (33) dA/d = 
Ay + f (29) d(«/l) beau 
(30)/(24) - (28) 

e+ In (24) 

In (4) 

(31) + (32) + (33) Sis 
Git. | ‘3 
aus H(L) (Tabelle 4) | H 
Hool(cfo)r= 70/21 27/36) Ol 
aus H(L) (Tabelle 4) | (D)r=T. 
(19) - (38)/(7) 

(13) + (36) 

(39) + (40) test 
(37) - (41) aap sO" /E 


Tabelle 6. Spezielle Werte im Rechenschema fiir laminar und turbulente Grenzschichten 


laminar 


0,9 
k wo—l 
p 2 
q 5+4 Hr 
m 6 
d 4 
e 0 
1 Dee 0 
Hoo 1,56 


(cfoo)T=Too | 1,328 Rep 


turbulent 


hydraulisch glatt 


1,12/log Rei 


0,0685 Re; — 0,486 


hydraulisch rauh 


0,0025 
| 0,0665 (I/ks) — 0,340 


aus H(L) fiir L = Le (Tabelle 4) 


0,455 (log Rez)—?*8 


| 0,418 [2 + log (I/ke)}- 258 
aq 
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Fiir L..(Re)) baw. L,o(I/ks) lassen sich in Anbetracht der sonstigen Vereinfachungen fiir praktische 
Rechnungen ausreichende Interpolationsformeln angeben (vgl. Tabelle 3 und Abb. 4). 


5. Rechenschema. In Tabelle 5 wurde auf Grund der hier abgeleiteten Gleichungen ein Rechen- 
schema zusammengestellt, das universell fiir laminare und turbulente Grenzschichten im hydraulisch 
glatten und rauhen Bereich bei ebenen oder rotationssymmetrischen Begrenzungen anwendbar ist. 
Bei ebenen Begrenzungen ist hierin iiberall R/J =1 zu setzen. Strémungsform und Wandbeschaffen- 
heit gehen durch zahlenmafig verschiedene Konstanten ein, deren Werte aus Tabelle 6 zu entnehmen 
sind. 

Die Zeilen 1 bis 5 enthalten die vorgegebenen Werte der Geometrie der Wandbegrenzung, der 
AuBengeschwindigkeitsverteilung und der Wandtemperaturverteilung. Wegen der erforderlichen 
Integrationen in Zeile 23 und 30 sind diese Werte zweckmafig fiir eine ungerade Anzahl aquidistanter 
Rechenpunkte x,/1 vorzugeben. . 

Fiir eine Grenzschicht mit konstanten Stoffeigenschaften (Ty = T,=T) kommen die Zeilen 6 
bis 19 sowie 25 und 39 bis 41 in Fortfall, wodurch sich das Rechenschema wesentlich verkirzt. Bis 
zur Zeile 34 ist die Rechnung unabhangig von der speziellen Reynoldszahl oder Wandbeschaffenheit, 
wenn man von dem sehr geringen Einflu8 auf die Exponenten p und q absieht. Im turbulenten Fall 
braucht also die Rechnung bei einer Anderung der Reynoldszahl oder der Wandbeschaffenheit nur 
ab Zeile 35 wiederholt zu werden. Die in Zeile 23 bzw. 30 auftretenden Integrationskonstanten 0, 
bzw. A, sind mit den in (8) bzw. (23) vor dem Integralglied stehenden Integrationskonstanten identisch. 
Man erhalt sie im Rechenschema durch Vorgabe der Anfangswerte der Grenzschicht 7) und Ly, und 
sinngemaBer Riickrechnung der Zeilen von unten her. Die Integrationen in Zeile 23 bzw. 30 sind bei 
Verwendung einer Tischrechenmaschine zweckmaBig mit Hilfe der Simpsonregel in Verbindung mit 
der 3/8-Regel leicht zu_gewinnen. 


6. Zusammenfassung. Aufbauend auf dem von E. Truckenbrodt angegebenen Quadraturverfahren 
fiir laminare und turbulente Grenzschichten werden folgende Erganzungen mitgeteilt : 

1. Durch eine Beriicksichtigung des értlichen Formparameters bei der Quadratur der Impuls- 
verlustdicke wird eine bessere Ubereinstimmung mit dem Experiment erzielt. 

2. Durch Einfiihren eines geeigneten Mittelwertes der Temperatur innerhalb der Grenzschicht 
(Bezugstemperatur) gelingt es, auch die Grenzschicht mit veranderlichen Stoffeigenschaften (Warme- 
tbergang, Reibungswarme) durch ein einfaches Zusatzglied bei der Berechnung der Impulsverlust- 
dicke zu beriicksichtigen. Fiir die Verdrangungsdicke wird eine Abschatzung angegeben. 

3. Der EinfluB der Reynoldszahl 1aBt sich soweit eliminieren, daB er in die erforderlichen Quadra- 
turen nicht eingeht. Durch Einfiihren einer der Wandrauhigkeit aquivalenten Reynoldszahl wird 
auch der, Fall der hydraulisch rauhen Wand eingeschlossen. Fiir einen vorgegebenen Druckverlauf 
langs einer Grenzschicht kann der gré8te Teil der Berechnung der Impulsverlustdicke und des Form- 
parameters unabhangig von der Reynoldszahl und der Wandrauhigkeit durchgefiihrt werden. 

4, Durch Elimination des Anteiles der Gleichdruckgrenzschicht im Truckenbrodtschen Grenzschicht- 
formparameter kann eine gewisse Unsicherheit in der Berechnung des Formparameters beseitigt und 
das Ergebnis fir die Gleichdruckgrenzschicht im Falle hydraulisch glatter und rauher Wand dem 
experimentellen Befund besser angepaBt werden. Hierdurch darf auch bei veranderlichem Druck- 
verlauf eine gréRere Treffsicherheit der Rechenergebnisse erwartet werden. 

5. Der Zusammenhang zwischen dem Truckenbrodtschen Formparameter L und den Grenzschicht- 
dickenverhaltnissen H und H wird fiir die turbulente Grenzschicht unter Verwendung der Rottaschen 
Ergebnisse neu berechnet. 

6. Fiir die numerische Durchfiihrung des Verfahrens wird ein einheitliches Rechenschema an- 
gegeben, das folgende Falle umfaBt: 

Grenzschicht laminar oder turbulent, 

Wandbegrenzung eben oder rotationssymmetrisch, 

Wandbeschaffenheit hydraulisch glatt oder rauh, 

Stoffeigenschaften konstant oder temperaturabhangig. 


Das rein digitale Berechnungsverfahren kann leicht tabellarisch von Rechenhilfskraften durchgefiihrt 
werden und gestattet ebenso wegen seines einheitlichen Aufbaues eine universelle Programmierung 
fiir elektronische Rechenanlagen. 


(Eingegangen am 16. Mai 1959.) 
Anschrift des Verfassers: Privatdozent Dr.-Ing. N. Scholz, Miinchen-Pasing, Backerstr. 57. 
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On the Transverse Flexure of a Semi-Infinite Plate with an Elliptic Notch 
By S. Shioya 


1. Introduction. The problems of the stress distribution in a thin plate having notches and sub- 
jected to plain bending were first studied by Lee! who treated the case of two symmetrically disposed 
hyperbolic notches. Neuber? also discussed similar problems. 

Recently, the author? gave the theoretical solution for an infinite strip containing two symmetri- 
cally placed semicircular notches under plain bending. 

The present paper gives a theoretical solution for a semi-infinite plate notched by an elliptic notch 
and transmitting a constant bending moment. The analysis is developed on the basis of the Poisson- 
Kirchhoff theory* of thin plates and by means of Muskhelishvilli’s complex variable method’, The 
method of perturbation is adopted for the determination of parametric coefficients involved in the 
solution. 

Numerical examples are also worked out in some detail, in order to clarify the effects of the notch. 


2. Preliminary results. When the middle plane of an isotropic homogeneous thin plate free from 
lateral load lies originally in the complex plane of z = x +7 y, its deflection w normal to the z-plane 
is governed by the biharmonic equation, provided the plate is assumed subject to the limitation of 
the Poisson-Kirchhoff theory of thin plates. 

As is well known, the solution of the biharmonic equation can be expressed in terms of two analytic 
functions @,(z) and w,(z) of the complex variable z: 

w= RefEqye) + 2] with 2) =f Ge) de. (a) 
Overbars are used to indicate conjugate complex quantities. : 
The bending moments M,, M,, twisting moment M,, and shearing forces Q,, Q,, all per unit of 
are length, are expressible by the formulas 
M, + My=—2D(1 +») [9,@) + 9,@)], 
M,— M,+2iM,, =—2 D(1—») [z¢,@) + ,@)]. (2) 
Q,—iQ, =— 4 D¢,(z) 
in which D = E h3/12 (1 —v?) is the flexural rigidity of the plate; E,» and h being, respectively, 
Young’s modulus, Poisson’s ratio and thickness of the plate. 
By means of the mapping function 
c= o(C eG 0 e° (3) 
curvilinear coordinates (g, ) are introduced into the z-plane. Substituting equation (3) into equation 
(1), we now have the deflection of the plate given by 


w = Re[a(6) o(0) + x0) (4) 
PS) = Pilo(2)] = vil2) » ll) = f pC) w'(6) al 
w(C) = vilo(d)] = vil) - 


‘The bending moments M, and Mp, twisting moment M,o, and shearing forces Q, and Qo, all 
per unit of arc length with reference to the new coordinates (9, @), are expressible by the formulas 


M, +My = —4D(1+>») Re[D(0)], 
» &(C) OO) + WE) w’(0) 
0? w’(¢) f ©) 


where 


and 


Ms, — M,—2%i Myo = 2D(1—v)o 


®'(C) 
bre = Ego @) 


G. H. Lee, J. Appl. Mech. 7 (1940) p. A-53. 
H. Neuber, Z. angew. Math. Mech. 20 (1940) S. 199. 
O. Tamate and S. Shioya, Trans. Japan Soc. Mech. Engrs. 24 (1958) p. 809 (in Japanese with English summary). 
S. Timoshenko, Theory of Plates and Shells, New York 1940. j ; 
E. g. G. N. Sawin, Spannungserhohung am Rande von Léchern, Ubersetzung aus dem Russischen ins Deut- 
sche von H. Neuber und O. Everling, Berlin 1956. 
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in which 
6 9O)2 U 5 NE 
2) = ©. oe) 
For a free boundary of a plate, the boundary equation is written as 
OMe 
M,=0, Pa eae hala 0 (6) 


where M,,, M,,,, and Q,, are, respectively, the bending moment, twisting moment, and shearing force, 
per unit of arc length of the plate boundary. Here n denotes the outward normal to the boundary 
and s the arc length. The boundary condition of equation (6) is transformed into? 

— 190+ ZOFO+HO=1G004G 
where 
a Way 
a 
and C, is a real constant and C,a complex one. When the 
region occupied by the middle plane of the plate is simply 
conuected, these constants may be taken to be zero. 

With these formulas, we shall seek the solution for the 
case where the semi-infinite plate, having an elliptic notch, 
is subjected to the plain bending of intensity M, = M, 
at infinity, about x-axis. In this instance, we assume that 
the straight boundary of the plate and the rim of the notch 
are both stress free. 

3. Method of solution. Let the middle plane of the plate 
occupy the right half of the z-plane, with the exception of 
the interior of an elliptic notch whose center is at origin 
(Fig. 1). 

With the aid of the mapping function 


o)=RE+F) O<|mi<i  @ 


and by varying the values of the real constants R and m; the straight boundary of the plate and 
various shapes of an elliptic notch in the z-plane may be mapped onto the line 6 = + 7/2 andp = 1 
in the ¢-plane, respectively. Denoting the depth and the half width of the notch by x, and y, re- 
spectively, the ratio %)/y) for several values of m are given in Table 1. 


Table 1. Ratio xo/yo for the Values of m. 


é 


Fig. 1. Shape of elliptic notch. 


l Foes 
ie. oO. B he ol eg ee | Bot el Scot | 0.3. \29048 


| 
Xol¥o | 0.333 | 0.538 | 0.818 10 | 1.222 1.857 | 2.636 


Then, the problem to be solved is reduced to the construction of two function y(C) and y(C) which 
satisfy the following boundary conditions: 


@) 9=+m2, — 29+ 22 FO +O =0, (0) 
@e=1,  —F 9+ S296) +90 =0. (10) 


Equations (9) and (10) correspond to the conditions of the straight boundary of the plate and the rim | 
of the notch, respectively, which are both stress free. While, the stress resultants at the infinity | 
part of the plate are M, = M, = constant, with M, = M,,y = Q, = Q, = 0, which may be easily 
transformed into 
(ili) @—>0O, 
M, +- Me = M,, (11) 
Mp — M, —21 Myo = M,e?'° . 


(5) = *p(C) + Mol) | : 
p(C) = “p(C) + Mw?) 


Consider the functions 


(12) 
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where 1y(C) and 'y(¢) correspond to the deflection of the semi-infinite plate, without notch, subjected 
to plain bending and can be written as 


I x M, R om 
| 
M, R m 

7 es 
Thus, it is evident that equations (13) satisfy the conditions (iii) at infinity. No(f) and "y(C) are 
the respective parts of the functions which are analytic in the entire region of the semi-infinite plate 
surrounding the notch and therefore include only terms of negative powers of ¢. Taking into account 
the requirement of the symmetry of the problem, we may write %(C) and “y(¢) in the form of series: 


cdg Re NT case bon : 
HAO) = — pay 2, sat t+ PEs Ea +»), 


(13) 


M,R oj 
"y(0) =— pas BE) Oy (tant OME bay OM tA) (14) 


, Can Scie : doy —(2n+1 
ON oe Ee =] 


where a2, b2,, C2, and dz, are real constants. It is easily shown that these functions leave no stress 
resultants at infinity. Accordingly, the problem reduces to determine the coefficients a2,, ben, Can 
and d;,, so that the straight edge of the plate and the rim of the notch are stress-free, as given by 
equation (9) and (10) respectively. 
Substituting equations (13) and (14) into equation (9) and equating to zero the coefficients of 
various powers of 9, we obtain 
eee > da, =—— bin. (15) 


In a similar way, equation (10) gives 
) es ° 5 x 43 ° 5 
5: (i +m) (er? Be e— #6) __ a @n+de (kn +2)0__ p—i(2n + 2) ay fe 


~ ben : ; s 
se a (Gh ciel) 0) es i(2 mn + 1)9) ae 
2 (2n + ye ( 


EMER Sg, care thee 0 As 
is n=0 


from which it follows that 


+ (1—m) +> —m!) cos26+4+ >) (1 +m)(1— gta) anos Qn +16 + 


n=0 


. . 1 
Fa +m)(1 + aeeye) amen @nt+3)6+ Dil mn) Sd x 


x ben cos 206 + dle +m) +m) a 7y_| bon cos 2m +2)6=0, i 
= (1 +'m)? sin 2 0— >a +m—C—mg ty, |ensin@n +06 i 

= >a $m) —0— mg 5 Jazysin@n +3)0— DU +m)(1 tare) x 

X be,sin2n Raper +-m)(1 gay; bensin @n +2)0=0. 


Moreover, sines and cosines of odd multiples of 0 can be expressed in Fourier series of sines or cosines, 
respectively, of even multiples of valid between 0 = —x/2 and 0 = 2/2 


4 —1y . (—Its(2n4+)) 
cos (2 n $004 | ett Dees en eed cos 2s 6], 


(18) 
. (—1)»+s2s 
sin 2n+ 1)O0=7 2 G@n—2stl@nt+is¥)) 


cout 


sin2s 9. 


| 
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Substituting equations (18) into equations (17) and then equating to zero the coefficients of | 3 
cos 2s 0 and sin2 s 6, we get i 
= (1 == 


1 4 
Fm) +4+m) +m); th+4 datm oat perEy t= 
1 \| 
aur m®) O10 + [(L +m) +— Rpereh de m) -+(1—m) @s—ayef ee? + (19) |} 
je apt fentnen+a(i—2)+ae(r+ 4} 
+ 2 + mapas pGaats FHOn Fasten tas pyr" 
si (s=)), 
1 1 
eek Tea (1 +m)(1 a renrice. a 
8 oo 
—$ 2 |) 0 enea x 
(—1)"»+s4s(2n-+ 2) 2% a) (20) | 


x @n—2s+DCn—2s+3)Qn-is+ 1)(2n + ee 
where 6;,, is Kronecker’s delta. The set of equations (19) and (20) solves completely the problem of 
determining the parametric coefficients included in the functions "g(¢) and "y(¢). Then, if we give 
the values for m and ¢ in equations (19) and.(20), a formal solution of the set of equations can, for 
instance, be obtained by means of the method of successive approximation. But, it seems rather 
convenient to solve these equations by the method of perturbation, in which m and é are the per- 
turbation parameter. We shall write a formal solution of equations (19) and (20) in the form 


ioe) lo @) co co 
== Dy 2 2a ant ePs bon =. & OP) m" eP (21) 


r—0 spi 2 r—0 op —- 


where (a(?) and (b(?) are independent of m and ¢. Substituting these equations into equations (19) 
and (20), and then equating the coefficients of the same powers of m and « on both sides, we arrive 
at the following results: 


Mf . = 
>) Do MalP) = ip oa 2h (= 5, .) pe. Lh A Malp—1) (1 65== 63,») 2 Bee Maip— 2) 


co 


U 
2» 3, Map) = 751,22, — (1 Op) Dy 292 Male + (1 —dn,p — 83,9) X (22) 
x > (2k—2s +3)"P2Mae—9 (5 >1), 
oe Oa) = 216 nea, eee ee (1—6 
— 2n Dp een sy l,r 9 2,7} ¥2,p 3 2, p) x 


x nop DP) (DBP) 4 DEP) 


= ql 02>) x 


oo 


x Pea Male—) =D mo (—Yalp) TIO ag] 4 (1—6e,» — 63,5) X 


«| Seog a + Spy.e— rage]. 


ce 


1 3 1 al 
2» ® @2! Mah) = ren Cor + 6;,s) + <i (Or-—2),2 + O04 ™) Ong + (L0u5,)- 5 (23) i} 
PRET ae ao 2s+1 
es lee eae te 
1 ; 
2 (7; SRA NT gee ye eet = 04,5) xX (r= 1) 


| 52 apes e Malp—1) + = ee (Na?) aes 2) ie ae peg ae (1 ons 62 Fy ee 2 x 


x[ Se k— 2s + 3) "B22 ae—2) + 2 (2k—2s +3) 3p {vat (s>1). 
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yi ee Syl 
( b(P) ari — 5 O08 O1,p =m T(z. a 51.1) 0ap (1 Ory 62, p) (2 s 1) 
x » D3?) af?) +. .S (2 des ait) 1Y2s + 2 og] (24) 
n= n=0 Ks zB 


1 
()b(P) = — = 90,591,701, an | 60.b1 = ae s +1) 


1 
x B (0 —2), 7 Or+ 1), *) a= Ops 1), s (i = 60, s) -b brah 02,p 


ae 4 “ 
seth 01,7, 02,5) pot 3 ve Ls oor ples +) | D3? Mal) 
n= 0 


=e oa (ae oer 2, (2m +25 + 7) p2++2 Malp—}) 


n=0 
a a, (2n+2s +7) are rag (r => 1) (25) 

where 

ps — (—1)n+2+1 ’ 

an (Qn—2s+1)\Qnp2s+3)’ 
ryyo Cen ae 

an 2(2n+3)Qn+5)’ 
Iys2s __ Caijerett 

2n  (Qn+2s+1)@n+2s+3)Qn+2s5+5) Ge), 
nysas (—1)p>tst1 [4n? + 8n—(2s + 1) pe) 

an (2n—2s—I1)2n—2s+3)\(Qn+2s+1)@n+2s+5)’ 
nyse (—1+1(2n—3) | 

an 2Qn—lQ@n+lent+s)’ ® 
myo (—1)pts412[4n? + 8n4+3(25+1)] 

an (2n—2s—1)Q@n—2s+ 1) Qn—2s+3)Qn+2s+3)Qn+es}5) S21.) 


With theserelations, the coeffici- 
ents a?) and ()b(?) in the Series (21) 
can be determined successively for 
PO bee anda = 1,253, ..2, 
respectively. In particular, the co- 
efficients a?) and (6), corres- 
ponding to the solution for the semi- 
tnfinite plate! notched by a single 
semicircular notch, were already 
obtained by the author. It may be 
noted that in the infinite sets of 
equations (22) and (23) the coeffi- 
cients of “a{?) are invariant for all 
values of r and p. Hence, if we solve 
these equations once, the solutions 
in successive cycles are obtained 
only by the calculations about the 
constant terms. 

_ For brevity, no proof of conver- 
gence is pursued here, but it will be 
illustrated by the numerical ex- 
amples which follow. 05 


0 40 40 30 40 50 6,0 40 80 
4, Numerical Examples. The 
foregoing solution is worked out in 
some detail. The calculations of 
a?) and “b(?) are carried out until r = 7 and p = 5, but their values are not shown here for brevity. 
By using these values, it is now rather straightforward to calculate the stress resultants in the neigh- 


Fig. 2. Stress resultant M,, on the minimum section (» = 0.3). 


oy 


1 Q. Tamate and S. Shioya, |. c. 
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bourhood of the notch. On working out the values of M, and (Q, — 0M,,/@s) ee es te ein 
on the rim of the notch, where s denotes the arc length of the notch, the greatest value ound is Pe 
than 0.04 M, for m = 0.45 and vy = 0.3 (M, for 0 = 0). Therefore, the values of parametric coeffi- 


0 40 40 30 40 50 60 40 80 


Fig. 3. Stress resultant M, on the minimum section (v = 0.3). - 
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Fig. 5. Stress resultant My on the straight edge (» = 0.3). 


cients obtained are thought to be sufficient to get the stress resultants with satisfactory approximation. 
When m is smaller, the approximation is more satisfactory. 


The values of the stress resultants M. and M,, along the axis of symmetry, 0 x, are shown for 
some values of m in Figs. 2 and 3. It is seen that the maximum values of M y/M, occur at the apex 
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of the notch and they vary almost linearly with the ratio %9/Yo as shown by dotted line in Fig. 2. 
Fig. 4 shows the values of the stress resultants M,/M, round the notch for several values of m. The 


stress resultants at the straight edge are shown in Fig. 5. 


In order to study the effects of the notch, we shall consider the stress concentration factor hee 


i. e. the ratio of the maximum bending moment at the apex of the notch to 
the rim of the notch the maximum bending moment is given by 


M,. Since M, = 0 on 


Ee ee ae NN) (0g? —l i 
ae Mr ee 
where 
")M() — Mal) + MBP) . (28) 
d( P) + lp?) 


With this equation, the stress concentration factor k, are calculated and plotted versus ¢ in Fig. 6. 
For the limiting case of m = 0, the results agree with those for a semicircular notch. 


o4 a3 Q2 0 ky 
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Fig. 6. Stress concentration factor kp versus e. 
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Fig. 7. Stress concentration factor kp versus 0’/x). 


The stress concentration factor are also plotted against 9’/x) in Fig. 7, where 0’ is the radius of 


curvature of the notch at its apex. 


In conclusion, the author wishes to express his hearty thanks to Professor Emeritus Seiichi Higuchi 
and Professor Osamu Tamate of Tohoku University for their kind advices and interests throughout 


the progress of the present investigation. 


(Eingegangen am 25. Mai 1959.) 


Anschrift des Verfassers: Shunsuke Shioya, Tohoku University, Department of Mechanical Engineering, 


Sendai (Japan). 
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Der Einflu8 der Kraftwagenfederung auf die Lenkstabilitat 
; Von G. Mitterlehner 


1. Einleitung. Diese Untersuchung ist die Fortsetzung meiner fritheren Arbeit ,,Die Lenkstabilitat 
des luftbereiften Kraftwagens gegeniiber kleinen Stérungen“.t Im folgenden wird vorausgesetzt, 
daB der Leser mit dem Inhalt der erwaihnten Arbeit und den dort verwendeten Bezeichnungen ver- 
traut ist; doch soll das darin behandelte Problem nochmals kurz erlautert werden. 

Ein Vorderrad eines Kraftwagens iiberfahrt eine kleine Unebenheit der Fahrbahn, wodurch die 
beiden durch die Spurstange miteinander verbundenen Vorderrader etwas um ihre Lenkzapfen ver- 
schwenkt werden. Auf diese Weise wird ein lineares System von drei Freiheitsgraden — seitliches 
Schieben und Gieren um die Hochachse des ganzen Fahrzeugs sowie Lenkausschlag der Vorder- 
rader — gestért, das auf diese Stérung ahnlich wie ein Kraftwagen reagiert. Zunachst wurde die 
Stabilitat bei festgehaltener Lenkung und elastischem Lenkgestange untersucht. Dann wurde das- 
selbe fiir losgelassene Lenkung durchgefiihrt, die allein dadurch rickgestellt wird, daB die von der 
Fahrbahn auf die eingeschlagenen Vorderrader wirkenden Seitenkrafte infolge des Nachlaufes Mo- 
mente um die Lenkzapfen erhalten.? In beiden Fallen ist zum Einhalten einer geraden Bahn eine 
kleinere Korrektur durch den Lenker unvermeidlich, was bei der Definition der Lenkstabilitat? 
beriicksichtigt wurde. Die Rechnung zeigt*, daB ein solches System vor allem infolge der Entlastung 
der Vorderrader durch das Moment des Luftwiderstandes um eine Querachse in der Fahrbahnebene 
mit zunehmender Fahrgeschwindigkeit unstabil wird. Die Stabilitatsgrenze liegt jedoch bei einer 
zu hohen Geschwindigkeit und ein vergleichbarer Kraftwagen wird erfahrungsgemaB schon wesentlich 
frither unstabil, was offenbar auf den dabei vernachlassigten Einflu8 der Federung zuriickzufiihren ist. 

Um der Wirklichkeit naher zu kommen, wird in der vorliegenden Arbeit ein System von vier 
Freiheitsgraden verwendet, wobei zu den drei vorerwahnten das Wanken des Fahrzeugs um eine 
Achse in der Langsrichtung hinzutritt. Damit lassen sich auch die Kreiselmomente der Rader be- 
riicksichtigen, die die Stabilitat eines schnellen Kraftwagens stark beeinflussen. Von den ,,unge- 
federten“ Bauteilen der Rader und Achsen werden nur die Drehmassen um Achsschenkel und Lenk- 
zapfen, nicht aber ihre Tragheit bei den annahernd geradlinigen Federbewegungen in Rechnung ge- 
stellt. Obwohl die Translationsschwingungen der Rader die Bodenhaftung und damit auch die Sta- 
bilitat beeinflussen®, tritt dieser EinfluB besonders bei hoher Fahrgeschwindigkeit und guter Fahr- 
bahn gegeniiber dem der Kreiselmomente vor allem der Vorderrader zuriick. Die aus dem Verschwinden 
der Systemdeterminante folgende charakteristische Gleichung wird mit Hilfe eines neueren in der 
Regelungstechnik gebrauchlichen Stabilitatskriteriums ausgewertet, das in diesem Fall recht rasch 
zum Ziel fiihrt. Die Stabilitatsgrenze liegt bei diesem System in der Nahe der Héchstgeschwindigkeit 
des vergleichbaren Fahrzeugs. Der Zustand von drei verschiedenen Kraftwagen wird bei der jeweils 
angegebenen Héchstgeschwindigkeit untersucht. Dabei ergibt sich Stabilitat gegeniiber kleinen, von 
der Fahrbahn ausgehenden Stérungen nur als Folge sorgfaltiger Abstimmung des Lenksystems und 
der Federung. Bei ausgefiihrten Kraftwagen liegen die Verhaltnisse infolge der progressiven Nicht- 
linearitat vieler Federungen etwas giinstiger, da diese zur Stabilisierung eines im kleinen unstabilen 
Systems fiihren kann. Dem wirkt jedoch der Kinflu8 der in Wirklichkeit degressiv nichtlinearen 
Reifenkennlinie® entgegen. SchlieBlich wird noch der Einflu8 der sog. Vorspur behandelt. 


2. Die Lage der Momentanachse des Wankens. Die Federung beeinfluBt die Stabilitat hauptsachlich 
dadurch, daB sie ein Wanken des Aufbaus um eine Achse in der Langsrichtung erméglicht. Eine 
solche Bewegung wird haufig als Rollen des Aufbaus bezeichnet, was aber hier zu Verwechslungen 
mit dem Rollen der Rader fihren kénnte. 

Die Lage der Wankachse fiir verschiedene Radaufhangungen hat zuerst S. de Lavaud? ermittelt. 
Bei gleicher Lagerung der Vorder- und Hinterrader zeigt sich in der Vorderansicht der Spurpunkt 


1 G. Mitterlehner, Ing.-Arch. 27 (1959) S. 88. 

2 S. FuBnote 1, dort Abb. 1 und 6. 

3S. FuBnote 1, dort Abschnitt 8. 

4S. Fu8note 1, dort Abb. 8. 

5 H. Kurz, Automobiltechn. Z. 60 (1958) S. 124. 
8 S. FuBnote 1, dort Abb. 2. 


? Sensaud de Lavaud, Le probleme de l’independance des roues d’arri av it i i 
1929 (im Selbstverlag des Verfassers). : SRDS HAM a te 
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der Momentanachse des Wankens mit der Bildebene als Momentanpol M (Abb. 1). Fiir kleine Seiten- 
neigung des Aufbaus aus seiner statischen Gleichgewichtslage und bei Vernachlassigung der ver- 
schiedenen Eindriickung der Luftreifen findet man M als Schnittpunkt der i @itadingelinie des 
Stiitzpunktes 4 mit dem Momentanpol M eines Rades und der Mittelebene des Aufbaus. H. Maruhn} 
hat auf die gleiche Weise den Momentanpol fir zahlreiche Arten der Radaufhangung aufgesucht. 
Fiir kleine Seitenneigung des Aufbaus ist das Riickstellmoment der Federung dem Noieiigesiikel 
annahernd proportional. R. Eberan v. Eberhorst® gibt die zugehérigen Konstanten fiir verschiedene 
Radaufhangungen an. Mit der so berechneten Seitenneigung und den daraus entstenden Anderungen 
der Raddriicke und der Schraglaufwinkel der Rader hat derselbe Verfasser die Stabilitat des Kraft: 
wagens beim Durchfahren einer Kurve behandelt. Fiir gréBere Seitenneigung des Aufbaus kann der 
Momentanpol nicht mehr auf die oben beschriebene einfache Weise gefunden werden. Die dann unter 
Umstanden recht komplizierte Berechnung der Seitenneigung zu gegebener Fliehkraft und Harte der 
Federung wurde von E. Marquard® fiir verschiedene Radaufhangungen durchgefiihrt. 

Wie vorher fiir das Schieben und das Gieren des Fahrzeugs, sowie fiir die Lenkschwingung werden 
auch fiir die Seitenneigung des Aufbaus kleine Amplituden vorausgesetzt. Damit kann der Momentan- 
pol des Wankens wie in Abb. 1 bestimmt und das Riickstellmoment der Federung der Seitenneigung 
proportional gesetzt werden. Der bei Personenkraftwagen vorherrschende Fall verschiedener Rad- 
aufhangung fiir Vorder- und Hinterrdder fiihrt auf eine 
Wankachse, die nicht parallel zur Langsrichtung des 
Fahrzeugs ist, sondern von vorne nach hinten ansteigt. 
Er wird in einem Berechnungsbeispiel behandelt. 


verkiirzte 
Parallelfuihrung Pendelachse Starrachse 


Abb. 1. Momentanpol des Wankens bei verschiedenen Radaufhangungen. Abb. 2. Kriafte am Fahrzeug beim Durchfahren einer Kurve. 


3. Verwendete Bezeichnungen. An dieser Stelle werden nur jene Bezeichnungen erklart, die in de 
vorliegenden Arbeit erstmals vorkommen. Beziiglich der anderen sei auf Abschnitt 5 der ersten 
Arbeit verwiesen. Ein GroBteil der erstmals auftretenden Gréfen ist aus Abb. 2 zu ersehen. Es 
bedeuten: 

B den Riickstellbeiwert des Wankens, 
cp die Federkonstante der Hauptfeder fiir ein Rad, 


CR die Federkonstante eines Luftreifens, ; 
c, . die resultierende Federkonstante von Hauptfeder und Luftreifen, 


d den Dampfungsbeiwert des Wankens, 

af die Federweite, 

F die Federkraft der Hintereinanderschaltung von Hauptfeder und Luftreifen, 

h den Hebelarm der Fliehkraft um die Wankachse, 

ky den Dampfungsbeiwert eines StoSdampfers, 

I, die Drehmasse des Aufbaus um seine Langs-Schwerachse, 

I,, die Drehmasse eines Rades um seine Achse, 

y die Koordinate der Seitenschiebung des Momentanpols vom ortsfesten System, 
y* die Koordinate der Seitenschiebung des Schwerpunkts vom ortsfesten System, 
A den Winkel der Seitenneigung des Aufbaus, 

Ap den Kippwinkel eines Rades, 

Lb =p; = den Faktor der Massenverteilung um die Hochachse. 


Die Indices i, bzw. a bedeuten ,,auf der kurveninneren“, bzw. ,auf der kurvenduBeren“ Seite 
des Fahrzeugs. ; 
1 H. Maruhn, Die Grundlagen der Federung von Automobilen, Berlin 1932. 


2 R. Eberan v. Eberhorst, Automobiltechn. Z. 55 (1953) S. 246. 
8 J. Marquard, Automobiltechn. Z. 51 (1949) S. 31 u. 52 (1950) S. 34. 
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4, Das Wanken des Aufbaus. Da das Wanken im allgemeinen nicht um die Langsschwerachse des 
Aufbaus mit dem Spurpunkt © in der Bildebene von Abb. 2, sondern um die Momentanachse mit 
dem Spurpunkt M erfolgt, vergréBert es die Seitenschiebung des Schwerpunkts. Wahrend beim 
ungefederten Fahrzeug die Seitenschiebung fiir jeden Punkt der Hochachse des Aufbaus dieselbe 
war, ergibt sich jetzt fiir 1 <1 die Seitenschiebung des Schwerpunktes aus der des ,.Momentanpols“ 
und der Seitenneigung des Aufbaus zu 
y* =v +h. (1) 
Im bisher behandelten System! waren y* und y identisch. Es erweist sich als praktisch, die Koordinate 
y beizubehalten, die aber jetzt nicht mehr Schwerpunktskoordinate ist. Dann mu im Tragheitsglied 
fiir die Seitenschiebung y durch y* = y + hA ersetzt werden, wodurch bereits eine Kopplung zwi- 
schen Schieben und Wanken auftritt. 

Da der Momentanpol des Wankens M gegeniiber dem ortsfesten System die Beschleunigung y 
hat, setzt man den Drallsatz am besten um den Schwerpunkt © an. Wenn man voraussetzt, daB 
die Wankachse eine Tragheitshauptachse des Aufbaus ist, kommt fiir 1 < 1 


LA Sw 282) 25 eae ee (2)2 


Das erste Moment rechts in (2) stammt von der Summe der Federkrafte unter den statischen Rad- 
lasten, die negativ gleich dem Gewicht des Aufbaus ist und lotrecht durch M geht. Das zweite ist 
das Riickstellmoment der Federung und das dritte das Moment der Summe der Seitenkrafte der 
Luftreifen, die in M an den Aufbau iibertragen werden. Das Vorzeichen des letzteren in (2) gilt fir 
den Fall, daB y und ¥ den gleichen Richtungssinn haben, wahrend in Abb. 2 y und y einander ent- 
gegengerichtet angenommen wurden. 

Wenn man die Elastizitat der Luftreifen beriicksichtigt, wird die resultierende Federkonstante 
fiir die Hintereinanderschaltung von Hauptfeder und Luftreifen fir das Wanken 

€, = Cp e,|(cr + €3)- 
Darin mu die Federkonstante des Luftreifens cg wegen seines gréBeren Abstandes von der Wank- 
achse entsprechend vergréBert werden: 
on = eg(sif)? 


Fiir vier gleichharte Hauptfedern und vier gleichartige Luftreifen mit demselben Innendruck 
wird die Differenz der Krafte in einer kurvenduBeren und in einer kurveninneren Federkombination 
F,— F; =c,fA und das resultierende Rickstellmoment fiir das Wanken ergibt sich wegen 


& Frat, = © 
zu (,f?@—Gh)A= Bi. (3) 
Damit wird aus (2) 
LASERS SS (2a) 


Bei Einzelradfederung sind die Achsen der Rader eines Radpaares oft tber einen Torsions- 
stab verbunden, der ungleichsinnigen Bewegungen der beiden Rader gegeniiber dem Aufbau 
entgegenwirkt und gleichsinnige nicht beeinfluBt. Die Wirkung eines solchen Querstabilisators 
kann in den Riickstellbeiwert des Wankens B einbezogen werden. t 

Wegen der Annahme, daf} die Seitenkrafte der Rader den Radlasten und den Schraglaufwinkeln 
proportional sind, hat die Verlagerung der Radlasten in der Kurve keinen Einflu8 auf die Summe 
der Seitenkrafte und man kann 2S der Gleichung (24) der ersten Arbeit entnehmen: 


rg oe ap 
ES=2h(54—2y +9 bo —69). 
Damit kommt als Bewegungsgleichung des Wankens 
eae 2. . ie: 
—LA= Bath 2k(75—2y +9—0 po —to) al. (2b) 


Darin wurde noch das dampfende Moment dA unter der Voraussetzung einer annahernd geschwin- 
digkeitsproportionalen Dampfung durch die hydraulischen Sto8dampfer eingefiihrt. Fiir vier gleiche 
StoBdampfer mit der Dampfungskonstante k,, die in der gleichen Entfernung f wie die Federn an- 
geordnet sind, wird d = k,f?. 


1S. FuBnote 1, S. 100. 
* Der Einfachheit halber sei vorangesetzt, da die Federn in Teleskophiil fihrt si . 
teile starr mit dem Aufbau verbunden sind.  epeigei acres ‘eee 
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Aus den Gleichungen (25) der ersten Arbeit? fiir das Schieben und Gieren des Fahrzeugs, sowie 
fiir die Lenkschwingung und der Gleichung (2b) fiir das Wanken folgt das System 


ey gee ae Cd oe 4 
C¥+—y—yp—F— w+ 20-0 +nCi=0, | 
= ies z el. 
wCly +—p+2ly—20—y + p(1—2) =0, 


il 4s aes (ne 
Goce) = Ra eae rare ee 


{) oo . Jo Deh Fees 2 
the eA +A + BA) +—y—2yp+p—lC—p—lCp=0. 


Darin wird im letzten Glied der ersten Zeile die Kopplung zwischen seitlichem Schieben und 
Wanken nach (1) beriicksichtigt. Im ersten Glied der zweiten Zeile tritt im dimensionslosen Faktor iv 
die Massenverteilung um die Hochachse auf, wahrend in (25) der ersten Arbeit mit = 1 eine be- 
stimmte Massenverteilung vorausgesetzt war und andere Werte fiir « erst in FuBnoten zu den End- 
formeln beriicksichtigt wurden. 

Um nicht zu unhandliche Ausdriicke zu erhalten, werden in (4) zunachst die Kreiselmomente 
der rollenden Rader weggelassen, die durch das Schwenken der Radachsen infolge des Wankens des 
Aufbaus, des Gierens des gesamten Fahrzeugs und der Lenkschwingung der Vorderrader auftreten. 
Sie werden spater gesondert eingefiihrt werden. Das System (4) fiihrt mit den Ansatzen y = y, e’', 
YF = Yo e°', PY = M e7! und A =A, e** auf die Determinante 


Co+—o betes o—1| ey) h Co? 
eG w Clot +o +20) (1233) 0 | 
me (0) 
1 l I 5 
ao re} pea et? 0 
2 I 1 | 
ee —foot3 a—s rap lieo* + do + B) | 


Die Koeffizienten der aus dem Verschwinden der Determinante (5) folgenden charakteristischen 
Gleichung achten Grades ergeben sich mit I,— mh? = I, und nach Kiirzen durch C/2hk zu 


bf ae re I, 
ao uGl nEG 5 I,, q%=(¥+1)— aE. _ I,+pCl kG ae 

I I l 

$+ ogy @ 
l ep I ae (6) 
= (b+ Dosey Bt ale told +4 |W Cl + Gata Pl + eel oly. 
a df l 

o= CIB +— qt BRE + eed —O)|+LA+O+"g- 7404. 


a= HS) Saath 4+d0to, G=BU+S, “=a =0. 


Darin sind die unterstrichenen Glieder zu vernachlassigen. Unendlich harte Federung und unendlich 
starke Dampfung verhindern das Wanken vollstandig, worauf eine Kontrolle der bisherigen Rechnung 
beruht. Wenn man alle Koeffizienten statt durch C/2 hk durch B/2 hk dividiert und B — oo gehen 
laBt, so verschwinden alle Glieder mit B im Nenner, es werden a) = a, = 0, und a, bis ag gehen 
iiber in die Koeffizienten a, bis a, nach (27) der vorigen Arbeit, wenn man wie dort = 1 setzt. 
Ebenso verschwinden fiir d > 00 a, und a, und a, bis a; gehen iiber in ay bis a, nach (27). 


1S. FuBnote 1, S. 100. 
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Da in (6) die beiden letzten Koeffizienten a, und a, verschwinden, kann die urspriingliche Glei- 
chung vom achten Grade durch o? dividiert werden, so daB sich eine Gleichung sechsten Grades ergibt: 
a, o® + a, 0° + a,o4 + a30° + a0? + a,0 + a, =0. (7) 
Das Verschwinden der beiden letzen Koeffizienten, das auch schon in der vorigen Arbeit beobachtet 
wurde, ist durchaus nicht zufallig. Die daraus entstehende Doppelwurzel 07,, = 0 erméglicht erst 
das Auftreten eines linear mit der Zeit zunehmenden Gliedes in der Lésung fiir eine beliebige Ko- 
ordinate. Dieses letzte Glied stellt den kleinen unmittelbaren Einflu8 der Stérung durch einen 
exzentrischen Sto auf das ganze Fahrzeug dar, der bei der Stabilitatsbetrachtung ausgeschlossen 
bleiben muB. Die Frage wurde in Abschnitt 8 der vorigen Arbeit ausfiihrlicher behandelt. 


5. Kontrolle der Stabilitat nach dem Kriterium von Michailow. Wegen der gréBeren Zahl der a; und 
ihres verwickelten Baues ist es nicht gelungen, mit Hilfe der Hurwitzdeterminanten analog dem in 
der ersten Arbeit beschrittenen Weg die Gleichung einer Grenzflache der Stabilitat in einem Para- 
meterraum aufzustellen. Die Stabilitat jedes Zustandes, der durch eine Anzahl physikalischer Kon- 
stanten bestimmt ist, muB hier gesondert tiberpriift werden. Dazu wird das Kriterium nach Michailow 
verwendet, das im folgenden kurz beschrieben werden soll. Ausfiihrlicher wurde es von M. W. Mejerow* 
behandelt. Abb. 3a zeigt die Lage der drei Wurzeln 0}, 02,03 einer Gleichung dritten Grades in der 
Gaufschen Ebene, die simtlich einen negativen Realteil haben sollen. Man kann die Gleichung als 
Produkt ihrer Linearfaktoren schreiben: 

£(e) = (6 — 94) (6 — 2) (6 — 95) - 

Setzt man fir o die imaginare Zahl iy, so wird 

f(iv) = (iy — 0) 7» — oy) (1¥ — 95), 
wobei die drei Linearfaktoren durch die Pfeile in Abb. 3a dargestellt sind. Der Betrag des Produktes 
dieser komplexen Zahlen ist gleich dem Produkt ihrer Betrage und das zugehérige Argument gleich 
der Summe ihrer Argumente. La8t man nun y von 0 bis + co wachsen, so ergibt sich zunachst bei 
y = 0 fir f(iv) das Argument 0, die Ortskurve beginnt auf der reellen, positiven Halbachse. Fiir 
y =0O wachst das Argument des ersten Linearfaktors um 7/2 und das der beiden anderen fiir die 
konjugiert komplexen Liésungen zusammen um 7, so daB fiir »y = 00 die Funktion f(i v) das Argument 
32/2 und den Betrag 00 hat (Abb. 3b). 


/mag. 


Abb. 3. a. Linearfaktoren einer Abb. 3. b. Ortskurve f(iv) einer Abb. 4. Ortskurve f(i) einer Gleichung sechsten Grad 
Gleichung dritten Grades. Gleichung dritten Grades. deren seatiice Woes negative Tesitedle iabane 


Allgemein mufs das Argument von f(iy) stetig von 0 auf nz/2 wachsen, wenn y von 0 bis + 00 
wachst und n der Grad der Gleichung ist. Das ist nur méglich, wenn alle Wurzeln links von der ima- 
gindren Achse liegen, so daB der Verlauf der Kurve f(iv) iiber die Stabilitat des durch die charak- 
teristische Gleichung f(c) = 0 bestimmten Systems entscheidet. ; 

In unserem Fall ist die charakteristische Gleichung (7) vom sechsten Grade. Wenn man fiir o 
wieder iy setzt, so wird aus (7) os 


fir) =— ay + ia, + ayy! —iasr?—a,v?+ia,v+a,=0. (8) 
Es ist nun nicht nétig, diese Kurve in der komplexen Ebene aufzuzeichnen, da man aus der Lage 
ihrer Schnittpunkte mit der reellen Achse bereits auf das Wesentliche ihres Verlaufes schlieBen kann. 
Das Nullsetzen des Imaginarteils von (8) fiihrt nach Division durch » ‘auf die biquadratische Glei- 
chung a,74— agv? +a, = 0. Damit wird 
ops ss Vai —4 a, a, 
= V1.2 ac a, . (9) 


1 M. W. Mejerow, Grundlagen der selbsttatigen Regelung elektrischer Maschinen, S. 64, Berlin 1954. 
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Der Realteil von (8) muB nun fir y = 0 positiv, fiir die kleinere Frequenzy, negativ und fir die 
- gréBere Frequenz y, wieder positiv werden, wenn f(iv) mit zunehmendem » die Quadranten der Reihe 


nach durchlaufen soll (Abb. 4). 
Da wegen a, > 0 die erste Bedingung stets erfiillt ist, bleiben als Stabilitatsbedingungen: 
— 4,0} + a,r{ —a,¥? + a, > 0,) 
— Ay V3 + A293 — a4 95 bag<0,f 
worin v7 bzw. v3 der Gleichung (9) zu entnehmen ist. 

Auf diese Weise la8t sich der durch die charakteristische Gleichung (7) mit den Koeffizienten 
nach (6) gegebene Zustand des Kraftwagens ziemlich rasch auf Stabilitat untersuchen. Ein gefederter 
Wagen wird bei niedrigerer Geschwindigkeit unstabil als ein ungefederter; doch ist die Abnahme 
der Stabilitat bei normaler Federung und Dampfung nicht allzu groB, solange man nicht den im 


ganzen sehr schadlichen Einflu8 der Kreiselmomente der Rader in Rechnung stellt. Bevor das ge- 
schehen kann, miissen noch einige Eigenheiten der Kraftwagenfederung behandelt werden. 


(10) 


6. Hindernisfederung und Kurvenfederung. Die obere Reihe der Skizzen in Abb. 5 zeigt das Ver- 
halten eines Rades beim Uberfahren eines einseitigen Hindernisses fiir verschiedene Radaufhangungen. 
Man sieht, daf bei Parallelfiihrung die Radachse nur verschoben, aber nicht hochgeschwenkt wird, 
so da8 kein Kreiselmoment geweckt wird. Dagegen entstehen bei den Pendelachsen starke Kreisel- 
momente, die um so gréfer sind, je mehr die Pendelachsen verkiirzt werden. Beim Uberfahren einer 
einseitigen Unebenheit ist in diesen Fallen die Stérung durch das Kreiselmoment des Vorderrades 
viel starker als die durch das Moment der horizontalen Komponente der Sto®kraft um den Lenk- 
zapfen, von der bisher ausgegangen wurde. Unter diesen Umstanden hilft es nicht viel, daB Kreisel- 
moment und Moment der StoBkraft entgegengesetzten Drehsinn haben. Wie aus Abb. 5c ersichtlich, 


sucht das Kreiselmoment I,, A, das linke Vorderrad des gegen die Bildebene fahrenden Wagens zur 
Wagenmitte einzuschlagen, wahrend das Sto{moment das Rad nach aufen einzuschlagen sucht. 
Deshalb hat schon de Lavaud' vorgeschlagen, bei starren Vorderachsen den Abstand des Lenkzapfens 
vom Vorderrad zu vergréBern und so das StoSmoment dem Kreiselmoment anzugleichen. Dabei 
ist jedoch zu bedenken, da die beiden Momente nicht gleichzeitig ihr Maximum erreichen, das beim 
Kreiselmoment mit dem gréBten es beim StoSmoment mit der gréBten Horizontalkomponente der 
Bodenkraft auftritt. Die Vorgange beim Uberfahren eines Einzelhindernisses, die fiir den zeitlichen 
Verlauf der beiden Momente von Bedeutung sind, hat E. Marquard? behandelt. 


Paralle/fuhrung ah Pende/achsen verkurzte Pendelachsen 


federung 


Ap/A negativ 


Abb. 5. Hindernistederung und Kurvenfederung bei verschiedenen Radaufhangungen. 


Die Entwicklung ist seit de Lavaud in Richtung der ebenfalls von ihm empfohlenen Kinzelauf- 
hangung der Vorderrader verlaufen und hat zu den sog. ,,flattersicheren“ Parallelfihrungen der Vorder- 
rader gefiihrt, die heute bei den Personenautos vorherrschen, wahrend die Lastwagen meist bei der 
Starrachse geblieben sind. Die Parallelfiihrung hat namlich beim Durchfahren einer Kurve den 
schweren Nachteil, daB die Rader die Neigung des Aufbaus voll mitmachen miissen (Abb. 5d), was 


“1S. de Lavaud, Causerie scientifique, Paris 1927 (im Selbstverlag des Verfassers). 
2 E. Marquard, Schwingungsdynamik des schnellen StraBenfahrzeugs, Essen 1952. 
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zu starken Kreiselmomenten der Vorderrader um ihre Lenkzapfen fihrt. Diese Kreiselmomente 
bilden bei hohen Geschwindigkeiten neben den seitlichen WindstoBen die Hauptgefahr fiir die Sta- 
bilitat, wie sich in den folgenden Abschnitten zeigen wird. Bei den Personenwagen kann man ihren 
Einflu8 durch leichte Rader und tiefe Schwerpunktslage verringern, bei den Lastwagen ist beides 
nicht méglich. Deswegen hat sich gerade bei schnellfahrenden Lastwagen die starre Vorderachse 
gehalten, die nur eine kleine Kurvenneigung der Rader infolge der Reifeneindriickung erlaubt. Auch 
die Eingelenkpendelachsen (Abb. 5e) zeigen in der Kurve ein ahnlich giinstiges Verhalten, das sich 
erst bei starkerer Verkiirzung der Pendelachsen wieder verschlechtert. Ungiinstig ist bei den Pendel- 
achsen, da8 durch Belastungsinderung oder durch Uberfahren einer Bodenwelle starke Anderungen 
der Spurweite und des Radsturzes! (letzteres bedeutet praktisch groBe Ap) entstehen. 


7. Die Kreiselmomente der Rader. Infolge der Gier- und Lenkschwingung um zwei praktisch 
parallele Achsen und des Wankens um eine dazu annahernd normale Achse iiben die rotierenden 
Rader Kreiselmomente auf ihre Lagerung aus, durch die diese drei Freiheitsgrade der Bewegung 
gyroskopisch gekoppelt werden. Es entstehen die Momente: 


durch das Kippen der Hinterrader 2 T,—Ag um die Hochachse, 
durch das Kippen eines Vorderrades —I, Ba um den Lenkzapfen?, 

; (11) 
durch Gieren aller vier Rader —4I, - y 


um die Wankachse. 
und durch die Lenkschwingung der Vorderrader 21 o~ yp 


Fir kleinere Schwingungen kann man den Kippwinkel/, der Rader der Seitenneigung des Aufbaus A 
proportional setzen. Damit wird Ap =A(Ap/A). Die Konstante Ap/A = MM,/AMp (Abb. 1); sie 
ist positiv, wenn A und M, auf verschiedenen Seiten des Fahrzeugs liegen, andernfalls ist sie negativ, 
d.h. Aufbau und Rader neigen sich in der Kurve nach verschiedenen Seiten (Abb. 5f). Durch die 
verschiedene Kindriickung der Luftreifen kommt zu dem so berechneten Wert noch ein kleiner po- 
sitiver Betrag hinzu, der fiir Parallelfiihrung verschwindet. Im folgenden Schema sind die Elemente 
zusammengestellt, durch die die Determinante (5) infolge der Kreiselmomente erganzt werden muB: 


as Side Ee aie ere ar tet (Fle 
1 (12) 


rE Po 9| Ee ee) 
Tr: 


2hk 


Den Einflu8 der verschiedenen Kreiselmomente sieht man aus der GréBe der Faktoren vor den 
eckigen Klammern. Danach ist das Element in der zweiten Zeile bedeutungslos, das von den Kreisel- 
momenten der Hinterrader infolge des Wankens stammt; es wird im weiteren Rechnungsgang ver- 
nachlassigt. Das Element in der dritten Zeile, das das Kreiselmoment eines Vorderrades infolge des 
Wankens enthilt, ist das weitaus wichtigste. Obwohl die Kreiselmomente der Vorderrader und der 
Hinterrader nahezu gleich sind, sind die ersteren, die die Lenkung beeinflussen unvergleichlich wirk- 
samer als die letzteren, die direkt auf das Fahrzeug wirken. Die Kreiselmomente der Vorderrader 
infolge des Wankens sind der Stabilitét sehr abtraglich und miissen stets beriicksichtigt werden. 
Die Elemente der vierten Zeile stammen von den Kreiselmomenten von Vorder- und Hinterradern 
infolge der Gier- und Lenkschwingungen. Diese sind im ganzen fiir die Stabilitat giinstig, aber von 
- geringerem Einflu8 als das Element in der dritten Zeile. Sie sind hier als direkt am Aufbau angreifend 
angenommen, wahrend sie tatsachlich an der vernachlassigten Zwischenmasse angreifen und Haupt- 
feder und Reifen im entgegengesetzten Sinn belasten. 


_ 1 Der Sturzwinkel, den die Radebene mit einer Normalebene zur Fahrbahn mit gemeinsamer Spur auf dieser 
einschlieBt, hat einen kleinen, hier vernachlassigten Einflu® auf die vom Rad aufgenommene Seitenkraft. 


* Der positive Drehsinn fiir ~ wurde dem fiir y entgegengesetzt angenommen (s. FuB8note 1, Seite 100, 
dort Abb. 3). 
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ie deen aa ean der Determinante (5) durch die Kreiselglieder nach dem Schema (12) 
ganzt und die Koeffiziente h ie : 
Sean n neu errechnet, so kommen zu den Koeffizienten a, bis a; nach (6) fol- 


at 1G: Ap v2 

ge a pala . 
ee TAs ieee 1) TEE RY 
2 eo AYU 40 (AR 

2 ie where (R\ + meray (a) 
A a 2 AT /Ap\ v2 (13) 

ws | R w R 

a4, iz ke 0) (B) + mEG =o (tle? 


y= 4h SATS + cae +9 (BS. 


Darin sind mehrere bedeutungslose Glieder bereits weggelassen. Die nicht unterstrichenen Glieder 
entstehen ausschlieBlich aus den Kreiselmomenten der Vorderrader um ihre Lenkzapfen infolge des 
Wankens und miissen stets beriicksichtigt werden. Die restlichen Glieder ergeben sich ie den 
Kreiselmomenten aller Rader um die Wankachse infolge des Gierens des Fahrzeugs und der Lenk- 
schwingung. Sie miissen bei sehr hohen Geschwindigkeiten (etwa ab 150 km/h) beachtet werden. 

Die Kreiselwirkung der Rader soll so klein wie méglich gehalten werden. Auf den ersten Blick 
scheint es am besten, den in allen Gliedern vorkommenden Faktor (Ap/A) zum Verschwinden zu 
bringen, d.h. die Vorderrader so aufzuhangen, daB sie das seitliche Wanken des ‘Aufbaus in der 
Kurve nicht mitmachen. In Abschnitt 6 wurde gezeigt, daB eine solche Radaufhangung durch 
Kingelenkpendelachsen beim Uberfahren eines Hindernisses zum Hochschwenken der betreffenden 
Achse und damit zu einer Stérung durch ein Kreiselmoment fiihrt. Daher sind bei den Vorderradern 
die Pendelachsen ziemlich selten und man findet sie nur bei einigen Kleinstwagen neuerer Bauart. 
Die meisten Kraftwagen haben Parallelfiihrung der Vorderrader mit Ap/A = 1. Sie vermeiden da- 
durch zwar die Stérung durch das Kreiselmoment eines Vorderrades beim Uberfahren einer Uneben- 
heit, bringen aber im Ablauf der Bewegung des durch das unvermeidliche StoSmoment um den Lenk- 
zapfen gestérten Fahrzeugs die Kreiselmomente vor allem der Vorderrader voll ins Spiel. 

Es ist deshalb zu erwagen, ob es nicht besser ware, beim Uberfahren einer Unebenheit ein kleines 
Kreiselmoment zuzulassen, das ja dem aus dem gleichen AnlaB auftretenden StoBmoment entgegen- 
gerichtet ware. Das letztere kann dadurch vergréBert werden, daB der Lenkzapfen keine Spreizung 
erhalt. Das wiirde zu einer Radaufhangung durch Querlenker fiihren, die nicht parallel sind, sondern 
gegen die Wagenmitte zu zusammenlaufen, wobei der Schnitt- 
punkt ihrer Mittellinien M, bei ungestérter Fahrt etwa in der 
Ebene des gegeniiberliegenden Rades liegt (Abb. 6). Dann ware 
mit Ap/A ~ 0,5 die Wirkung der Kreiselmomente auf den Ablauf 
der Bewegung stark herabgesetzt. Auch die Hebung des Momentan- 
pols M, und die sich daraus ergebende Verkiirzung des Hebelarms 
der Fliehkraft h wirken im gleichen Sinne. Die Hinterrader kénnten 


an Eingelenkpendelachsen mit tiefliegendem Drehpunkt M, ge- bb. 6 si Se aan 
MItA R, = U,0- 


lagert sein, so daB die Wankachse M,M, nahezu horizontal verlauft. 

Bei der iiblichen Parallelfiihrung der Vorderrader liegt M, etwa in Fahrbahnhéhe. Um den Hebels- 
arm der Fliehkraft zu verringern, wird M, durch die Verwendung verkiirzter Pendelachsen fiir die 
Hinterrader hoch gelegt, wodurch die Wankachse in einem Winkel bis 10° gegen die Fahrbahn von 
vorne nach hinten ansteigt. 

Aus den Ausdriicken (13) sieht man, da® die Wirkung der Kreiselmomente der Rader auf das 
bereits gestérte Fahrzeug vor allen von Ag/A, h und der Drehmasse I, eines Rades um seine Achse 
abhangt, die alle in den urspriinglichen Koeffizienten nach (6) nicht vorkommen. Zumindest eine 
dieser drei GréBen mu8 so klein wie méglich gehalten werden, wenn nicht bereits bei erreichbaren 
Fahrgeschwindigkeiten Unstabilitat des Fahrzeugs eintreten soll. 


8. Uberpriifung der Lenkstabilitat eines langsamen Personenkraftwagens. In diesem und den folgen- 
den Abschnitten wird die Lenkstabilitat von drei ausgefiihrten Kraftwagen bei der jeweils angegebenen 
Héchstgeschwindigkeit tiberpriift. Die Wagen unterscheiden sich in Bauart und Verwendungszweck 
grundlegend voneinander. Alle Daten, die von den Herstellerwerken der Fachliteratur tibergeben 
wurden, sind in der Berechnung verwendet. Leider geniigen diese Angaben keineswegs. Der Verfasser 
war bemiiht, den restlichen zur Berechnung notwendigen Gré8en technisch mogliche Werte zu geben. 

. i 


108 G. Mitterlehner: Der Einflu8 der Kraftwagenfederung auf die Lenkstabilitat Ingenieur-Archiv 


Das erste Beispiel behandelt einen sehr leichten, aber geraumigen Personenwagen mit Kurbel- 
achsen vorne und hinten, was etwa einer Parallelfiihrung der Rader durch horizontale Querlenker 
entspricht (A/a = 1). Die Wankachse liegt in diesem Fall in Fahrbahnhohe, so daB als Hebelsarm 
der Fliehkraft die ganze Schwerpunktshéhe iiber der Fahrbahn auftritt. Von einem solchen F ahrzeug 
kann man keine Stabilitat bei hohen Geschwindigkeiten verlangen. In der Tat ist es mit einem fiir 
seine GréBe ausnehmend schwachen Motor ausgeriistet, mit dem es nur eine Héchstgeschwindigkeit 
von rd. 80 km/h erreichen kann. Vor allem die geringe axiale Drehmasse der verwendeten Rader 
bewirkt, daB das Fahrzeug bei dieser Geschwindigkeit noch stabil bleibt. Bei der Wahl des Ruck- 
stellbeiwertes des Wankens wurde die auffallend weiche Federung beriicksichtigt. 


Angaben: Fahrzeuggewicht G = 495 kp (mit Fahrer G = 600 kp), Radstand | = 2,4m, Spur- 
weite s = 1,26 m, Héchstgeschwindigkeit v,,,, = 80 km/h. 


Annahmen: 
Wirksamer Hebelsarm der Fliehkraft Osos 
Drehmasse des Aufbaus um die Langsschwerachse I, = 30 kpms?, 
Drehmasse eines Rades um den Lenkzapfen I, = 0,02 kpms?, 
Drehmasse eines Rades um den Achsschenkel I,, = 0,03 kpms?, 
wirksamer dynamischer Halbmesser eines Rades r-= 053hms 
wirksamer Nachlauf der Vorderrader n = 0,02 m, 
Federkonstante einer Hauptfeder cp = 4700 kp/m, 
Federkonstante eines Luftreifens Cr = 10000 kp/m, 
Federarm/Radarm (Abb. 5a) Lil, = 0,6, 
Dampfungskonstante eines Sto8dampfers ik, = 1% kpws, 
Beiwert der Massenverteilung um die Hochachse Tae 
Beiwert der Radlastanderung durch den Luftwiderstand C2015 10F a 
Seitensteifigkeit eines Rades unter der statischen Radlast k = 1200 kp/rad. 


Damit werden die Federkonstante einer Hauptfeder, bezogen auf die Spurweite 
Cy, = Cp (I/ plz)? = 1690 kp/m, 
die resultierende Federkonstante jedes Rades 
ta 
C= oe. 1450 kp/m , 
und der Riickstellbeiwert des Wankens 
B, = ¢,s°?—Gh = 2000 kpm. 
Weiter ergeben sich mit k} = ky (Ip/lg)? = 63 kpm—s der Dampfungsbeiwert des Wankens zu 
Fl aes ky s* = 100 kpms, 
der Tragheitsbeiwert zu C = m/4 k = 0,0125 m7!s?, 
die RechengréBe I) zu I) = I, — mh? = 15 kpms? 
sowie fiir v = 22,2 m/s (= 80 km/h) der Beiwert der Radlastanderung durch den Luftwiderstand zu 
€ = 0,0369. 
Damit kommen die Koeffizienten nach (6) samt ihren Erganzungen infolge der Kreiselmomente 
nach (13), wenn jeweils die unterstrichenen Glieder weggelassen werden, 
ao = 9,780 -10, a, = 0,820-10, a,=—0,970, a= 9,13, a, =114,8, a, = 558, a, = 2074- 
(Die a; sind nicht dimensionslos, a, hat die Einheit kpm, a; kpms usw.). 
Nach (9) wird 
vi = 1048, 12 = 64,8 
und nach (10) 


— Go} + ari — av? + a, = 5,0-10!> 0, 
— a, 8 + a. v5 — ay v2 -+ ag = —15,0-10°< 0. 
Der Zustand des Fahrzeugs ist demnach bei v = 80 km/h stabil. Wenn der Rechnungsgang fiir eine 


Geschwindigkeit von 100 km/h wiederholt wird, was wegen ¢ = C(v) eine Neuberechnung sadmtlicher 
a; nétig macht, ergibt sich bereits Unstabilitat. 


9. Darstellung des Einflusses von Parametern. Der Einflu8 eines Berge cree tT, der nur in der Art 
Qiv) +7 Riv) =0 © (14) 
in die Gleichung (8) f(iv) = 0 eingehen muB, kann mit 

t= ing Qu ») 
R(iv) : 
in der Ebene der komplexen Zahlen dargestellt werden. 
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2 sae man = iia die Drehmasse eines Vorderrades um den Achsschenkel + = I,,, was 
nicht zu grofe Geschwindigkeiten zulassig ist, da dann die unterstrichenen Glieder in (13) weg- 
gelassen werden kénnen, so sind in " 
6 
ir) = ae 
n=0 : 


und 
6 
Riir) = S a,.0°—" 
n=0 


mit den Werten des obigen Beispiels die Koeffizienten 


y A. (=a wr = = 
tia 0,0 0° 13.0 sets 2,9 =o, 4,9 = % 
a2 = 0, i 0, a. = 0, a,7 = 0 
ae Gy: 9 = 0005 -saz 5 102;3 @; 9 = 328, 
a,,=111l, a,=415 a; , = 7670 
elnzusetzen. 


Die Kurve 7t(i ) schneidet die reelle Achse ein- oder mehrmals und legt so die Grenzwerte fest, 
innerhalb der man 17 (= I,,) verandern kann, ohne da® der Zustand des Fahrzeugs unstabil wird. 
Wo der stabile Bereich liegt, mu8 man durch Einsetzen eines Wertes fiir I,, in dem betreffenden 
Bereich nach dem Michailow-Kriterium ermitteln. 

Die Kurve 7(i y) fiir v = 0 bis —oo ist das Spiegelbild der Kurve 7(iy) fiir y = 0 bis + 00 be- 
zuglich der reellen Achse. Die Punkte auf den Kurven 7(i v) sind in Abb. 7 und Abb. 8 mit dem Wert 
der Frequenz  bezeichnet, fiir den sie errechnet wurden. Die beiden Aste von 7(i y) zerlegen in Abb. 7 
die komplexe Ebene in mehrere Gebiete. Wenn man die Begrenzung des stabilen Gebiets innerhalb 
desselben durch Schraffur 
auszeichnet und _ diese Z m2 (io) 
Schraffur auch dort fort- 
setzt, wo die Kurven nicht 
mehr die Grenze des sta- 
bilen Gebiets bilden, kann eR 
man die sog. Schraffur- \ * 
regel von Neimark} an- \ 

wenden: ‘ 
" 


Bei jedem Durchgang 
durch die Kurve von 
der schraffierten Seite her 7 R 
geht eine Wurzel der cha- Ly j a 

| 


A 005 Q7kpm sek? 
INS 
Rechnungswert 


aon : . 
_rakteristischen Gleichung 


mit negativem Realteil de 
in eine Wurzel mit posi- 
tivem Realteil iiber. Z. B. Sorel 
kommt man bei Vergré- Abb. 7. EinfluB der axialen Drehmasse 
Bern des Parameters Te eines Vorderrades auf die Stabilitat. 
der nur reeller positiver 4 : 
Werte fahig ist, aus dem “7? saGie Hachachse auf die Stabilitat. 
stabilen Gebiet I gleich 
in das Gebiet III mit zwei Wurzeln mit positivem Realteil, die natiirlich konjugiert komplex sind. 
An Stelle von lauter gedampften Frequenzen tritt jenseits des Grenzwertes von I, auf der reellen Achse 
eine angefachte Frequenz auf, wahrend die iibrigen (zunachst) gedampft bleiben. In Abb. 7 ist auch der 
ausgefilhrte Wert von I,, eingetragen. Man sieht, daB schon eine verhaltnismaBig geringe VergréBerung 
der axialen Drehmasse der Rader den Fahrzustand bei der Héchstgeschwindigkeit unstabil machen 
wiirde. Der Verlauf der Funktion 7(i v) links von der imaginaren Achse ist praktisch bedeutungslos. 
Wenn man als Parameter statt I,, das Verhaltnis des Kippwinkels der Vorderrader zur Seiten- 
neigung des Aufbaus /p/A wablt, erhalt man einen dbnlichen Verlauf von T(iv) wie in Abb. 7, da 
Ag/4 in denselben Koeffizienten auftritt wie I. Nur der MaBstab ist dann ein anderer. Der wirk- 
same Hebelsarm der Fliehkraft h kann nicht als Parameter gewahlt werden, da er wegen I, = I, — mh? 


A, 
coo-s_ 


stabiler 
Bereich 


1S. M. W. Mejerow, Grundlagen der selbsttatigen Regelung elektrischer Maschinen, S. 68, Berlin 1954. 
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nicht in der Art (14) in die Gleichung (8) eingeht. Davon abgesehen, gilt fiir h dasselbe wie fiir 
I, und Ap/A. 


Dagegen kann der Einflu8 der Massenverteilung um die Hochachse bei gegebenem Radstand 
: ml? 
und gegebener Gesamtmasse dargestellt werden. Es ist dann t = = sb. n 


Mit den Werten des obigen Beispiels, bei dem 4 = 1 gesetzt wurde, werden 


Ayo = 0, a, 9 =0,00281, aay HO, go = 5,94, My = 54,8, 5,9 = 440, 

dy, = A, a, , = 0,005397 a, = 0,975. a, = 3,19, ag, == 000, a5, = Ue 
sowie 

Ge.0 = Se» de, = 0. 


In diesem Fall ist der Nenner von 7(i), d. i. die Funktion R(i 1) in (14) vom sechsten Grade und 
7(i v) hat deshalb in der komplexen Ebene einen ziemlich verwickelten Verlauf (s. Abb. 8). Fury = 0 
ist T(iv) = —oo und fiir y-> 0o geht T(iv) > 0. Die mit unrunden Zahlen bezeichneten Punkte 
der Kurve sind fiir die Nullstellen des Realteils oder des Imaginarteils des Nenners errechnet. 

Da der Parameter yz wieder nur reelle positive Werte annehmen kann, sind nur die Schnittpunkte 
der Kurve mit der reellen positiven Achse interessant. Die Grenze des stabilen Bereichs liegt hier 
bei = 0,85. Eine VergréBerung des Faktors y, d. h, der Drehmasse des Fahrzeugs um die Hoch- 
achse, bei gleichbleibendem Gewicht und Radstand wirkt stabilisierend. Diese Folgerung gilt ganz 
allgemein, da der Charakter der Funktion 7(iy) auch fiir andere Fahrzeugdaten erhalten bleibt. Sie 
stimmt auch mit den Erfahrungen iberein, die mit Rennwagen gemacht wurden. 


10. Uberpriifung der Lenkstabilitat eines schnellen Sportwagens. Die fiir dieses Fahrzeug der Lite- 
ratur entnommenen Angaben sind: Gewicht G = 1300 kp, Radstand 1 = 2,4 m, Spurweite s = 1,4m 
(Mittelwert) und Héchstgeschwindigkeit v,,,, = 240 km/h (66,7 m/s). Es handelt sich also um einen 
sehr schnellen Kraftwagen, weshalb die Koeffizienten der charakteristischen Gleichung nach (6) 
durch die vollstandigen, von den Kreiselmomenten stam- 
menden Ausdriicke (13) ergénzt werden miissen. Abb. 9 


lenkern, wahrend rechts die Hinterradaufhangung durch 
Eingelenkpendelachsen mit tiefliegendem Drehpunkt zu 
sehen ist. Beides ist an dem betreffenden Fahrzeug ausge- 
fiihrt, doch sind die MaBe in Abb.9 angenommen. Die 
Momentanpole des Wankens ftir die beiden verschiedenen 
Nae OV order usd diintat chatenacbuoe Radaufhangungen sind M, und M,. Bei einer Schwerpunkts- 
eines schnellen Sportwagens. lage in der Mitte des Radstandes liegt der Schnittpunkt der 
Wankachse M,M, mit der Querschnittsebene durch © in 
Abb. 9 in der Mitte zwischen M, und M,. Damit wird h = 0,275 m. Nachdem /,/A = 1 ist und die 
Rader bei einem so schnellen Fahrzeug eher gréBer als normal sein miissen, so daB man I, nicht 
sehr klein halten kann, tragt ein kleines h am meisten zur Stabilisierung bei. 
Fiir die Rader wurde gewahlt: r = 0,324m, 
I, = 0,045 kpms? , I, = 0,035 kpms?, n = 0,02 m, k = 1500 kp/rad. 


Wenn man bedenkt, daB das schon in der ersten Arbeit behandelte ungefederte Fahrzeug eine 
besonders hohe Lenkstabilitat aufweist, wird man fiir sehr schnelle Fahrzeuge eine ziemlich harte 


zeigt links die Vorderradaufhangung an parallelen Quer- | ii 


Federung wahlen. Mit cz = 20000 kp/m und c7, = 8000 kp/m, wobei der letztere Wert fiir f=s5 | 
gilt, wird c, = 5720 kp/m, und eS 10850 kpm. Mit einer spezifischen Dampfung der Relativbewe- | 


‘“Cp= 
4 F 

= 255 kpm™s und d=k)s?=500kpms. Fiir die Parallelfihrung vorne sind c+ (Ip/lp)®, baw. 
k aR|Ee)? (Abb. 5a) die Konstanten einer Hauptfeder, bzw. eines StoBdampfers, wahrend diese Werte 
tir die Pendelachsen mit (s/f)? multipliziert werden miissen. Die Drehmasse des Aufbaus um seine 


Langsschwerachse wurde zunachst mit I,, = 75 kpms? angenommen. Obwohl hier die Wankachse 
M,M, zur Langsschwerachse nicht genau parallel ist, gilt fir die RechengréBe I, annahernd 

o = [,— mh’ = 65,2 kpms®. Als Beiwert der Massenverteilung um die Hochachse wurde ee) 
gewahlt, was bedeutet, daB man die F ahrzeugmasse in 2 Punkten zwischen den Radpaaren konzentriert 
denken kann. Der Tragheitsheiwert ergibt sich zu C = m/4k = 0,02165 m—s?. 


1 L. Kraus, Automobiltechn. Z. 59 (1957) S. 122. 


gung zwischen Achse und Aufbau bei reiner Hubschwingung D = 0,25 wird ke == 0:25. 2)/% 


XXIX. J i - c 
Band 1960 G, Mitterlehner: Der EinfluB der Kraftwagenfederung auf die Lenkstabilitat ee 


Die Luftwiderstandszahl k; wurde aus der Bedin 
yee gung errechnet, da der Luftwiderstand d 
iat des Rollwiderstandes bei einer Fahrgeschwindigkeit von etwas iiber 60 km/h ee a 
it dem Rollwiderstandsbeiwert fr = 9,02 wird dann k; = 0,09 kpm-*s?. Die durch den Luft- 
widerstand bewirkte bezogene Radlastanderung wird 
BePAN: "9-hok 
é == = EE vt = 0,26 - 10-1 o?, 


worin h; = 0,5 m die Héhe der Wirkungslinie des Luftwiderstandes iiber d i u 
po 66,7 m/s, das sind 240 km/h, wird € = 0,155. aap dew 
a Mit den oben angegebenen Werten wurden die Koeffizienten der charakteristischen Gleichung a, 
is @, nach (6) errechnet, wobei a, bis a, um die Glieder infolge der Kreiselmomente ay, bis a;;, nach (13) 
6 


vergroBert wurden. Damit ist die Funktion (8) f(iv) = > 4,(i v)°—" festgelegt. Ihr Verlauf in der 
0 


komplexen Ebene muB fiir 0 <_» < oo durch Bestimmung ihrer Schnittpunkte mit der reellen Achse 
untersucht werden. Aus dem Nullsetzen des Imaginarteils von (8) ergeben sich aus (9) die Werte 
von », bei denen f(i7) die reelle Achse schneidet. Damit f(i) die Quadranten der komplexen Ebene 
der Reihe nach durchlauft, wie es zur Stabilitat notwendig ist, mu8 nach (10) fiir den gréBeren Wert » 
der Wert if (y,) > 0 und fiir den Kleineren Wert y, der Wert f(iv.) < 0 sein. 

Wie sich aus der ersten Zeile von Tabelle 1 ergibt, ist das System fiir die zunachst gewahlten 
Werte der Parameter unstabil. Es wurde deshalb zuerst die Drehmasse J,, eines Vorderrades um 
semen Lenkzapfen vergréfert (Zeile 2), dann auBerdem die Dampfung vergroBert (Zeile 3) und weiters 
zusatzlich die Drehmasse um die Langsachse I, verringert (Zeile 4). SchlieBlich wurde mit dem 
Faktor 4 noch die Drehmasse um die Hochachse vergréBert (Zeile 5) und endlich die Drehmasse 
um die Langsachse abermals vermindert (Zeile 6). Im Ergebnis dieser Zeile wurden erstmals sémtliche 
Kreiselglieder nach (13) beriicksichtigt, wahrend vorher die kleinen zweiten Glieder rechts in ag, 
und a,, weggelassen wurden. Die zweimalige Verkleinerung von I, und die VergréBerung von d 


dienen offenbar zur verstarkten Dampfung des Wankens. Dieser Freiheitsgrad soll — nach 
Weglassen der Koppelglieder — nahezu aperiodisch gedampft sein. 
3 Tabelle 1. 

i, kpms? | dkpms | I. kpms? uw (1) ” Sa ” Su* fi») f (is) | Zustand 
0,035 500 15 | 1 955 38,8 —70,2-104 | — 71,7-10?| unstabil 
0,070 500 715 | 1 | 469 39,9 —12,6-104 | — 65,5: 10?| unstabil 
0,070 7150 15) | 1 430 36,4 — 5,6-104 | — 58,2-102)| unstabil 
0,070 750 60 1 437 35,4 — 1,9-104 | — 57,7: 10? | unstabil 
0,070 7150 60 | 1,5 409 28,0 — 0,1-104 | — 92,2 - 10? ? 

0,070 750 50 | 1,5 416 29,0 + 0,7: 104 | —102- 10? stabil 


Man sieht, daB die Werte von mehreren Parametern erheblich geandert werden muften, um das 
Fahrzeug in einen stabilen Zustand zu bringen. Die Anderung zahlreicher anderer Parameter inner- 
halb technisch méglicher Grenzen hatte geringen EinfluB auf die Stabilitat. Andere Parameter mit 
starkem EinfluB, wie z. B. h, waren von vornherein mit einem recht gimstigen Wert angenommen, 
der bei der gegebenen Art der Radaufhangung kaum verbessert werden konnte, oder waren festgelegt, 
wie Gewicht, Spurweite, Radstand und Héchstgeschwindigkeit. 

Die Stabilitat des hier behandelten linearen Ersatzsystems fiir den Kraftwagen gegeniiber kleinen, 
von der Fahrbahn ausgehenden Stérungen ist bei hohen Geschwindigkeiten jedenfalls nicht leicht 
oder u. U. auch garnicht zu erreichen. Dennoch sollte beim Entwurf eines Kraftwagens zunichst 
diese ,,lineare Lenkstabilitat“ angestrebt und die herkémmlichen Stabilisierungsmittel wie progressive 
Federung, Vorspur etc. erst zusatzlich angewendet werden. 

Bei den Vorderradern kann die Masse der nicht rotierenden Teile, wie Bremsbacken samt La- 
gerung, erhéht werden, um die Drehmasse I, des Vorderrades um seinen Lenkzapfen gegeniiber der 
Drehmasse um den Achsschenkel J,, zu vergréBern. Das kann wirksamer auch durch VergréBern 
des Abstandes des Lenkzapfens von der Radebene geschehen. J[,,, das in allen Kreiselgliedern auf- 
tritt, soll dagegen so niedrig wie méglich gehalten werden. Durch VergréBern von I, und eventuell 
durch Verkleinern von k wird die Nachgiebigkeit der Lenkung vermehrt und der Nachgiebigkeit 
der Federung angepaBt. Diese ist keine feste GréBe, sondern wachst mit der auslenkenden Fliehkraft 
und damit mit der Fahrgeschwindigkeit. Es mu8 in diesem Zusammenhang erwahnt werden, daB 
de Lavaud schon anlaBlich des Uberganges von Hochdruck- auf Niederdruckreifen nach dem ersten 
Weltkrieg die Ansicht vertreten hat, daB die Einfiihrung einer weicheren Federung eine mehr nach- 
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giebige Lenkung bedingt. Das ergab sich aus der Untersuchung des von ihm verwendeten Systems 
einer starren Vorderachse mit um die Lenkzapfen schwenkbaren rollenden Radern und sollte auch 
heute bei der Einfiihrung der hydraulischen Lenkhilfen beachtet werden. Ks ist denkbar, daB eine 
solche Servolenkung gegeniiber einer von der Fahrbahn ausgehenden Stérung zu wenig elastisch ist 
und deshalb der Zustand des Fahrzeugs bei hohen Geschwindigkeiten unstabil wird. Eine vollstandige 
Fixierung der Achsschenkel in der Geradeaus-Stellung gegeniiber StoBen der Fahrbabn verhindert 
zwar alle periodischen Bewegungen des Fahrzeugs, der Lenker verliert aber das »Strabengefithl . 
weshalb die Korrektur der auch dann auftretenden kleinen Seitenschiebung und des Gierwinkels 
erst dann erfolgen kann, wenn die Kursabweichung sichtbar wird. 


11. Der Momentanpol des Wankens und der Kippwinkel der Rader bei der Starrachse. Es wurde 
bereits erwahnt, daB bei Lastkraftwagen vielfach Starrachsen auch zur Vorderradaufhangung ver- 
wendet werden. In diesem Abschnitt sollen die Vor- und Nachteile der Starrachse fiir die Stabilitat: 
beschrieben werden. 

Zunachst werden beim Uberfahren eines einseitigen Hindernisses die Achsen beider Rader um 
den gleichen Winkel verschwenkt, so da8 zu Beginn der Stérung neben dem StoSmoment auch ein 
Kreiselmoment auftritt. Ein weiterer Nachteil ist, da bei der Trampelschwingung der Starrachse 
samt Radern zwischen Hauptfedern und Luftreifen und der Lenkschwingung der Rader, die mit- 
einander gyroskopisch gekoppelt sind, recht unangenehme Resonanzerscheinungen eintreten kénnen. 
Das hat de Lavaud seinerzeit bewogen, die Starrachse fiir die Vorderrader abzulehnen, die bei Per- 
sonenkraftwagen auch nicht mehr vorkommt. 

Der Vorteil der Starrachse besteht vor allem darin, daB beim Durchfahren einer Kurve die Rader 
die Seitenneigung des Aufbaus nur im Ausmaf der Reifeneindriickung mitmachen und da8 daher im 
Ablauf der Bewegung des gestérten Fahrzeugs die gefahrlichen Kreiselmomente der Vorderrader um 
die Lenkzapfen klein bleiben. AuBerdem liegt der Momentanpol des Wankens wesentlich héher als 
bei Parallelfiihrung. 

Abb. 10 zeigt den Einflu8 der stark iibertrieben gezeichneten Reifeneindriickung auf die Lage 
des Momentanpols des Wankens. Fiir starre Reifen wurde der Momentanpol M, in der Héhe der 
Achse angenommen (w =r), er kann aber noch hoher liegen, 
wenn die eingespannten Hauptfedern biegesteif sind. Fiir un- 
endlich harte Hauptfedern und elastische Luftreifen ist der 
Momentanpol M ,.. Fir elastische Hauptfedern und Luftreifen 
liegt er in M. 

Aus dem Sinussatz im Dreieck M, M,, M folgt fiir kleine 
Winkel 


— = “R | (15) 


Der Kippwinkel der Rader Ap verhalt sich zur Seitenneigung 


des Aufbaus J wie die Reziprokwerte der zugehérigen Federkon- 
Abb. 10. EinfluB der Reifeneindrtickung auf die Lage AR cA 3 4 2 B 
des Momentanpols des Wankens bei Starrachse. stanten: — = —, worn c, = und Q= sind. 


$s 
A ch Ra (F) Ch + ch 
Bei den iiblichen Verhaltnissen von cy, zu cg und f zu s wird Ap/A = 0,1 bis 0,2 gegeniiber Ap/A = 1 bei 
Parallelfiihrung. Damit sind bei Fahrzeugen mit starrer Vorderachse die Kreiselmomente der Vorder- 
rader um ihre Lenkzapfen beim Durchfahren einer Kurve erheblich kleiner als bei Parallelfiihrung. 
Dazu kommt, dafs wegen der hohen Lage des Momentanpols auch der Hebelsarm der Fliehkraft kleiner 
wird. Fiir Lastkraftwagen, bei denen man meist von vornherein mit hoher Schwerpunktslage und 
groBen Drehmassen der Rader rechnen muB, erscheint deshalb die Starrachse besonders geeignet, 
-wenn man nicht auf héhere Fahrgeschwindigkeit verzichten will. 


* 
Cr eR 


12. Uberpriifung der Lenkstabilitat eines Lastkraftwagens.1 Fiir das behandelte Fahrzeug wurde ein 
Gewicht von 6000 kp (bei 2/3 der héchsten Nutzlast) ein Radstand / = 3,1 m, eine Spurweite s = 1,6 m 
und eine Héchstgeschwindigkeit v,,,, = 80 km/h angegeben. 

Die Hohe des Schwerpunktes iiber der Fahrbahn wurde mit 1 m angenommen. Fiir die restlichen 
Konstanten wurden die nachstehenden Werte gewahlt: 

Aufbau: I, = 300 kpms?, k; = 0,975 kpm~s?, h, = 1,2 m, 

€ =0,0635 fir v= 80km/h, wv =1. 


1 Die Voraussetzung, daB der Schwerpunkt in der Mitte des Radstandes liegt, ist hier nur beim unbeladenen 
Fahrzeug halbwegs erfiillt. Mit zanehmender Beladung verlagert sich der Schwerpunkt nach hinten und oben, 


‘ 
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Hauptfedern und StoBdampfer: c, = 28000 kp/m, f = 1,2 m, 
ky = 935 kpm='s. 

Rader: r =0,46m, I,, = 0,4 kpms?, I, = 0,3 kpms?, n = 0,02 m, 
Ce = 75000 kp/m, k = 7500 kp/rad. 

Damit werden: C = m/4k = 0,02 m7s?, 

Cp = 133000 kp/m, c, = 23300 kp/m, Ag/A = ec, = 0,175, 

Aw = w(dp/A) = 0,08 m fir w=r, h=1m—(w—Aw) 0,6m, 

I, = 1,— mh? = 84 kpms*, B= 30000kpm, d = 1345 kpms. 


Mit diesen Werten werden wieder die Koeffizienten der charakteristischen Gleichung a, bis a, 
nach (6) errechnet, worin a, bis a; um die in (13) nicht unterstrichenen Glieder aus den Kreiselmomenten 
a3), bis a, vergroBert wurden. Aus dem Nullsetzen des Imaginiarteils der Funktion (8) 

6 


fG)y = Dd a,(vr)°~ 
n=0 
folgt mit (9) 
Vea ot ye =e 51,5 
und mit (10) fi) = + 37,0-104> 0, 
fir.) = — 388 -10?< 0. 

Der Zustand des Fahrzeugs bei Héchstgeschwindigkeit ist demnach stabil. Wenn fiir die Vorder- 
rader statt der Starrachse eine Parallelfiihrung gewahlt wird, wird bei im iibrigen unveranderten 
Annahmen A;/A = 1, h = 0,8 m und J, = — 84 kpms?. Die in diesem Fall um ein Vielfaches gréBeren 
Kreiselmomente der Rader und der gréBere Hebelsarm der Fliehkraft um die Wankachse bewirken 
Unstabilitat des Fahrzeugs bei derselben Héchstgeschwindigkeit. 


13. Die Vorspur der nicht angetriebenen Rader. Bei Hinterradantrieb sind die Radebenen der beiden 
Vorderrader nicht parallel, sondern schneiden sich in einer annadhernd lotrechten Geraden, die in 
der Fahrtrichtung weit vor dem Fahrzeug liegt. Diese Eigenschaft nennt man Vorspur der Vorder- 
rader. Die stabilisierende Wirkung der Vorspur hangt nun so eng mit dem Wanken des Fahrzeugs 
zusammen, dafs sie hier kurz gezeigt werden soll. Wenn man in Abb. 1] die Reifeneindriickung ver- 
nachlassigt und die Neigung der Querlenker J, der Seitenneigung / fiir A 
kleine 1 proportional setzt, so bewirken die beim Wanken auftreten- 
den Krafte in Hauptfeder und StoBdampfer eine Radlastanderung 


AN = ep lpg (Inlln) + Bale dg (lelPn) = (5/2) (Cp +h 4), (16) 
worin Ap =(s/2 Ip), ¢; = Cp (Ip/ly)? und ki =k, (Ip/lp)? ist. Der 
EKinflu8 der Kreiselmomente der Rader auf die Radlastverteilung 
ist in (16) nicht beriicksichtigt. Die Wirkung der Vorspur besteht 
nun darin, daB das kurvenduBere Rad mit der gréBeren Radlast 


auch den gréBeren Schraglaufwinkel hat. Die Summe der Seiten- 
' krafte der Vorderrader wird 


G s * *4 
Sint Suhel + 3 GA + hid) 


G * *4 
+ laa +52) 


X14 


a 


und mit Ore = Oy -+ Xoo On; es Oy — X% (Abb. 11), sowie k = ky(G/A) 
kommt 


ES, = 2 boy + hy s (CA + 5A) a - (17) 


Durch die Vorspur wird die Summe der Seitenkrafte der Vorder- oe by ee oben) tne Sette elaat. 
rader um das zweite Glied rechts in (17) vergréBert. Dieser Aus-  Vorderrider beim Durchfahren einer Kurve. 
druck mu8 in das System der Bewegungsgleichungen (4) eingefithrt 
werden, wodurch alle iibrigen Freiheitsgrade mit dem Wanken des Aufbaus gekoppelt werden. Die 
Ausrechnung der um die Vorspurglieder erweiterten Determinante (5) ergibt Erganzungsglieder zu 
den Koeffizienten der charakteristischen Gleichung nach (6), vor allen zu a, und a,. Der weitere 
Rechnungsgang mit den Koeffizienten nach (6) und (13) zeigt, da die Vorspur die Stabilitat ver- 
groBert. Sehr wirksam ist dieses Mittel allerdings nicht, da der Winkel a, unter einem halben Grad 
bleiben mu, wenn nicht die Reifenabniitzung stark zunehmen soll. Auch in anderer Beziehung 
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ist die Vorspur ein recht bedenkliches Stabilisierungsmittel. Ihre Wirkung beruht namlich auf der 
Verschiedenheit der Radlasten der kurvenduBeren und kurveninneren Rader, die man zur Erzielung 
einer méglichst groBen Gesamtseitenkraft in der Kurve gerne klein halt.’ Die Vorspur wird wohl 
auch deshalb so haufig angewendet, weil sie vielfach am betriebsbereiten Fahrzeug in der Werk- 
statte leicht verandert werden kann. ; , 

Wenn die Rader nicht parallel gefiihrt werden, kann man sie so lagern, dab die Vorspur mit der 
Einfederung wachst und dadurch ihre stabilisierende Wirkung stark erhéhen. Die genaue Kinstellung 
der Rader ist in diesem Fall besonders schwierig, weil von der Belastung abhangig. AuBSerdem ent- 
steht beim Durchfahren enger Kurven wegen der zu groBen Vorspur des auferen Vorderrades ein 
beachtlicher Kurvenwiderstand. 

14, Zusammenfassung. In der vorliegenden Arbeit wurde der Kraftwagen als lineares System von 
vier Freiheitsgraden — seitliches Schieben, Gieren um die Hochachse, Lenkschwingung der Vorder- 
rider und seitliches Wanken — behandelt. Dabei wurde zunachst der Schraglauf der Luftreifen, 
der Nachlauf der Vorderrader, die Radlastanderung durch den Luftwiderstand, die Elastizitat der 
Hauptfedern und Luftreifen und die Dampfung des Wankens durch StoSdampfer beriicksichtigt. 
Danach wurden auch die durch das Wanken, das Gieren und die Lenkschwingung geweckten Kreisel- 
momente der rollenden Rader bei verschiedenen Radaufhangungen in die Rechnung einbezogen. 
Es zeigt sich, daB diese Kreiselmomente bei Windstille die Hauptursache einer bei hohen Geschwin- 
digkeiten auftretenden Unstabilitit des Kraftwagens sind. SchlieBlich wurde noch der EinfluB der 
Vorspur besprochen. 

Die Uberpriifung der Lenkstabilitat von drei verschiedenen Kraftwagen an Hand des oben be- 
schriebenen Ersatzsystems zeigt, daB alle drei bei der jeweils erreichbaren Hichstgeschwindigkeit 
zwar stabil sind, aber sich in der Nahe der Stabilitatsgrenze befinden. DaB dies fiir die meisten Kraft- 
wagen zutrifft, 1aBt sich aus der weiten Verbreitung der Vorspur schlieBen, wobei oft die Stabilitat 
erst durch die richtige Einstellung der Vorspur zustande kommt. Ohne zwingenden Grund wiirde 
man sonst kaum zu diesem bedenklichen Mittel greifen. 

Das lineare Ersatzsystem zeigt das Stabilitatsverhalten eines Kraftwagens und 14Bt den EinfluB 
der physikalischen GréBen der einzelnen Bauteile auf die Stabilitat des Ganzen deutlich erkennen. 
Der Fehler, der durch die Vernachlassigung bestehender Nichtlinearitaten, vor allem in der Beziehung 
zwischen Schraglaufwinkel, Radlast und Seitenkraft entsteht, diirfte bei kleinen Schwingungen nicht 
allzu gro sein. Bei der praktischen Auswertung der Ergebnisse ist zu bedenken, da man stets von 
einem gegebenen Aufbau samt Motor und Innenausstattung ausgehen muB. Fir den Aufbau eines 
schnellen Fahrzeugs muB eine niedrige Schwerpunktslage, eine kleine Drehmasse um die Langsachse, 
eine groBe Drehmasse um die Hochachse, eine windschnittige Form und ein Radstand gefordert 
werden, der gro8 genug ist, um eine zu starke Entlastung der Vorderrader durch den Luftwiderstand 
zu verhindern. Das alles sind bereits bekannte Grundsatze fiir den Konstrukteur, die in der Arbeit 
nur bestatigt wurden. 

Wenn die obigen GréBen fiir den ausgefiihrten Aufbau bekannt sind, kann dazu eine Federung 
entworfen werden, die einen hohen Riickstellbeiwert des Wankens, eine starke Dampfung und eine 
nicht zu tiefe Lage der Wankachse haben soll. AuBerdem soll die Drehmasse eines Vorderrades um 
seinen Achsschenkel und damit der Raddurchmesser so klein sein, wie es die Reifenbeanspruchung 
bei hohen Geschwindigkeiten zulaBt. SchlieBlich sollen die Vorderrader so aufgehangt sein, daB sie 
schon bei kleinen Schwingungen die Seitenneigung des Aufbaus nicht voll mitmachen. 

Zuletzt wahlt man noch die KenngréBen der Lenkung, d.s. die Drehmassen eines Vorderrades 
um seinen Lenkzapfen, den Nachlauf der Vorderrader und die Seitensteifigkeit eines Luftreifens 
und kann jetzt den Zustand des ganzen Systems berechnen. Ergibt sich Instabilitat, so wird man 
zunachst die Gréfien der Lenkung — innerhalb technisch méglicher Grenzen — in einem die Stabilitat 
fordernden Sinn zu andern suchen. Wenn das nicht ausreicht, kénnen auch noch die GréRen der 
Federung — z. B. durch Annahme einer Zusatzfederung — verbessert werden. Die Kenngréfen 
des Aufbaus, die nur am ausgefiihrten Stiick verlaBlich gemessen werden kénnen, diirfen dagegen 
nicht variiert werden. 

In der vorangegangenen und in dieser Arbeit wurde versucht, den modernen Verfahren zur Sta- 
bilitatsbestimmung das Gebiet des Automobilbaues zu erschlieBen, auf dem von Beginn an haupt- 
saichlich empirisch vorgegangen wurde. 


(Eingegangen am 2. Juni 1959.) 
Anschrift des Verfassers: Dr. Gert Mitterlehner, Wien II, Halmgasse 1 (Osterreich). 


1 H. Scherenberg, Automobiltechn. Ztschr. 57 (1955) S. 155. 
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Kine Verallgemeinerung des Prinzips vom Minimum 
der potentiellen Energie zweidimensionaler elastischer Kontinua 


Von D. Riidiger 


1. Einleitung. In der vorliegenden Arbeit werden die Gedanken, die zu einer Verallgemeinerung 
des Prinzips vom Minimum der potentiellen Energie elastischer Kérper fiihren!, auf zweidimen- 
sionale elastische Kontinua iibertragen. Die durch Einfiihrung der Randersatzkraft entstehenden 
Schwierigkeiten kann man durch Aufteilen der Randkurve C in differentiierbare Kurvenstiicke 
beheben. Da Randbedingungen des Variationsproblems nicht auftreten, sind fiir beliebige Rand- 
wertaufgaben der Schalen-, Platten- und Scheibentheorie die Verfahren von W. Ritz und E. Trefftz 
hinsichtlich der zu verwendenden Funktionen ohne Einschrankungen anwendbar. 

Zur Ableitung des Variationsprinzips werden die Grundgleichungen der technischen Schalen- 
theorie? in allgemeinen Koordinaten 


NET, Pp =0, 
Q*|, + bag Ne? + p? =0, (1,) 
My O° =0, 
Nae = =" (De? 4 ye%7 2 9D,,), 
mee = — Eh OM (yoo 4 y on 669 W),,) Cs) 
" lv 12 78 


verwendet. In (1,) unf (1,) stellen N*? den Langskrafttensor, M*’ den Momententensor, Q* den 
Querkraftvektor, p*, p? den Vektor der Flachenbelastungen, 


Ds6 = > (V*? + Vels) — bee (14) 
den Deformationstensor der Schalenmittelflache, V*, W den Verschiebungsvektor und b,, die 


Koeffizienten der zweiten Fundamentalform der Flachentheorie dar. Sind diese Null, so zerfallen 
die Schalengrundgleichungen (1,) und (1,) in die Grundgleichungen der Plattentheorie 


O|,+p?=0, 

Me*|,— QF =0, (2,) 
Eh wh a 

Lee pear er (W\a8 +y7e 7 666 Ws) 


und in die der Scheibentheorie 
N wes + pP =0, 


1 Oo 
eae | YD (24) 
Nob — ae (D«68 + ve? gh Ds) : ; 


Legt man den Tangentenvektor der Randkurve C durch e*, ihren Normalenvektor durch n, =, 8 is 
fest, dann ist die am Schalenrand angreifende Langskraft mit der Randnormalkraft N = n, ng N* 
und der Randschubkraft S = e, ng N*? durch 


Ké — n, N*? = Nné + Se’, (3;) 
der am Schalenrand angreifende Momentenvektor mit dem Randbiegemoment B =n, ng M8 
und dem Randtorsionsmoment T = — n, eg M*? durch 

Me = fn” My, = Be’ + Tn? (32) 


egeben *. Fir die Randquerkraft gilt 
= Q=n, Q. (33) 
1 D. Riidiger, Ing.-Arch. 27 (1960), S. 421. 
2 Vel. ae D. Ridiger, Festschrift Richard Grammel. Ing.-Arch. 28 (1959), S. 281. 
3 H, Neuber, Z. angew. Math. Mech. 29 (1949), S. 144. 
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2. Das Variationsprinzip. Eine elastische Schale mit der Flaiche F unterliege den Flachen- 
belastungen p*, p® und sei beliebigen Randbedingungen unterworfen. Dabei sollen an den Rand- 


stiicken a z und 5 die Randnormalkraft N, die Randschubkraft S, die Randersatzkraft Q* ae ioe 
2 — — 
—@T/ds und das Randbiegemoment B gegeben sein, wahrend an den Randstiicken c, ¢, ¢ und c 
die Randverschiebungen n* V,, e* V,, W und die Randverdrehung n* W|,, vorgeschrieben sind. 
1 a 2 a 3 3 Ve ee 
Es ist dann C=c+¢, C=ce+c¢,C =e+e und C=c+e. 
Die Ableitung des Variationsprinzips erfolgt mit der Formanderungsarbeit — 6l7 der im 


Gleichgewicht befindlichen Schale bei dem Verschiebungszustand dV*, dW. Diese ist mit dem 
Flachenelement df der Schalenmittelflache : 


et odd MN a BS) OL g hese: Lf (Qa =f Dy NUP SP pt) OW df ¢ 
+ [fA Ml, + OSH, f= 0. (t) 
Fiihrt man in (4,) die Identitaten 
Ne?|, dV, = (N*? 6V 5) |g — N*? OV ela » 
Q*|,0W = (Q OW) |. 08 OW, ’ 
M*6|,OW|5 = (M*? 6W|5)|,— M8 dW\g, 


ein, dann ergibt sich bei Beachtung von 
1 


dD xp => (OValp + OV |x) — bap DW 


die potentielle Energie 6/7 zu 
OlT = ff N*? 6D,,df — ff M*8 dW|,, df — 
F F 


— Jf (N° V p)]a af + Sf (Me? 9W|p)|5 df — Lf (Q° 0) |. af — 
— Jf (p* 0Vq + pt OW) af = 0. (4,) 


Fiir die beiden ersten Flachenintegrale in (4,) kann mit dem Elastizitatsgesetz (1,) 


H Ne" Dy df— | | M*? OW |pa df 
F F 


= [| Eh (D*? 5Dyp + €*? 6° D, 4 5D,,) of — 


1] — 7? 


i 
| Bly WA hag ae? 
A = gage (VP OW |, cee lg OW Ie) of 

F 
ripe | Eh 5 (D*8 x? fd 
=e fp ee OO Dag tv Ve Dig Dys)| af — 
F 
1 Bh ht 
—$[ [PA BOPP I Wig treo Wl. Wp af 
J 
I | 


| O(NF Dao) af — | [ A(M** Wap) af 


geschrieben werden. Beachtet man noch die wegen (1,) giiltigen Beziehungen 
N*8 Dyg =(N*? V,)|g + p* Vz — byp NP W, 
— MP? Wap = — (M*? Wadia + (Q* Wis + bag N86 W + pW, 
so folgt schlieBlich 


[ [27 era—[ revamped} | four raiar— 


— yf foarewine art > | foe Wyle af +3 | | Oe" V+ p? W) df. 


Ria Sexe 
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Damit erhalt man fiir die potentielle Energie 


sm = 5 | | oe V,)l — a {| Mer Widipar +5 [fp (O° W))y df — 
=i (NF OV pla af + | | (Me? §W7|,) a | f (Q° OW)|, df + 


1 3 3 Aeon o 3 SS 
+z [re V, + pW) df [Je dV, + p> OW) df =0. (45) 


Die Anwendung des Gaufschen Integralsatzes 


Lf Ole =$ 0) meds 


auf die ersten sechs Integrale in (4,) ergibt 
1 1 1 
CUTS r5 pans V.) ng ds ee W\,) mg ds + Pre W) n, ds — 
—fovvor, n, ds -paeronin ufo ds + 


+4 >| foe V+ p> W) af — [Jo oVq + pW) df =0, 
und nach Kinfiihrung der Randkrafte und Randmomente (3) wird 
: ieee epunny EY dsc rps yds 


+1 goo W) ard ga (T e* W\,) ) ds — 3h (Bm Wi, \ds 
Cc Cc ¢ 
<P Nw Va ds —G S & bY, ds — 
3 . 
— P90 ds —P To OW), ds + GB m dW, ds -- 
Cc G Cc 
oo | | 5 (p* Va + p? W) df — | | (p" Va + p?OW) df =0. (44) 
F F 


1234 =i2it 
Ersetzt man noch die Randintegrale in (4,) durch die Linienintegrale langs ¢, c, c, c und ¢, ¢, ¢, ¢, 


so folgt 
iy ee ey 4) da eg. x | (Vm Va) es Hee! jase y | 80S @ Fa) ds + 


8 
° Rees 


ce 


+z [Ow + Te wi) ds +z [QW + Te W),) ds— 


c 


° rer 


—; | n® Wa) yds fae n* W],) ds — 
A 


° nage 


— f now V_) ds— | Now on EAR Als a) ds— [ SOY, ane 


i 


© ||to 


c c 


zal [QoW + Td(e* W\,)] as— [ [QoW + Tle W|,)] ds + 


9 lle 
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2 
[ Te |. ds = TW] { 
y! 


118 . Riidige r 
se | B6(n* W|,) ds + | BO6(n* W\,) ds + 

= zs 

4 [Joe V+ p?W) v—| Jo 5V, + p?6W) df =0 
Die Linienintegrale, die das Randtorsionsmoment enthalten, werden durch partielle Integration 
Man erhalt auf einem Randstiick zwischen zwei Punkten J und 2 
OG i 
an W ds 3 
i 


umgeformt. 
i 


3! 


~ 


ld ee 


3 
Randt Rand © mom77£77, 


J 
Abb. 1. Aufteilung des Randes 
so daB zum Beispiel fiir eine Schale bzw. Platte mit den vier vorgegebenen Randstiicken entspre- 


5 ines 


ces 


° ee 


ea ; 
ae Pag dss | Vm Y, ) ds +> gt V,) ds + 7 [Se V,)ds + 


ce 


+ [Or W) ds +5 IT WH +3 IT WIE 


1 
z 
1 
f [oe W) ds +3 o(T WR +4 d(T W] 


ae V,)ds— | Soe Va) \ ives 


ae n° W,) )ds— 5 | (Bm Wh) ds 


° |lee 


@ |le 
a ieee 


— [ o V,) ee 


a Q* OW ds — (TowR— [Tow}s — 
(rw + 
(4) 4 


—f orowa— trom 
ag [ Bo ne Wa) ds+ | Bow Wa) ds + 
[ (ov. + pom) af =o 


9 lle 
© |= 


+4 [ [oo ¥. +p Map 
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geschrieben werden kann. Fiihrt man in (4,) die Identitaten 


— Nd(n*V,) =—d(Nn* V,) + dN nv V,, 
—S6(eV,) =—6(Se*V,) + dS e* V,, 
—Q*6W =—6(Q*W)+6Q*W 

— Tow =—d6(TW) + oT W, 


B0(n* W|,) = 0(B n* W\,) — OB n* W\, 
ein, dann wird zunachst 


oI = + | (Nm V,) eee “ag Vs watt 5 ae 


© 


+ O(T WR + oT WI —y | Be W|,) ds + 
+z [ (Ne V,) Gee sel (SeV a)ds+ 3 [ 0(Q* W) ds + 
i EB 
+5 O[T WY +5 O(T WE <y | Bm Wy) de — 


ude 
¢ 


— [ Noor V,) ds— | SO(e* V4) ds — [.Q* 9m as — 
1 2 BSS 


c 


—d[T WR —o[T Ww} a | B 8(n* W|,) ds + 


o |= 


+ [6T W]? + [6T W}i— 


(Nm V,) do — | 0(S e* Va) ds — { ar) Vee 


i 2 i 
+ [ON mw Vids + [ dS 0 Vas + [ 00" Was — 
Fé i 3 
— d5[T W]¥ —6[T WY + [ (Bm W\,) ds + 
+ OT WE + OT WIE — | OB mW ds + 


3 ||» 


+z [for Va+p? W) af — | | (pV, + pow) df= 0 
“P F 
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und nach Zusammenfassung 
; : if 
OT = 5 [ AN w Va) do + > | OS e Va) ds +3 | (Q* 7) ds — 


nS ee 
c c c 


1 1 
— 7 AT WE 7 AT WIs— a | Bm W),) ds — 
4 


=| ae ray [90¥ V,) ds — x | rH) as— 


1 , 1 / Tee off 
—4 oT WIE —4 OIF WE +> "1 0(B n° Wq) ds — 


© || 


1 


— [ Nom V, s as— | $0( S d(e | Q* 6Wds + 


+ [oT We + een | B 0(n* W\,) ds + 


a [ oN ms Vads + [ OSeVeds + | dQ*Wds + 


3 


a 
c 


9 ||r 


+ [OT WE + [OT WY — | 8 n° W), ds + 


+z [ [ow Vat pm) a— | | (pr V+ pow) df =0. (5) 


Die Beziehung (5) ist die modifizierte Energichauptformel, aus der sich sofort das Prinzip vom 
Minimum der potentiellen Energie zweidimensionaler elastischer Kontinua ableiten 1aBt. 


tamer oat) eee ee a 

namlich in den Randstiicken ¢, c, ¢ und = die Randkrafte N, S, Q* und das Randmoment B, in den 
ie ee oe == = 

Randstiicken c, ¢, ¢ und ¢ die Randverschiebungen’n* V,, e* V,, W und die Randverdrehung n* W,, 


in den Randpunkten 1, 2, 3, 4 die Verschiebung W und auf der gesamten Oberflache die Flachen- 
krafte p*, p® fest vorgegeben sind, gilt 


: | 
emer pe ee ae W ds —— [TW g (TW 3 [ BrW|,ds— 
3 


c Cc €c 


=}ver, ds— 7 [SeV, ds— 5 [ Q* Wds— (TWH TW set 


= aa a3 
Sein ae ioe ra+ viel els [Be W),ds + 
fret ds + [SeV, eecncre ATW [Bn n* WI, ds + 
zs ae, A 


Hffenstemas |e V,+ PW) df|=0. 
Fr F 


Da 


Ser | 


OEE eerste Or Ry yh Oe 
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Hieraus ergibt sich, daB unter allen méglichen Verschiebungen eines zweidimensionalen elastischen 
Kontinuums, die in allen Randpunkten beliebige Werte annehmen kénnen, diejenigen 
eintreten, die den Ausdruck > 


= e & 
Het (N—2N) w Vds + | (S—23) Vads +3 | (Qr—2 Q*) W ds — 


1 


c c c 


=a (oe? B) n* W\,ds —5-[T(W—2W)R}—2 [r(w—2 WE — 


1 = — 2S 
<2 [Nw (2 V,) ds—z [Se(v,—2 V4) ds—4 [or(r—2 W) ds + 
1 == 1 = ay ay 
+z [Be (Wl,—2 W\.)—[T (W—2 W)% — [7 (w—2 Wy + 
1 fe nt E. 
+a [ | le —2 7) Ve + 27) Wa (6) 
F 


zu einem Extremum werden lassen. Das Extremalprinzip (6) hat gegeniiber der klassischen Formu- 
lierung den Vorteil, daB an die infinitesimal benachbarten Systeme von Verschiebungen hinsicht- 
lich der Randbedingungen keine Forderungen gestellt zu werden brauchen. Man kann also zur 
numerischen Darstellung der Lésung V* = V*%(x!, x?), W = W (x1, x”) eines vorgegebenen Rand- 
wertproblems der Schalen-, Platten- und Scheibentheorie im Sinne von W. Ritz einen Ansatz 

n nN ya» n Tea 3) 

Ves > 9", LO aA . (7,) 

= 


ol 
v v ° . . ° e 
mit willkiirlichen Funktionen g*, «@ machen. Die Funktionen sind dabei nur so allgemein auszu- 


wahlen, daB bei hinreichend groBem n und passend gewahlten Koeffizienten Oh ¢ die Lésung der 
vorgegebenen Randwertaufgabe mit beliebiger Genauigkeit dargestellt werden kann. 
Das System von Verschiebungen (7,) liefert fiir den Deformationstensor, den Langskrafttensor 


und den Momententensor 


n n 1 V,0 v v,B v y 
Det 5 Ea pee pr) — HP ea] ; 


vy=1 
n Eh nm n 
1 te (Dep + vet” oD), | (72) 
n Eh he ( . ; ) 
ME org (ce v POW], 5 


und fir die Randkrafte, Randmomente und Flachenbelastungen 
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Die Koeffizienten c, C (VES, 2,06, hy pos Ansatzes (7,) werden nun so bestimmt, daB die 
mit (7) gebildete potentielle Energie 


n _ 


fr=1 {(w¥—25) n* Vous +4 | (S—23) e* V,, ds + > (o*—2 0°) Wise 


2 
A 2 RS 
ed (B—2 B) me Was — 3 |t (2 w)\,—+|r(w—2 w)|, — 
aay 
—_ [Ne Gacy V,)ds—1{ Se(v,—27,) que L[ do(—2 W) ds + 


+3 | Be(,—2 7) a—+[r(7_2 Ww), 4 [tlw —2 W)|y + 


+2 | | (2m) 0. +02) (8 


dem Extremum méglichst nahe kommt. Die Bedingungsgleichungen fiir das Eintreten des Extrem- 
wertes lauten 


SE ee ey epee ND th ee (9) 
dc dc 


und haben die Galerkinsche Form, 
3. Das Verfahren von Trefftz. Beim Trefftzschen Verfahren wird der Ansatz 
= Deg, W= Seo (10) 


verwendet, in dem die Funktionen ", ow Partikularlésungen der homogenen Schalengrundgleichun- 
gen (1) sind und entsprechend dem Variationsprinzip (6) in den Randpunkten beliebige Werte 
annehmen kénnen. Bildet man die dem Verschiebungszustand (10) entsprechenden Randkrafte 


nn n n n n 


N, S, Q, Q* und Randmomente B, T, so ergeben sich mit 
a: 1 » le va ed 
d aw +985) — bP wo, 


nxe = i (axe a py 6x7 gBd ‘ia 


y Eh cae | eee 
i oi epee er ale? + ven? 69 ol), 
v v - v v y 
go = m8, ne], = 0, 1a! + beg n® ? 8 
und 
v ps » y 
n=n, ng n?, Sn een , 
Y y v 5 
bean igi t=—n,e,m*", 
v at ) 
v v v t 
= 71 Ta no ee ee Me’ 
q a ? q q Os 


eas 


PPE anew AN eye be oo, 


/ 
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die Ausdriicke 


v= 5 
n Loe ay (11,) 
Cs = c q* 4 

y=1 
n Ce tpl! 
Bet ee: bs 

y=1 

Ny 

ie ct 

y=1 


Da die Funktionen (10) voraussetzungsgem4B Partikularlésungen der homogenen Schalengrund- 


gleichungen sind, wird 


n 


om 0) (11,) 


Die n linearen Gleichungen fiir die Konstanten lauten mit (11) 


cea bem a tees 


1 


y= 0=5 [ (nw V, + N nt y—2 nt oy) ds + 


wy 
c 


1 Pas) n n bee * 
th | Geb, +S eG, —25eGa)ds + 
oo 
1 an fia, n io ae Arh 
Pair ealigeas 
3 


3 | 


wet | (bn Wl, + Bre oh —2 Bowe 


TE 
ae 
yon ny y 2 von By ’ 4 
Sw + To—21W,—+ hw + To —21 P|, — 
sodef (ive Vi + Nn, Q,—2nn Vy ds — 


ei parentet= v7 7)a— 
2 
Lape [Gr 4 Or 24°) ds + 


pi 


+4 [lm Wi, + Bre o|,—2bn° W),) ds — 


4a 
Pe Ny (De RA yn ny | fe 
—| [leo +R e)a=0 Gaon eal (12) 
F 


g* 


f 
is 
| 
if 
| 
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Bei Anwendung des Bettischen Theorems! fiir eine Schale mit Randbedingungen entsprechend P 
Abb. 1 gilt 


Pe joins [Seine] pia— | Bednaslah 4 [Tol + 


ce 


Le y POS) Ue a a y Nn y 13/ nm vii’ 
+ iiiuee (gee, a4 [Ga bwol.ar+ [tab +[tol 
an 
lees V, ds + 
= 


tne peas |e Vis | q* Was | bw Was +(e Wy +|e wl. 


= 
n yp n]2 yn ]4 
set Vs J qe Was—| bm Whe as + [ew + lew) + 


c 


° Jooe ° |looe. 


+ 


0 [le 


© | {ro 
@ {|e 
© |= 


y : Hh 
und damit kénnen die n linearen Gleichungen fiir die Konstanten ¢ (vy = 1, 2,..., n) in der Form | 


[(w— N) hed + | (8—s) o* Py ds +[(e—@) 5 a— 


c ¢€ 


— | (61) w dull wh _-[@ wk 


V 


iw (7), —,) as —[e(w— w)],—|t (w— )], — 


—| | (e+ pa) a=o isd eee | (13) 


angegeben werden. Lat man in (10) die Anzahl der zur Approximation verwendeten Funktionen |} 


tiber alle Grenzen wachsen, dann liefert (13) die Koeffizienten des gegen die exakte Lésung konver- 
gierenden Ansatzes. 


4, Zusammenfassung. Es wird ein Variationsprinzip zweidimensionaler elastischer Kontinua 
aufgestellt, das als Sonderfalle das Prinzip vom Minimum der potentiellen Energie und das Castig- 
lianosche Prinzip enthalt. Da Randbedingungen des Variationsproblems nicht auftreten, kénnen }} 


die bei den Verfahren von Ritz und Trefftz verwendeten Funktionen in allen Randpunkten be- | 
liebige Werte annehmen. 


(Eingegangen am 12. Juni 1959.) 
Ansehrift des Verfassers: Professor Dr.-Iug. D. Riidiger, Freiberg (Sachsen), StraBe der Einheit 12. 


+ Vel. z.B. D. Riidiger, Z. angew. Math. Mech. 40 (1960) S. 114. 


Wibbiaiti eas ncn 3 0005 


XXIX. Band 1960 H_. Schaefer: Analogie zwischen den Verschiebungen und den Spannungsfunktionen 125 


Die Analogie zwischen den Verschiebungen und den Spannungsfunktionen 
in der Biegetheorie der Kreiszylinderschale 


Von Hermann Schaefer 


1. Einfiihrung. In friiheren Arbeiten’? habe ich darauf hingewiesen, da® in der klassischen 
Elastizitatstheorie des zweidimensionalen ebenen Kontinuums (Scheibe-Platte) eine vollstandige 
Analogie besteht zwischen den 3 Verschiebungen 1, v2, v; und den 3 Spannungsfunktionen ®,, ©,, D, 
in der angegebenen Reihenfolge. Dabei sind v, und v, die Verschiebungen in der Scheibenebene, ®, und 
®, die Spannungsfunktionen der Plattenbiegung. Altbekannt ist die Analogie von Timpe-Wieghardt 
zwischen v, der Durchbiegung der Platte, und ®,, der Airyschen Spannungsfunktion der Scheibe. 

Beim gekriimmten zweidimensionalen Kontinuum, der Schale, ist eine solche Analogie bislang 
noch nicht bemerkt worden. Nun stehen bei einer sinnvollen Ubertragung der Elastizitatstheorie 
von Scheibe und Platte auf die gekriimmte Schale mancherlei Méglichkeiten offen. An Biegetheorien 
diinner Schalen herrscht deshalb kein Mangel. Den meisten unter ihnen sieht man jedoch auf den 
ersten Blick an, da in ihrem Rahmen unsere Analogie sicherlich nicht gilt. 

In dem Wust der Schalengleichungen scheint mir unsere Analogie wertvoll als Ordnungsprinzip. 
Natiirlich muB gezeigt werden, daB dieses Ordnungsprinzip nicht den Rahmen der oben erwahnten 
Méglichkeiten sprengt. DaB die Analogie zwischen Verschiebungen und Spannungsfunktionen in 
sinnvoller Weise auf die gekriimmte Schale tibertragbar ist, soll am Beispiele der Kreiszylinder- 
schale dargelegt werden. 


2. Die Gleichgewichtsbedingungen. Das zweidimensionale Kontinuum der Kreiszylinderschale 
_ habe den Radius a. Die Koordinate in Richtung der Erzeugenden sei x,, in Umfangsrichtung %2. 
Die Schale sei nur an ihren Randern belastet. 


Abb. 1. Schnittkrafte und Verschiebungen Abb. 2. Schnittmomente und Drehungen 


Nach Abb. 1 lauten die Gleichgewichtsbedingungen fiir die Schnittkrafte 


O,Ny1 + 02Nox =0, ‘ (2,1) 

1 
0,Ny2 + 02No2— Fi Q, =0, (2,2) 

1 
0,02 a 050; ee No» =0, (2,3) 

und fiir die Schnittmomente nach Abb. 2 
0,Mz; — 0,.M,;—Q =0, (2,4) 
1 é 
0,M22 a 0.My» a Q. aca R, =0, (2,5) 
1 

— i (0,R, aig Rs) =F: (Nie as Np) =F et My» = 0 (2,6) 


1 H, Schaefer, Osterr. Ing.-Arch. 10 (1956), S. 267. 
au. a Abh, Braunschweig. Wiss. Ges. 5 (1956), S. 142. 
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(es bedeuten 0; = 0/0x;; 1 = 1,2). R, und Rg, die Momente in Richtung der Schalennormalen, 
werden nur voriibergehend eingefiihrt und spater Null gesetzt. : 


3. Die Spannungsfunktionen. Die 6 Gleichgewichtsbedingungen fiir die 12 SchnittgréBen erfillen 
wir identisch in den 6 Spannungsfunktionen Q,, 2,, Q, und ®,, B,, Ps: 


Niq= 00835 We leon (3,1) 
Ny = 0.2 + = Qz, No 240s (3,2) 
6:0, 0708 Or208 ) (3,3) 
AL ee pat, M,, = 2,9, , (3,4) 
Myy = 0,0, + ~ ®;, i pee (3,5) 
R, = 2,—2,0,+—®,, Rp = 2, + 4%. (3,6) 


4, Die Randwerte der Spannungsfunktionen. Aus der Darstellung (3,1) bis (3,6) der Krafte 
und Momente durch Spannungsfunktionen folgt 


— N,, dx, + Ny, dx, = dQ, , (4,1) 

Nye day + Nyy dy = dQ, + — O, dry, (4,2) 

Q, diy + Qy day = dO, — — Qy diy, (4,3) 

M,, dx, + (My, — Q;) dx, = d®, , (4,4) 
(My, + Q,) dx, + My, dxz = d®, + — ©, de, (4,5) 
(Ry — ,) dx, — (Ry — Q,) dx, = d®, — * ®, dry. (4,6) 


Wir betrachten einen einfach zusammenhangenden Bereich B auf der Schale, berandet von der 
Kurve C. Die Bogenlange s von C sei so orientiert, da bei wachsendem s der Bereich B auf der 
linken Seite von C liegt. C ist in der Parameterdarstellung 


‘ X, = %,(s), Xq = X9(5) (4,7) 
gegeben. 


Abb. 3, Randkrifte Abb. 4. Randmomente 


Die am rechten Rand von C angreifende Randbelastung halt sich am Bereiche B das Gleich- 


gewicht. Die Zusammenhange zwischen Randbelastung und SchnittgréBen entnimmt man den 
Abbildungen 3 und 4. 
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Demnach gelten auf C die fol 


genden 6 Differentialbezichungen zwischen Spannungsfunktionen 
und Randbelastungen: 


= Sey hee (4,8) 
oo Oe ke (4,9) 
=... “a0 
wi a 220, es (4,11) 
See Ot Sm, + 0; (4,12) 
a eee ay 


Sind k,, k,, k, die Randkrafte je Langeneinheit und m,, my, mz die Randmomente je Langeneinheit 
als Funktionen von s gegeben, so lassen sich die Randwerte der 6 Spannungsfunktionen 8 Fra. 0 03.8 Be 
®,, ®,, D, nach (4,8) bis (4,13) durch Integrationen langs C berechnen. 

Ohne uns hier mit der statischen Deutung der Randwerte unserer Spannungsfunktionen allzu- 
lange aufzuhalten, wollen wir nur vermerken, da der Vektor Q = (2Q,, Q,,.Q;) die Kraftesumme und 
der Vektor = (,, ®,, D,) die Momentensumme der Randbelastung eines endlichen Stiickes E 
von C ist, ® bezogen auf den Endpunkt von E. 


5. Die Deformationen und ihre Vertriglichkeitsbedingungen. Um die infinitesimalen Deforma: 
tionsgréBen unseres zweidimensionalen Kontinuums zu definieren, bedienen wir uns des Prinzips 
der ,,virtuellen Krafte“, das wir hier in der folgenden Form postulieren: 

Das iiber den Bereich B erstreckte Flachenintegral 


[ (2 Nnea + FO8. + Y Man + ¥ Rx) dF @=1,2; k=1,2) (5,1) 


soll fiir simtliche Kigenspannungszustinde N;;,, Q;, M;,, R; verschwinden. 

Diese Nix usw. sind die SchnittgréBen im Bereiche B bei unbelastetem Bereichsrande C. Sie 
miissen den Gleichgewichtsbedingungen (2,1) bis (2,6) geniigen. 

Natiirlich ersetzt man in (5,1) die ,,virtuellen Krafte“ N, , usw. durch die ,,virtuellen Spannungs- 
funktionen“ Q, We OR und @,, ®,, ®,, wie in (3,1) bis (3,6) angegeben. 

Dammit sind die Gleichgewichtsbedingungen abgegolten. Wegen des unbelasteten Randes C 
kénnen nach (4,8) bis (4,13) samtlichen 6 virtuellen Spannungsfunktionen die Randwerte Null 
erteilt werden. (Beilaufig bemerken wir, daB an inneren Lochrandern mehrfach zusammenhangen- 
der Bereiche die oo® Lésungen der homogenen Gleichungen (4,8) bis (4,13) eine Rolle spielen.) 


Nach Einfihrung der virtuellen Spannungsfunktionen in (5,1) wird partielle Integration an- 


gewendet (die Randintegrale verschwinden), nach 2,, 25, 23, ®,, ®,, ®;, geordnet und von der 
Unabhiangigkeit dieser 6 Funktionen Gebrauch gemacht. 

Die den SchnittgréBen zugeordneten Deformationen sind dann durch die folgenden 6 Gleichungen 
definiert : 


O4€1 — Oo811 + %1 ==/0, (5,2) 
Be Oia + ve —— ba sip (5,3) 
— (0xBr + G62) — G12 — an) += Ha = 9 (5,4) 
: SES ee oie Te (5,5) 
a oa Se (5,6) 


1 
— 0%, + 0% + paar is 0. (5,7) 


Ingenieur-Archiv 
128 


H. Schaefer: Analogie zwischen den Verschiebungen und den Spannungsfunktionen 


igli i i i i Gestalt wie die Gleichgewichtsbedin- 
Diese Vertraglichkeitsbedingungen haben praktisch dieselbe talt 
aioe @ 1) bis (2,6). So, wie wir diese nach (3,1) bis (3,6) identisch in den saponin 7, 
Hetriodigtent kénnen wir die Deformationen durch 3 »,Drehungen™ 4, Og» (3 und 3 ,9 Verse a a 
gen“ we vz, also durch 6 voneinander unabhangige VerschiebungsgréBen ausdriicken. Wir be- 
5 9 “2? 2 
kommen die Darstellung: 


— 11 — 050, , a7 %o1 es 0,0, ? (5,8) 


— Hyg = 090 + 2 Ws » — Xg9 = — Wp, (5,9) 
— Yo = 0,03 — 2 Wo 5 —7%, = 003, (5,10) 
Eo, = 00, + 5, Ej4 == Oy, (5,11) 
ban = 22% + — 05, e149 = 01% — a (5,12) 


By = @, — 0y03 + 2 V2» Ba = Wy + 040s - (5,13) 


i ifizi ir mi schi idimensionalen 
Den Vektor v = (v,, V, V3) identifizieren wir mit dem Verschiebungsvektor unseres zwei 
Kontinuums. 


6. Die Analogien zwischen Verschiebungen und Spannungsfunktionen, Deformationen und 
SchnittgréBen. Erfahrungsgema8 sind die Normalenmomente R, und R Null. (Das bedeutet 
streng genommen eine Aussage itiber eine Materialeigenschaft unseres zweidimensionalen Kontinu- 
ums.) Um zu einer Analogie zu kommen, setzen wir auch die Schubdeformationen /, und f, Null. 
Dabei entarteten die Querkrafte Q, und Q, zu Reaktionskraften. Die den Normalenmomenten 
zugeordneten Deformationen 7, und 7, werden uninteressant, weil sie keinen Beitrag zur elastischen 
Energie der Schale liefern, und die Gleichungen (5,10) werden iiberfliissig. Nach (5,13) kénnen wir 
die Drehungen @, und m, durch Verschiebungen ausdriicken: 


1 


@, = 0903 — = > Wo = — 0403 5 (6,1) 
und nach (3,6) wird analog 
2, =2,0,—— G,, 2, = —3, D. (6,2) 
Es stehen nun gegeniiber: 
E11 = Or, , | M,, = 0,9, , (6,3) 
Era = 040, — Ws, M,, = 0,2, — Q, , (6,4) 
Eg, = 020, + Ws, M,, = 0,9, + 25, (6,5) 
1 1 
99 == O05 a ae V3 3 | My. = 0D, ++ ae ®, > (6,6) 
— #11 = 0,0, , Nyy = 0222, , (6,7) 
1 1 
— Hq = O50. + Pe Ny2 = 0082, + fa Q, , (6,8) 
— Hy, = — 0,0, , | Nyy = — 2,2, , (6,9) 
me X20 = 010. . N35 =—— 0,25 . (6,10) 


Dieser kinematisch-statischen Analogie muBten allerdings die Querkraftsverformungen 8, und f, 
geopfert werden, ein Umstand aber, der bei diinnen Schalen nicht ins Gewicht fallt. 


7. Das Elastizitatsgesetz. Unsere Kreiszylinderschale habe iiberall dieselbe Dicke h, und h sei 
so klein gegen a, daB h?/a® gegeniiber 1 immer vernachlassigt werden darf. Das oben betrachtete 
zweidimensionale Kontinuum identifizieren wir mit der Mittelflache der Schale. 

Bei der Kinfiihrung des Elastizitatsgesetzes, des Zusammenhanges zwischen SchnittgréBen 
und Deformationen, haben wir darauf zu achten, daB unsere Analogie nicht zerstért wird. 
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Nach dem Beispiel der ReiSner-Meifnerschen Biegetheorie der axialsymmetrisch beanspruchten 
Rotationsschale setzen wir zundchst 


&,=L (Nix — v Nop) , My, = K (%4, + Y a9) (7,1) 
Egg = L (Nog — Nj3) > Mg, = K (Hg. + ¥ 41) 
mit 
= mis ae E he 
und der Umkehrung 
%, = B(M,,—y Ms) ’ Ny, = D(&y + Y Ego) (7,3) 
Hao = B(M,,—e M,,), Noe = D (E92 + ¥ &3) ‘ 
mit 
ip Eh 
B= ER und D:= fis (7,4) 


(E Elastizitatsmodul, » Poissonsche Querkontraktionsziffer). 
Die elastische Energie je Flacheneinheit der Schalenmittelflache ist 
1 1 1 1 1 
A > (Nyx &11 + Noo 22) + = (Nye + No:) 2 (E12 + &1) + =" (Ni2— Noi) oa (é12 — 21) 
1 M. 1 1 1 1 
ashe (My1 11 + My %0) + 2 (My. ++ Mg;) oF (12 + 1) + > (M,,— M31) > (12 — #21) « 


2 
(7,5) 
Die DeformationsgréBe 


1 1 
pe (E12 — €1) = ” (0,v2 — 0,0,) — ws (7,6) 


ist beim ebenen Kontinuum Null. Die Festsetzung 
ip — ©31 (7,7) 


bedeutet fiir unser gekriimmtes zweidimensionales Kontinuum eine kinematische Zwangsbedingung, 
derzufolge der unsymmetrische Teil des Tensors N;,, zu einem Reaktionsmoment entartet. Gleichung 
(7,7) entspricht den Festsetzungen (6,1), die der Schale die Querkraftsdeformationen verbieten. 


Um bei der Festsetzung ¢,, = €,, unsere Analogie zu wahren, miissen wir nach (6,4) und (6,5) 
M,, = M2, (7,8) 
‘fordern. Auch die Symmetrie des Tensors M;, ist in der Klassischen Plattentheorie vorhanden. 


Wir fiihren dre Abkiirzungen 


As 1 nz 1 
Ni. = a (Ny2 + Noi); 12> > (412 + 21) (7,9) 


ein und kénnen nach den Festsetzungen (7,7) und (7,8) der elastischen Energie die Gestalt geben: 
1 ~ A 
A > (Ny1 €11 + Noe 22 + 2 Myo 12 + Muy 41 + Mae %o2 + 2 Myo %12)- (7,10) 


Unser Elastizitatsgesetz vervollstandigen wir nun durch die ebenfalls dem ebenen Kontinuum 
entlehnten Beziehungen 


Nig = D(1—») 42, My, = K (1—») %o- (7,11) 
Die Gleichungen (7,7) und (7,8) fiihren zu den analogen Gleichungen 
1 
Dae = (0,2 — 05%) » 23 = - (2,2, — 0,P,) . (7,12) 


Damit sind samtliche Deformationen durch die Verschiebungen v,, v,, v; und simtliche Schnitt- 
gréBen durch die Spannungsfunktionen ®,, ®,, D, ausgedriickt. 
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Wir betonen, daB nach wie vor simtliche Gleichgewichtsbedingungen exakt erfillt sind, auch 
die wegen R, = R, = 0 algebraische (2,6) 


' 
Ny2— Noy + Mie =0. (7,13) 


8. Die dualen Variationsprobleme. Zufolge unserer Elastizitatsgesetze (7,1) bis (7,4) und (7,11) 
kénnen wir (7, 10) entweder durch die Deformationen oder die SchnittgréBen ausdriicken. Wir 
erhalten also die beiden Darstellungen 


Ay = Gs {len + #20)? — 20 —») (61s ra — eh] 

+ [bar + 22)? —2(1 —») a1 %a2— Mall} (8,1) 
Ay = ers {{(Max + Mya)®— 2(1 +») (My, Mog — M3,)] 

+ $5 [Mis + Noo)? —2(1 +) (Nix Nos— N3,)I}- (8,2) 


Um zu Variationsproblemen zu kommen, hat man A, durch die 3 Verschiebungen und A, durch > 
die 3 Spannungsfunktionen auszudriicken. Die Variationsprobleme lauten 


f A,[v,, Vg, 3] dF —> Min. , (8,3) 
| A,[®,, By, B,] dF > Min. (8,4) 
Zufolge unserer bisherigen Analogien unterscheiden sich aber die Euler-Lagrangeschen Differential- 
gleichungen beider Variationsprobleme — wie (8,1) und (8,2) erkennen lassen — nur durch das 


Vorzeichen von y, was zu beweisen war. 


9. Das Variationsproblem fiir die Spannungsfunktionen. Nachdem wir das gesteckte Ziel erreicht 
haben, bleibt noch die Genauigkeit unserer neuen Biegetheorie zu erértern. Dazu miissen wir eins 
der beiden Variationsprobleme weiter verfolgen. Wir wahlen (8,4). 

Zunachst eine kleine Vereinfachung in (8,2): Es ist nach (7,13) 

a 1 
Ny No2— Nia = Nyy No2— Ny2 Noi — zy M2. (9,1) 


4 a2 12 


h2 
Wegen a <1 darf das letzte Glied gegen M?, in der ersten eckigen Klammer von (8,2) vernach- 


lassigt werden. Entsprechendes gilt fiir (8,1). 


In (8,2) sind einzufiihren: 


My: + My = M = 0,0, + 2,0, ++ ©,, (9,2) 

Nyy + Nog =4®,—— 0,0, (A = 0,0, + 2,0,), (9,3) 

My; My, — M}, = (8, 2, — 2,0, 0,0,) + * Dy 2,0, — 22, (9,4) 
= 1 @0,—9,0)), (9,5) 

Nix Naa — Nia Nox = (0,2; 22% — 2,2, 2,2) + + 2, (212%, os ay.) : (9,6) 


Die ersten geklammerten Glieder auf den rechten Seiten von (9,4) und (9,6) sind Divergenzen 
und liefern deshalb keine Beitrage zu den Euler-Lagrangeschen Gleichungen eines Variations- 
problems. Den Integranden von (8,4) kénnen wir somit bis auf einen konstanten Faktor schreiben: 


1 
4a =|z MP+ (1 +») Q@—(1 +9) 5% 0,| 


h? fl 1 2 1 
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Nach den Regeln der Variationsrechnung erhalten wir die 3 Differentialgleichungen 


1 2 
2,M — (1 + ») 8,25— (I+) — 2,0, =—™ (1 +») = 0,2, (220, —+ %,), (9,8) 


he 1 2 
QM + (1 +) 0.0, =7, (1 +) 5~ 0,0, (240, —< ®,) ee ets (49, = ax.) (9,9) 


1 1 he 1 aS 
~M—(1+%)— 20, +5 4 (49, — : ay,) — (l-f9)= 2,20, 0), (9,10) 
Aus (9,8) und (9,9) lassen sich M und Q, sehr einfach eliminieren: 


I h? 
AM = (1 +b v) ae 0,0,®, aL 2a 0205 (4%, ee = ayb,) > (9,11) 


(1 +») 42, = — (1 +») — 0,0,.0, + & +t” as (2.0 — 0,) 


h2 1 
+ yyq 91% (4%, eae 20,). (9,12) 


Die stérenden Glieder mit ©, auf den rechten Seiten sind von der GréBenordnung h2/a2 gegen- 
iiber den ihnen entsprechenden Gliedern auf den linken Seiten. Wiirden wir sie einfach weglassen, 
dann kénnten unsere Differentialgleichungen nicht mehr die cof Lésungen 0, OO, @®) fiir die 
am Rande unbelastete Schale enthalten. Diese werden, wie schon im Abschnitt 5 erwahnt, bei 
einem mehrfach berandeten Schalengebiet benétigt. (Im dualen Variationsproblem (8,3) haben 
wir das Anlogon in den oo® Verschiebungszustanden v®, v©, v) der undeformierten Schale.) 
Eingehende Betrachtungen itiber die Rolle der ,,Nullspannungsfunktionen“ des ebenen Kontinuums 
findet man in meiner friiheren Arbeit*. .Nach (3,3) und (3,5) gelten fiir die uns hier interessierende 
Nullspannungsfunktion @Y die Gleichungen 


AO + | OY = 0, (9,13) 


0,09 — 20 =0 
Demnach diirfen wir in (9,11) 


0,0, = — = P, (9,14) 


setzen, das gleiche in der zweiten Klammer auf der rechten Seite von (9,12). In dem vorangehenden 


Gliede darf 
AD, = — — 0,0, (9,15) 


gesetzt werden. Damit aber tritt auf den rechten Seiten von (9,11) und (9,12) nur noch ®, auf. 


Um die Elimination von ®, und @, weiter zu treiben mit dem Ziele, eine einzige Gleichung 
fiir ®, zu erhalten, machen wir die Substitution 


(ig a be ea 


2 
(9,16) 
PD, = iH: Gear he 
und erhalten 
M=~t* 464-9, 
(9,17) 


* Siehe FuBnote 2 von Seite 125. 
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Die Gleichungen (9,11) und (9,12) lauten jetzt 


Ly 1 1 h? 1 : 
LE” AAG = —~ AD, + (L +») = 18,Ps + ay-g 00 (4%, tng ®,), (9,18) 
1 1 h? 3+ 4 1 


Die Gleichungen (9,16) fiihren wir ean in (9,10) ein, wenden den Operator 44 auf diese Gleichung an 
und ersetzen die hier auftretenden AAG und AAH durch (9,18) und (9,19). Damit ist die Elimination 
von ®, und @, vollzogen, und es bleibt nach einiger Rechnung (h?/a? <_1) 


o8 AAAAD, + 208 A 5a (Ay + 2S) 
2 a2@, 12a ad 
ob a‘ i (4%, + 2 =) ++ ie (1 eras y?) ae | a 0 é (9,20) 


Fiir v,, die Verschiebung in Richtung der Schalennormalen, besteht dieselbe Differential- 
gleichung; denn (9,20) ist invariant gegeniiber der Vertauschung des Vorzeichens von y. Die GroBe 
@, entspricht aber der Airyschen Spannungsfunktion der Scheibe und v, der Durchbiegung der 
Platte. Somit treffen wir auch bei der Schale auf die Timpe-Wieghardtsche Analogie. 


Kin Priifstein fiir die Genauigkeit einer Biegetheorie des Kreiszylinders ist die charakteristische 
Gleichung der Partikularlésung 


*%4 


@, =e * cosn (n 0, 1, 2;3,61%)5 (9,21) 
die bei uns 
78 — 4 n2 78 + |6 n? (x? —1) +2 (1—va)| 28 
—n? (4 n4— 8 n? + 3) 22 + nt (n?— 1)? = 0 (9,22) 


lautet. 

In der nach ganz anderen Gesichtspunkten aufgebauten Biegetheorie der Kreiszylinderschale 
von W. Fligge!, einer Theorie groBer Genauigkeit und logischer Geschlossenheit, die allerdings 
keine Spannungsfunktionen benutzt, gibt der Ansatz 


foe 


vy =e * cosn 2 (n = 0,1, 2,3,...) (9,23) 
die charakteristische Gleichung 


#8 —2 (2 n? —») 18 + [6 n4(n?— 1) + 4 a — pay] 28 


—n® [4 nt — (8 — 2) n? + 4 —29] 22 + nf (n2?@— 12 = 0. (9,24) 


Der Vergleich mit (9,22) zeigt in den Koeffizienten von A° und J? Unterschiede, die wegen 0 <» < 


nur bei n = 0 und n = | erwahnenswert sind. Aber dann macht der ttberwiegende Faktor von /4 
(GréBenordnung 10) gerade die Glieder mit 7° und 72 bedeutungslos. Es kann somit festgestellt 
werden, da unsere Analogie nicht mit einem Verlust an Genauigkeit erkauft zu werden brauchte. 


10. Randbedingungen bei belastetem Rande. Bekanntlich kénnen den 3 Verschiebungen Vy, Vay Vg 
4 Bedingungen am Rande auferlegt werden. Zufolge unserer Analogie kénnen die 3 Spannungs- 
funktionen auch dann nur 4 Bedingungen erfiillen, wenn am Rande die 5 Komponenten der Be- 
lastung vorgeschrieben sind. In diesem Falle haben wir so vorzugehen: Nach (4,8) bis (4,13) werden 
die Randwerte der 6 Spannungsfunktionen Oe ere, @, durch Integration ermittelt. 
Anschliefend kann nach (6,2) die Normalableitung von ®, berechnet werden. Somit haben wir 
4 Randbedingungen fiir D,, ®,, ®, gewonnen, und die Randwerte von Q,, Q, und Q, miissen im 
Laufe der weiteren Rechnung ignoriert werden. Im Endergebnis der Schalenberechnung kann man 


aus den bekannten ®,, ®,, ©, die Q,, Q,, Q, berechnen und dann riickwArts aus (4,8) bis (4,13) 


1 W. Fligge, Statik und Dynamik der Schalen, S. 152. 2. Aufl. Berlin-Gottingen-Heidelberg 1957. 
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eine Randbelastung ermitteln, die zwar im allgemeinen nicht mit der vorgeschriebenen iiberein- 
stimmt, ihr jedoch an jedem Bogenelement der Berandung statisch aquivalent ist. 


Die bei mehrfach berandeten Bereichen auftretenden Komplikationen sind bereits oben be- 
sprochen worden. 


li. Zusammenfassung. Fiir die Kreiszylinderschale wurde eine Biegetheorie aufgestellt, in der 
die Gleichgewichtsbedingungen (unter Voraussetzung der Symmetrie des Momententensors M;,) 
durch drei Spannungsfunktionen ©,, ©, ©, exakt erfiillt sind. Bei der Definition der Deformations- 
gréBen und der Einfiihrung der Elastizitatsgesetze war die ReifSner-MeiBnersche Theorie der 
symmetrisch belasteten Rotationsschale das Vorbild. Die drei Differentialgleichungen fiir die Ver- 
schiebungen v,, v,, v; unterscheiden sich von den drei Differentialgleichungen fiir die Spannungs- 
funktionen ®,, ®,, ®, formal nur im Vorzeichen der Poissonschen Querkontraktionsziffer y. Die 
beiden Differentialgleichungen achter Ordnung, die man nach Eliminationsprozessen sowohl fiir 
v3 als auch fiir ©, erhalt, unterscheiden sich nicht mehr voneinander. So trifft man bei der Zylinder- 
schale die Timpe-Wieghardtsche Analogie zwischen Durchbiegung v, der Platte und Airyscher 
Spannungsfunktion ®, der Scheibe wieder. 

Es konnte ferner gezeigt werden, daB unsere neue Biegetheorie der bekannten Fliiggeschen 
Theorie an Genauigkeit nicht nachsteht. 

Es ist wohl nicht zu bezweifeln, daB auch bei Schalen beliebiger Gestalt unsere Analogie vor- 
handen ist. Sie scheint uns wertvoll als Ordnungsprinzip inmitten der Fiille von Gleichungen, die 
nun einmal zu einer Schalentheorie gehéren. 

Die Formulierung des Schalenproblems mit Hilfe der drei Spannungsfunktionen @,, ®,, D; wird 
sich immer dann empfehlen, wenn die Randbelastung vorgegeben ist. Denn dann lassen sich die 
Randbedingungen in den Spannungsfunktionen iibersichtlicher formulieren als in den Verschiebun- 
gen. Auch die Gewi®heit, daB selbst durch radikales Streichen lastiger Glieder in den Differential- 
gleichungen der Spannungsfunktionen die Gleichgewichtsbedingungen nicht verletzt werden, mag 
manchem Rechner angenehm sein. 


(Eingegangen am 22, Juni 1959.) 
Anschrift des Verfassers: Prof. Dr.-Ing. H. Schaefer, Braunschweig, BernerstraBe 9. 
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Eigenschwingungszahlen zusammengesetzter Schwingungs-Systeme 
Von E. Hiibner 


1. Einleitung. Der Einsatz moderner Rechenautomaten erméglicht es in neuerer Zeit, Schwin- 
gungsgebilde mit einer hohen Zahl von Freiheitsgraden zu berechnen. Besonders pradestiniert 
fiir diese Verfahren sind stab- bzw. plattenférmige Gebilde, bei denen der Berechnungsgang mit 
Hilfe des Matrizenkalkiils auf einfache, stets wiederkehrende Rechenoperationen (Matrizenmulti- 
plikationen) zuriickgefiihrt werden kann!. Dabei lassen sich auch erschwerende Randbedingungen 
(beispielsweise elastisch gelagerte Stiitzen, oder in besonderen Fallen die elastischen Eigenschaften 
eines vorhandenen Fundamentes) in die Rechnung einbeziehen. 

Bei komplexen, aus kontinuierlichen Teilgebilden zusammengesetzten Systemen ist die An- 
wendung dieser Berechnungsverfahren im allgemeinen nicht ohne weiteres méglich. Hier kann eine 
Verallgemeinerung des bekannten Verfahrens® ? der Aufspaltung eines Schwingungssystems in ein- 
facher aufgebaute, und damit leichter zu berechnende, Teilgebilde wesentliche Vorteile bringen. 
Wird dabei die Zerlegung so vorgenommen, daB die entstehenden Teilgebilde stab- oder platten- 
férmig sind, wie dies bei den meisten technisch wichtigen Anordnungen méglich ist, dann kénnen 
die vorerwahnten Verfahren zur Berechnung der Eigenwerte und Higenformen der Teilgebilde an- 
gewendet werden. Die Eigenwerte des Gesamtsystems lassen sich dann in einem anschlieSienden 
Rechnungsgang, in dem die vorliegenden Verknipfungsbedingungen zum Ansatz kommen, aus- 
rechnen. ; 

Wesentlich ist aber, daB der Ablauf dieses Rechnungsganges méglichst rationell gestaltet wird. 
Dafiir eignet sich in ganz besonderer Weise die Matrizenrechnung; durch deren Anwendung werden 
die zur Verknipfung der Teilgebilde-Schwingungsgleichungen erforderlichen algebraischen Rechen- 
operationen weitgehend ihres abstrakt-mathematischen Gehalts entkleidet. Der Ablauf der 
Rechnung wird dabei in eine rein schematisch durchfiihrbare ,,Verzifferungsaufgabe“ iiber- 
gefiihrt, die von angelernten Kraften — ein wesentlicher Gesichtspunkt im Zeichen des heutigen 
Mangels an ausgebildeten Fachkraften — durchgefiihrt werden kann. 


2. Frequenzgleichung zusammengesetzter Systeme. a) Einfach gekoppelte Gebilde. Wir 
betrachten die beiden in Abb. 1 symbolisch dargestellten Gebilde 1 und 2, die durch gleichfrequente 
harmonische Wechselkrafte mit der Amplitude Q, bzw. Q, zu erzwungenen Schwingungen mit den 


1 2 
3 i 
% % 


Abb. 1. Symbolische Darstellung Abb, 2. Einfach gekoppelte wechselkrafterregte Gebilde. 
wechselkrafterregter Schwingungsgebilde. 


Amplituden q, bzw. q, angeregt werden. Unter der Voraussetzung, dafs die Gebilde linear und un- 
gedimpft sind, gelten die beiden Amplitudenbeziehungen 


© =4,Q,, q2 = d,Q, (1) 


mit den im allgemeinen frequenzabhangigen Proportionalitatsfaktoren d, und d,, die man als 


dynamische EinfluBzahlen (im angelsachsischen Schrifttum als Receptance bzw. Admittance ge- 
nannt) bezeichnet. 


* Siehe etwa H. Fuhrke, Ing.-Arch. 23 (1955), S. 329, und E. Pestel, G, Schumpich und S. Spierig, VDI- 
Berichte 35 (1955), S. 11. 


foe * Bishop, Proc. Inst. Mech. Engrs, 169 (1955) S. 1031, T.G. Sofrin, Journ. Aeron. Sciences, 13 


3 R. E. D. Bishop and D. C, Johnson, Vibration Analysis Tables, Cambridge 1956. 
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Werden nun die beiden Gebilde gemaB Abb. 2 an den bezeichneten Punkten starr miteinander 
gekoppelt, dann miissen die beiden Gleichungen gelten 


: aia de = 9) Q+Q=Q. (2) 
Mit (1) folgt daraus 


(a ta)1=@- (3) 


Die Amplitude der Erregerkraft Q ist frequenzabhangig. Im Resonanzfall ist Q = 0; damit folgt 
die Resonanzbedingung unmittelbar zu 


L 1 
eae bzw. d,+d,=0. (4) 


Diese Gleichung stellt die Frequenzgleichung des gekoppelten Gebildes dar, aus der sich die 
Eigenfrequenzen ausrechnen lassen. 


Abb. 3. Schwingungsgebilde mit mehreren Erregerkriften. Abb. 4. Mehrfache Kopplung zweier Schwingungsgebilde. 


b) Mehrfache Kopplung. Greifen an einem ungedémpften linearen Gebilde mehrere, bei- 
spielsweise n, gleichfrequente harmonische Erregerkrafte an (siehe Abb. 3), dann werden die 
Schwingungsamplituden der n Kraftangriffspunkte durch die n linearen Gleichungen 


qi = 4;,Q, + din Qe ++++ din Qn (« = 1,2,...n) (5) 
beschrieben, die sich durch die Matrizengleichung 
q= 2-2 (6) 


ausdriicken lassen. D ist dabei die Matrix der dynamischen HinfluBzahlen, die wir kurz De- 
formationsmatrix nennen wollen. 

Koppelt man nun zwei Gebilde dieser Art in der in Abb. 4 symbolisch angedeuteten Form, dann 
laBt sich fiir jedes Teilgebilde eine derartige Gleichung anschreiben. Mit den Setzungen 


nh Q; 
G1 =| 4% |> Oy, =| QO, |}, 
3 Qs (7) 
q4 Q, 
Go =1 5 }> De. =| Qs |> 
I 0, 
erhalten wir dann 
q = 2,4 : (8) 
qo = Dz De 


Bezeichnet man die Amplituden der von auBen her einwirkenden (gleichfrequenten) Erreger- 
krafte mit Q;, Qu, Qu, die Schwingweg-Amplituden der Kraftangriffspunkte mit qr, qi, qi» 
und setzt man weiterhin 

q1 Qr 
q=l|qu |, N%=1On |, (9) 


quit Quit 


136 E. Hiibner: Eigenschwingungszahlen zusammengesetzter Schwingungs-Systeme Ingenieur-Archiv 


so lassen sich die aus der Abb. 4 ablesbaren Verkniipfungsbedingungen durch die beiden Matrizen- 
gleichungen 


eee (10) 
, Dy + DO, = 0 
ausdriicken. 
Resonanz liegt wieder dann vor, wenn {) = 0, also gemaB (10) 
Dy =o mee D. (lia) 

ist. Wegen q, = q. = q ergibt sich damit aus (8) die Resonanzbedingung 

D, O — — Dz Oi (11b) 

oder 
(D, + Ds) O,=0. (12) 


Diese Bedingungsgleichung besitzt dann nichttriviale 
Lésungen fiir D, bzw. O., wenn 


det (D, + D,) =0 (13) 


ist. Die Auswertung dieser Determinantenbeziehung 
liefert die Frequenzgleichung des zusammengesetzten 
Abb. 5. Mehrere mehrfach gekoppelte Gebilde Gebildes, aus der die Eigenschwingungszahlen be- 
rechnet werden kénnen. 


c) Kopplung mehrerer Gebilde. Werden mehrere Gebilde miteinander gekoppelt, dann 
bestehen analoge Beziehungen. Betrachten wir beispielsweise das in Abb. 5 symbolisch dargestellte, 
aus drei linearen Teilgebilden bestehende System. Fiir jedes Teilgebilde gilt eine Gleichung der 
Form (6), die sich, umgeformt, durch 

OH =D1q (14a) 
darstellen 148t. Die Kehrmatrix einer Deformationsmatrix hat die Bedeutung einer Steifigkeits- 
matrix. Wir setzen deshalb 


© =D = Matrix der dynamischen Steifigkeitszahlen (14b) 
und kénnen damit die fiir die drei Teilgebilde geltenden Beziehungen durch 
Oy, = G41, 
Dy = Ge qe, (14¢) 
Ds = Gs qs 


ausdriicken. Die den Amplituden der Schwingungswege bzw. Erregerkrafte durch die Kopplung 
auferlegten Verkniipfungsbedingungen lassen sich unmittelbar dem Kopplungsschema (Abb. 5) 
entnehmen: 


n1=%=—47=— 41> 


q2 = 45 = Qs (15) 
= We = Gr > 

Q, +0; +07=01. 
Q. +0; = Qn, (16) 
0, +0, = Om - 


Geht man nun mit (14c) und (15) in (16) ein, so erhalt man ein System linearer Gleichungen, das 
mit den Setzungen 
qt Qi 


q=(qn |, O=10n |, ALT) 
qtr : Qin 
durch die Matrizengleichung 


D=Cq (18) 


dargestellt werden kann. ( ist dabei die auf die Kopplungspunkte bezogene Matrix der dynamischen 
Steifigkeitszahlen des Gesamtsystems, die wir kurz Steifigkeitsmatrix nennen wollen. 
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Resonanz tritt bei Q = 0, also 


CqG=0 (19) 


auf; nichttriviale Lésungen gq bestehen dabei nur dann, wenn 
det © =0 (20) 


wird. Durch ist di i r si i i 

< ch (20) ist die Frequenzgleichung gegeben, aus der sich die gesuchten Eigenfrequenzen des 
gekoppelten Gebildes berechnen lassen. 
* oe ee Rechnungsgang erfordert, insbesondere dann, wenn ein Gebilde mit einer 
has zahl von Kopplungen vorliegt, einen erheblichen Anteil algebraischer Berechnungsarbeit. 

er Aufwand an geistiger Arbeit zur Bewaltigung dieser Aufgabe kann, wie schon eingangs 
ae ; ‘ ; : : ; : 
trwahnt, durch Anwendung dex Matrizenrechnung in Form einer Matrizen-Transformation be- 
merkenswert verringert werden. 

Dieses Verfahren, das beispielsweise in der Elektrotechnik bei der Berechnung elektrischer Netz- 

werke allgemein angewendet wird!, soll hier kurz erlautert werden. 


3. Verkniipfungsmatrix und Steifigkeitsmatrix. Wir gehen von dem vorstehend behandelten, 
in Abb. 5 symbolisch dargestellten gekoppeltem Gebilde aus, und ziehen zunichst die fiir die drei 
Einzelgebilde geltenden Matrizengleichungen (14c) zu einer einzigen Matrizengleichung 


Oy qy 
Ds qe (21a) 
Ds qs 
zusammen, die mit den Setzungen” 
Q, | 
Dy Q, qh I C, aa 
B=(O.J=|G)> G=(qa)=|a|- C= [Ge] | (21b) 
Os : q3/ : es le 
Q, {7 
in die Form 
N=Cq (21c) 


gebracht werden kann. 

Wir stellen uns nun die Aufgabe, eine Matrizen-Transformation aufzusuchen, mittels der die 
durch einfaches Nebeneinandersetzen der Teilmatrizen gefundene Matrix © in die dynamische 
Deformationsmatrix © des gekoppelten Gebildes [siehe Gleichung (18)] tibergefiihrt wird. Zu 
diesem Zweck schreiben wir die zwischen den Elementen der Spaltenvektoren g und q bestehen- 
- den, bereits in (15) dargestellten Verkniipfungsbeziehungen in einem etwas abgednderten Schema 
an: 

a= 7 ee Us 
= 0+ gu +0, 
q3 = I Sp Oseter O's 
G0, 420 din (22a) 
I =0 +qu+?, 
nde = Co 0 + aus 
depot er 


Dieses Gleichungssystem la8t sich mit der Setzung (17) durch die Matrizengleichung 
C= 2G, (22 b) 


die eine auf den Spaltenvektor q angewendete Matrizen-Transformation darstellt, ausdriicken. 


1 G. Kron, Tensor Analysis of Networks, New York 1949. 
2 Die leeren Felder in der Matrix © sind mit Nullen besetzt zu denken. 
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Die Transformationsmatrix ®, die als Verkniipfungsmatrix bezeichnet wird, kann unmittel- 


bar aus (22a) abgelesen werden ?: 


(22c) 


Eine kurze Betrachtung der zwischen den Erregerkraftquellen und dem Koppelgebilde aus- 
getauschten Energie fiihrt zu der Feststellung, daB die von den Erregerkraften in das gekoppelte 
System eingespeiste Blindleistung (mit w = Erregerkraft-Kreisfrequenz) 


== 1 
Ngi = 5 O (Qr gr + Qu qu + Qur git) » (23a) 
bzw. in Matrizenschreibweise (mit (17)) 
Ne = : oy Dd, (23b) 


wobei qj’ in der iiblichen Schreibweise den zu q transpornierten Vektor bedeutet, gleich der Summe 
der in die Teilgebilde eingespeisten Blindleistungen 


1 
Np => (Qi G1 + Qs G2 + +++ Q7 47); (24a) 
oder mit (21b) 
Nar =o (24b) 
sein muB. 
Daraus folgt sofort 
qa (25) 


Nun wenden wir auf den durch (21b) definierten Spaltenvektor q die in (22b) dargestellte Trans- 
formation an, erhalten 


ON =COq (26a) 
und multiplizieren diese Gleichung von links her mit der Transponierten @’ : 
PDO=PCOG. (26b) 
Mit der zu (22b) transponierten Form 
Gee he (27) 
gehen wir in (25) ein 
7Oo=q OD. (28a) 
Mit (26a) folgt daraus 
. qQX=q0'CGq, (28b) 
oder 
N= PS CO.G. (28c) 


Kin Vergleich dieses Ergebnisses mit (18) liefert sofort die gesuchte, auf die Matrix © anzuwendende 
Matrizen-Transformation 


C=O CO. (29) 


Die Verkniipfungsmatrix ® ist denkbar einfach aufgebaut. Sie besteht nur aus Kinsen und 
Nullen. Damit wird auch die in (29) dargestellte Transformation leicht durchfiihrbar und lauft 


letzten Endes auf eine Summenbildung der Elemente der Teilmatrizen €; in den Elementen von € 
hinaus. 


1 Die leeren Felder in ® sind mit Nullen besetzt zu denken. 
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Die Durchfihrung der Rechnung kann durch eine geschickte Anordnung der die Verkniipfungs- 
bedingungen beschreibenden Gleichungen unter Umstanden wesentlich vereinfacht werden, wenn 
es namlich gelingt, die Verkniipfungsmatrix ® so aufzubauen, daB méglichst viele Einheitsmatrizen 


in ihr enthalten sind (siehe Abb. 9e). 


© ist eine quadratische Matrix, deren Ordnungszahl n (= Anzahl der Kopplungsstellen) véllig 
unabhangig ist von der Zahl der Freiheitsgrade der beteiligten Teilgebilde. 

Beispielsweise ist fiir das in Abb. 5 dargestellte System n = 3; dabei kann es sich bei den beiden 
gekoppelten Teilgebilden um Gebilde mit beliebig vielen Freiheitsgraden handeln (im Extremfall 
unendlich viele, wenn es kontinuierliche Gebilde sind). 

AbschlieBend ware noch zu erganzen, daB die Amplituden q, bzw. Q; verallgemeinerte Koordi- 
naten bzw. verallgemeinerte Krafte sind; d:h., sie kénnen sowohl Wege bzw. Krafte, als auch 
Drehwinkel bzw. Drehmomente (durch passend gewahlte Bezugslangen homogen gemacht) darstellen. 

Die durch (20) zum Ausdruck gebrachte Resonanzhedingung wird sich nur in den seltensten 
Fallen in allgemeiner Form geschlossen auswerten lassen. ZweckmaBigerweise wird man hierzu 
das bekannte RestgréBenverfahren anwenden, indem man eine passende Kreisfrequenz w wahlt 


und den dazu gehérenden Wert det © ausrechnet. Man wiederholt die Rechnung mit abgeanderten 
Werten von w, trigt die Werte det € in einem Koordinatensystem tiber ~ auf und bestimmt die 
Nullstellen der sich ergebenden Funktionslinie. 
Bei der Ausrechnung von det © kann man sich mit Vorteil, insbesondere bei groBen Ordnungs- 
zahlen n, der bekannten, aus dem Gaufschen Algorithmus hergeleiteten Verfahren! bedienen. 


4, Bestimmung der dynamischen EinflufBzahlen. a) Kontinuierliche Gebilde. Bei der 
vorliegenden Betrachtung wollen wir uns aus Raumgriinden auf reine Biegeschwingungsgebilde 
beschranken. Die hierbei gewonnenen Ergebnisse lassen sich analog auf Langs- und Drehschwin- 
gungsgebilde, bzw. auf Gebilde mit gekoppelten Biege- und Verdrehungsschwingungen, iber- 

tragen. 
: Greifen an einem Biegestab zeitlich harmonisch veranderliche, verteilte Belastungskrafte und 
-Momente gleicher Frequenz an, deren Amplitudenfunktionen p(x) und M(x) seien, dann kann 
die Gleichung fiir die Schwingungsform in bekannter Weise quellenmaBig (mit u als Parameter) 
~ durch : : 
ye) = f Gols w) [plu) +2 gw) y(u)] du + [Gru(, uw) M(u) du (30) 
i 


ausgedriickt werden (die hierbei zugrundegelegten positiven Zahlrichtungen der GréBen y(x), p(x) 
und M(x) sind aus Abb. 6 zu ersehen). Die in obiger Gleichung verwendeten GréBen bedeuten: 


Gp(x, u) die Greensche Funktion fiir die Belastungskrafte, |x |r Qe 
Gy(x,u)die Greensche Funktion fiir die belastenden Biegemo- 
mente, 
o(x) die Massenbelegung des Biegestabes, 2 
A= ow (m = Erregerkraft-Kreisfrequenz). Abb. 6. Biegungstrager mit Zahlpfeilen 


Aufgrund des Entwicklungssatzes laBt sich die Amplitudenfunktion y(x) als Summe 
y(x) = D &, *y(x) » (30a) 
eal 


mit den k. als Entwicklungskoeffizienten (die man als Normalkoordinaten oder generalisierte 
fe erdinaien bezeichnet), und den *y(x) als Kigenfunktionen, ansetzen. Mit diesem Ansatz gehen 
wir in (30) ein; dann multiplizieren wir diese Gleichung mit o(x)’y(«) und erhalten daraus nach 
einer Integration iiber die gesamte Lange des Gebildes 


[ (Sk ve) ete) 910) de = ff Cals 1) 910) ole) de] ne) de 


A | [Shs f Gols 1) ele) (uo) de + fff Gals) ols) 9) es] MU) de (300) 
§ stot 


Unter Beachtung der Orthogonalitatsrelationen 


ay J ol) ‘y(x) ya) de =0  (r #5), (30) 


1 Siehe ene R. Zurmiihl, Matrizen, Seite 62, 2. Aufl., Berlin, Gottingen, Heidelberg 1950. 
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und der fiir den Fall der r-ten Eigenschwingung geltenden Gleichung 


‘y(u) = A, J Gp(x, u) o(x) "y(x) dx (30d) 


sowie der aus der Elastizitatslehre her bekannten Beziehung 
Gy(x, u) = Gp(u, x) , (30e) 
mit Gp(u, x) als EinfluBfunktion fiir die Neigung der elastischen Linie eines Biegetragers infolge 


einer Belastung durch eine Einzelkraft, erhalten wir nach teilweiser Integration schlieBlich, wenn 
anschlieBend an Stelle des Parameters u wieder x gesetzt wird, 


ky [ ole) 9%) de = 7 | ye) pla) de +7 e ole) 9) de 5, ['9’) M@) dx. Gof) 
i i i i 
Durch den hochgestellten Akzent bei ‘y’(x) soll hier in bekannter Weise die erste Ableitung nach 
der unabhangigen Veranderlichen x gekennzeichnet werden. 
Mit den Setzungen 


Ll, = { 0(x)'y?(x) dx als r-te generalisierte Masse, 
1 
"Kp = { "y(x) p(x) dx als r-te generalisierte Kraft, herriithrend von den Belastungskraften, + (30g) 
i 


"Ky = | "y’(x) M(x) dxals r-te generalisierte Kraft, herriihrend von den Belastungsmomenten, 
i 


erhalten wir nach Auflésung nach den Entwicklungskoeffizienten 


ee Sa its 5 eS 
7 fr Gr —A) (r¥ = 1, 2,3,...). (31) 


Gehen wir mit diesen Lésungen in (30a) ein, so ergibt sich die gesuchte Amplitudenfunktion y(x) 
und ihre Ableitung nach x zu 


(32) 


und a5 


Fur den Fall, da8 das Gebilde nicht von stetig verteilten Kraften bzw. Momenten, sondern von 
n harmonischen gleichfrequenten Einzelkraften bzw. m Einzelmomenten erregt wird, nehmen die 
in (30g) dargestellten Ausdriicke fiir die generalisierten Krafte, wenn wir die Lagekoordinaten 
der Kraft- bzw. Momentenangriffspunkte mit x, bezeichnen, die Formen 


Kp = > Prva) 


m 


‘Khiy = 2 My, "y' (xz) 


(r= 1,°2, 35...) (33) 


an. 
Die in einem durch die Lagekoordinate x; gekennzeichneten Punkt herrschende Schwingungs- 
amplitude erhalt man durch Einsetzen von (33) in (32). Werden dabei die Summanden in den ent- 
stehenden Doppelsummen umgruppiert, so folgt, wenn gleichzeitig durch Einfiihrung einer passend 
zu wahlenden Bezugslange /, auf homogenisierte Krafte iibergegangen wird, 


as Sy (x4) y(xj) NO Me) So ty'(xx) ty(xp) 
UO) TO ey tT Oh Diane See 


Si 


Analog dazu erhalt man die Neigung der Kigenschwingungslinie an der Stelle x; (in homogener 


Form) zu 
Ae 3 rye) Fy'(i) NOV Me yy Wr y’(ou) Ty'(Hp) 
oy = >, Py, a ee pian tre (34b) 


Ein Vergleich mit (5) zeigt unmittelbar die physikalische Bedeutung der bei den P, baw. M,/1y 
stehenden Summen: sie stellen dynamische EinfluBzahlen dar. 
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Ist x; die Lagekoordinate jenes Punktes, in dem die Schwingungsamplitude y(x;) herrscht und x, 


die des Angriffspunktes einer Erregerkraft P;, bzw. M,/I, dann ergibt sich die auf diese Punkte 
bezogene dynamische EinfluBzahl zu 


i= > eo) (fir die Kraft P,) ; 
baw = (35 a) 
dig =e a ve a (fiir die homogenisierte Kraft M;,/I)) . 
Analog erhalten wir aus (34b) die durch 
pee Ee eet 8d, = we (35) 


definierten, auf die Punkte x; und x; bezogenen dynamischen EinfluBzahlen fiir die Neigung der 
Eigenschwingungslinie: 


Op = eas 
r=l1 


Une 
(350) 
Alera ae 
Aus den entwickelten Gleichungen kénnen die bekannten Reziprozitatsrelationen 
dj, == dij ’ dj, = dj, ’ dj, =e d',,j ’ (35 d) 


leicht abgelesen werden. 

Im allgemeinen konvergieren die Summen auf den rechten Seiten der obigen Gleichungen ziem- 
lich rasch, so daB man sie, abhangig vom Grad der geforderten Genauigkeit, bei Gliedern verhaltnis- 
maBig niedriger Ordnungszahl, diese werde mit f bezeichnet, abbrechen kann. Physikalisch be- 
deutet dies, da®B das der Berechnung vorliegende Gebilde durch ein aus diskreten Massen und 
Federn bestehendes Ersatzgebilde von f Freiheitsgraden approximiert werden kann. 


«) Ersatzgebilde. Fir das vorstehend erwahnte approximierende Ersatzgebilde von f 
Freiheitsgraden lassen sich f Eigenvektoren 


ry (x) 
y= ye) (enh 9y agehh (36a) 
ry (iq) 


angeben, deren Elemente die Schwingweg-Amplituden der n Bezugspunkte x, in den f Eigen- 
formen sind, die man zu einer (nf) Matrix (Matrix der verkiirzten Eigenvektoren) zusammen- 
ziehen kann: 


re (OT AY oT) (36b) 


Dabei muB aus physikalischen Griinden die Anzahl der (an sich nur durch die geforderte Genauig- 
keit festzulegenden) Freiheitsgrade f stets gréBer als (oder zumindest gleich) nm sein. 

Mit (36b) kann beispielsweise die fiir die Krafte P, geltende Deformationsmatrix Dp, deren 
Elemente durch (35a) bestimmt sind, als Matrizenprodukt in der fiir das numerische Rechnen 
bequemen Form 


Dp = "| : 


*/ 36 
ne=a|* (36¢) 


angeschrieben werden. 


Es ist dabei 
Dp eine (nn) Matrix, 
X* eine (nf) Matrix (X*’ ist ihre Transponierte), | ; 
[ ] eine diagonale (ff) Matrix, die also mit Nullen in den nichtdiagonalen Elementen besetzt ist. 
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fp) Steifigkeitsmatrix. Die Matrix der dynamischen Steifigkeitszahlen ist gemals (145) die 
Kehrmatrix der Deformationsmatrix. Fiir den meist immer vorliegenden Fall f +n ist eine ge- 
schlossene Darstellung von & in allgemeiner Form aus (36c), was fiir die weitere Rechnung sehr 
niitzlich ware, wegen der Rechteckform von X* leider nicht méglich. Man mu deshalb zu jedem 
gewihlten Wert A die Deformationsmatrix berechnen und daraus erst die Kehrmatrix bestimmen. 
Die dadurch bedingte Erschwernis des Berechnungsganges ist jedoch im allgemeinen nicht allzu- 
groB, da fiir gewéhnlich die Anzahl der an einem Teilgebilde vorhandenen Kopplungsstellen, 
und damit die Ordnungszahl n der Deformationsmatrix, immer verhaltnismaBig niedrig ist. 


b) Sonderfalle. Die Vorteile des beschriebenen Verfahrens treten dann besonders augen- 
fallig in Erscheinung, wenn die Zerlegung so vorgenommen werden kann, daB Teilgebilde entstehen, 
deren Eigenwerte und Eigenschwingungsformen sich leicht berechnen lassen. Von den vicious: 
Méglichkeiten seien nachstehend die wichtigsten Falle kurz behandelt. 


a) Homogene stabférmige, biegungsbeanspruchte Gebilde. Das Schwingungs- 
verhalten dieser Gebilde wird in fast allen Lehr- bzw. Handbiichern behandelt, so daf die zu den 
verschiedenen Randbedingungen gehérenden Lésungen leicht aufgefunden werden kénnen. Die 
Berechnung der durch (30g) definierten generalisierten Massen gestaltet sich fiir homogene Biege- 
stabe besonders einfach gemaB1 


= My 270-9" +90), ee 


wobei M die Masse des Stabes, und "y(1) die Ordinate der r-ten Eigenschwingungsform an der 
Stelle x =I ist. Die hochgestellten Akzente bei "y(I) bedeuten hierbei Ableitungen nach den 
Argumenten o, 2 x/I der in den Eigenfunktionen (siehe Tabelle 1) vorkommenden transzendenten 
Funktionen, den bekannten Rayleighfunktionen. 

Beispielsweise ist fiir den frei-freien Biegestab "y’’(l) = 0, "y 
Masse durch 


ttt 


(1) = 0, so daB die generalisierte 


br = VAD) (37b) 


gegeben ist. 

Der Vollstandigkeit halber sind in Tab. 1 die Eigenwerte, die Eigenschwingungsformen und die 
generalisierten Massen homogener Biegestabe fiir die drei wichtigsten Lagerungsarten bis zur 
fiinften Ordnung angegeben. Zu bemerken ware hierzu noch, daB die in den Lésungen vorkom- 
menden transzendenden Funktionen (Rayleighfunktionen) tabelliert vorliegen.* 


f) Starre Gebilde. Ist ein Teilgebilde sehr viel steifer als die iibrigen Teile des Gesamt- 
systems, dann kann man es, ohne einen allzu groBen Fehler zu begehen, als vollkommen starr be- 
trachten. 

Wir gehen von einem starren plattenférmigen Korper aus, der frei von Bindungen an das 
feste Bezugssystem sei. Er besitze drei Freiheitsgrade, zufolge deren er Bewegungen in der y-Achse 
und Drehungen um die x-Achse bzw. z-Achse eines dreiachsigen raumfesten Koordinatensystems 
(mit dem Ursprung im Schwerpunkt des Kérpers in der Ruhelage) ausfiihren kann. 

Unter diesen Voraussetzungen besitzt der Kérper nur den einzigen, allerdings dreifach zu 
zahlenden, Eigenwert A = 0, mit den zugehérigen EKigenfunktionen 


! et y( xe) == 1, O2y (25 2) == wile OSe(a, 2) == 2)... (38 a) 
und J, als passend zu wahlender Bezugslange. Die zu diesen Eigenfunktionen gehorenden generali- 
sierten Massen ergeben sich gemaB (30g) zu 

Moi =M, Hos = 9,/I5 ; Hos = 9,/Tp , (38b) 
mit M = Masse des Kérpers ; 


0, = Massentragheitsmoment des Korpers, bezogen auf die z-Achse 
O,, = desgl., bezogen auf die x-Achse des Koordinatensystems. 


u Lord Rayleigh, Theorie of Sound, Bd. 1, Seite 256, London 1894. 
2 Siehe K. Hohenemser u. W. Prager, Dynamik der Stabwerke, Seite 323, Berlin 1933, sowie die unter 
FuBnote 3 auf Seite 134 angefiihrte Arbeit, in der auch die Eigenschwingungsformen mit ihren Ableitungen 


fiir homogene Biegungsstabe bis zur Ordnungszahl 5 fiir die technisch wichtigsten Lagerungsarten tabelliert 
angegeben sind. 
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Tabelle 1 


ry(x) = A,( B(x) — C, D(x)) 


eee LOR 
ee RS SD Sa ee 
| Ls Tare My C, | A, Schwingungsform 
a 
Frei-freier Biegetrager 0, | 0 | M ane 1/2 G (a ai 
il _ el 
z=0 z= 05 0 M/3 one So 0 7 
F 1 | 4,73004 | M 0,982 502 1 jo aS eee 
Oo, % | | 
Bix) = Goj—j—x + cos x | 2 | 7.85320] M | 1,000777 | 1 iPS ee 
| | 
=O; Grau z | ee a 0,999:967 pS 
D(x) = Gin he eer sin Wee | 
| 4 | 14,137 2 M 1,000 002 1 ONS 
c _ Sei orm — cos or | ‘ 
r=" Gino, 7 — sino, x 5 172788 | M | 0,999999 | 1 eee 
| 0 ; Za if 
Einseitig eingespannter | 1 ~ | 1,875 10 | M 0,734096 | 1 | ENE 
Bi t - i | | 5 
iege ae | | | | aft 
xz=0 x=l | ; 
FE e094 09 Mole deoed, 2) gr 
Bg) = Gof x ee “2% | 3 | 7,85476| M | 0.999225 1 ee 
| | 
D(x) = Gin x — sin ae x| 4 10,9955 M | 1,000 034 ee. i 
| | 
Gin o, m — sin 0; a 
gt ae 5 | 14,1372 | M | 0,999999 | 1) poe 
aes | 


Beiderseits gelagerter 


| 
- - 4 
Biegetrager . A, = ow? = rt (7) = Ty(x) = sin + x 
@=0 ea) . 
oy re a | HM, = M/2 


M = Masse des Biegetragers; EJ = Biegesteife ; | = Lange des Biegetragers; m= M/l 


Greift in dem durch die Koordinaten %,, 2, gekennzeichneten Punkt des Kérpers eine in Richtung 
- der y-Achse wirkende harmonische Wechselkraft P,, an, und wird die in einem Punkt mit den Koor- 
dinaten x;, z; gemessene Schwingweg-Amplitude mit y(x;, z;) bezeichnet, dann ergibt sich die durch 


TE 
dip, = ¥(% jp 2;)/ Pr (39a) 
definierte dynamische EinfluBzahl gemaB (35) zu 
= jovi Xj Xk Zi Zk 
n= —F (rz a eee (39b) 


Analog dazu ergeben sich die EinfluBzahlen fiir harmonisch zeitlich veranderliche (homogenisierte) 
Momente M,,/l, baw. M,/I,, die um die x-Achse bzw. z-Achse drehend wirken, 


din RELA Oui und d;, 10h,’ (39 c) 
wobei dann definitionsgemaB = ke: 
y (x; 2) = Dd dy Pi, + 4px d) Ma/lo + 45, 5 Milly (40) 


wird. Selbstverstandlich hatte sich dieses Ergebnis auch aus den dynamischen Grundgleichungen 
herleiten lassen. 
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‘ ; ; ; ‘ , es 
y) Massefreie elastische Gebilde. Ist die HKigenmasse emes steifen ede tt: 
vernachlassigbarer GréBenordnung gegeniiber den iibrigen Massen meer. ee sa 
i u astisc 
kann man dieses Gebilde als masselos betrachten. Bei der Rechnung werden n 


Eigenschaften beriicksichtigt. 


Edy 21% 


Abb. 7, Stabférmiges Gebilde, durch Lingskriifte beansprucht. 


i ii i gestellte, als masselos voraus- 

1. Zug-Druckbeanspruchte Gebilde. Fir das in Abb. 7 dargestel , als 
gesetzte santa Gebilde, das durch die beiden Krafte P, und P, in seiner Langsachse ee 
wird, gilt mit den in die Abbildung eingetragenen Bezeichnungen, wenn die Federkonstante des 


Gebildes mit c bezeichnet wird, 


P,+ P,=0, P, = (42 — 41) 3 (41a) 
umgeformt erhalt man 
Py = eee (41b) 
Py = —¢q + C4, 
woraus sofort die Steifigkeitsmatrix 
1 —l ne 
Ce oer ata ( c) 


des massefreien Stabes entnommen werden kann. 

2. Biegungsbeanspruchte Gebilde. Zwischen den an einem Biegungstrager angreifenden 
n statischen Belastungskraften Q, und den an den Kraftangriffspunkten auftretenden Verfor- 
mungsgréfen y; bestehen lineare Beziehungen, die mit den statischen Steifigkeitszahlen h;, in 


bekannter Weise durch 
Q; = hy Yj (k = es oe ples n) (42) 
‘ 


ausgedriickt werden kénnen. 
Nun entstehen bei der Unterteilung eines Systems zumeist Gebilde, die frei von Bindungen an 


das feste Bezugssystem sind. Es liegen dann frei-freie Biegungstrager vor, fiir die im Sinne der 
Statik keine Einflu8zahlen, und damit keine Steifigkeitszahlen, angebbar sind. 
Bedenkt man jedoch, da8 es sich im vorliegenden Falle um schwingende Gebilde handelt, die 
auftretenden Krafte und VerformungsgréBen daher WechselgréBen sind, die wegen der voraus- 
gesetzten Masselosigkeit (und Dampfungsfreiheit) der Gebilde alle in Phase liegen, dann 1a8t sich 
sehr einfach eine Beziehung zwischen den Kraften und Schwingwegamplituden wie folgt herleiten 


Abb. 8, Frei-freier masseloser Biegungstriger. — a) Belastungsschema; b) Momentaniage des verformten Gebildes. 


Wir betrachten einen frei-freien Biegungstrager, an dem n gleichfrequente verallgemeinerte 
Wechselkrafte Q,, angreifen mégen (Abb. 8a). Die Schwingweg-Amplituden der Kraftangriffspunkte 
seien q,, die VerformungsgréBen y,. Da die am masselosen frei-freien Gebilde angreifenden Krafte 
in jedem Augenblick ein Gleichgewichtssystem bilden, kénnen wir zwei vollig willkiirlich heraus- 
gegriffene Krafte als Auflagerreaktionen, deren Angriffspunkte als (mitschwingende) Auflager- 
punkte, betrachten. 


ylidererenne pipes 
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Die Abb. 8b zeigt beispielsweise ein Gebilde dieser Art im ausgeschwungenen Zustand. Wenn 
wir die Punkte 1 und n als Auflagerpunkte wahlen, dann liegt, im Sinne einer Momentanbetrach- 
tung, ein zweifach gelagerter Biegungstrager vor, fiir den die Matrix der statischen Einfluzahlen 
@,; und, durch Kehrwertbildung, die Matrix der statischen Steifigkeitszahlen k;;, mit Hilfe der ele- 
mentaren Elastizitatsbeziehungen bestimmt werden kénnen. Fiir die n — 2 Krafte Gre ga rae ee die 
lassen sich dann die n — 2 Gleichungen 


P, =) kj Yj (k = 2,3,...n—1) (43 a) 
anschreiben. 


Die Verformungswege y; kénnen gemaiB Abb. 8b als Linearkombinationen der Schwingungs- 
wege 41, qn und q; , 


¥ =G—G/) a —Gl) an » (43b) 
dargestellt werden. Substituieren wir diese Gré®en in (43a), so ergibt sich 
n—l = n—l1 n—1 
Prat Shei Gl) + DS bei ti — In ey (LI)) » (43 ¢) 
= j=2 j=2 


woraus sofort die (n — 2)? dynamischen Steifigkeitszahlen 


C= Oe hy, 


ee n— 
eas 7 (43 d) 


Pus ene (k = 2, 3, ...n—1) (43 e) 
1 


erhalten. 
Zwei weitere Gleichungen folgen aus den Gleichgewichtssatzen der Statik, die wir mit 


n—l1 n 


P=— TPG), 3 P,=0 ames 


ansetzen wollen. 
Damit erhalten wir mit (43e) und (43d) ein System von n linearen Gleichungen in den q;, 
das durch 


n 


Py => Ce 0; Citas rlet 252. 71) bzw. B=Cq (43 g) 
J = 1 
ausgedriickt werden kann. =. 
Die Matrix € ist, wie sich leicht verifizieren laBt, symmetrisch zur Hauptdiagonale; die in (43d) 
noch nicht angegebenen EinfluBzahlen betragen dabei 


= Sai,  j=— Sil) i =23...—1), 
Sate oy (43h) 
Cee ve Chn (J,/2) ? eet saan 235 Ckn > 
k=2 k= 1 


und es gilt auBer der Symmetrieaussage cj, = c ; noch die Rechenkontrolle, daB die Zeilen- bzw. 
Spaltensummen von © verschwinden; d. h., 
S = = Se: paelelh =. A). A3i 
au le atl 2p ocsn) baw. 2 ci O26 n) (431) 
Wie man aus der Herleitung ersieht, enthalt die Matrix © in ihrem Kern die (n — 2, n — 2) Matrix 
der statischen Steifigkeitszahlen des zweifach gelagert gedachten Gebildes, wahrend die Rand- 
elemente als Linearkombinationen der statischen Steifigkeitszahlen erscheinen. 
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Diese Ergebnisse lassen sich sinngema8 auch fir den Fall anwenden, dal anstelle der Krafte 
von auBen her einwirkende Belastungsmomente vorhanden sind, wenn dabei beriicksichtigt 
wird, daB die VerformungsgréBen die Verdrehungswinkel der Stabquerschnitte sind. 


5. Dynamische Steifigkeitsmatrix einfacher Gebilde. Bei einfach aufgebauten Gebilden, die 
beispielsweise aus wenigen Einzelmassen und -federn bestehen mégen, bzw. durch solche Gebilde 
approximiert werden kénnen, 1aBt sich die auf die Kopplungspunkte bezogene dynamische Steifig- 
keitsmatrix unmittelbar aus den fiir das Gebilde geltenden Schwingungsgleichungen herleiten. 

Dieser Rechnungsgang kann mit Hilfe eines Eliminationsverfahrens in Form einer fir das 
numerische Rechnen vorteilhaften Matrizenoperation einfach und iibersichtlich durchgefiihrt 
werden. Es mége hier kurz angedeutet werden. 

Die Schwingungsgleichungen eines aus diskreten Massen und Federn aufgebauten Gebildes von 


f Freiheitsgraden kénnen durch 


(S* Bes) Q*) q* = ‘pis (44a) 
mit ©* = Steifigkeitsmatrix!, Y(* = Tragheitsmatrix, ausgedriickt werden. Der Klammeraus- 


druck hat dabei im Sinne unserer Betrachtungen die Bedeutung der dynamischen Steifigkeits- 
matrix des Schwingungsgebildes, so daB wir 


Ox =O 1 9% (44) 

setzen kénnen. 
Greifen am Gebilde n verallgemeinerte Krafte, wobei n < f sei, an, dann J4Bt sich durch eine 
passende Numerierung der Kopplungspunkte erreichen, da der Spaltenvektor 9% die Form 


Q, 0, 
Ox =| 0, =) aib ayee| Rea ks (4c) 


annimmt. Das Gleichungssystem (44a) kann in zwei Teilsysteme aufgeteilt werden; in Matrizen- 
schreibweise bedeutet dies ein Unterteilen der Matrizen gemaB 


(D\ _(CE* Gx) fq 
(0) (ee ee) (a) ws 
mit 2 
ies n fe 
n Qn +1 j f 
{2 Qn+-2 Cr™ | OF 
g ae : ? Qi = : und (ox x =n . (44e) 
In Uf os" re | 


Das Aufspalten liefert nach Ausfiihrung der Matrizenmultiplikation die beiden Matrizengleichungen 


O = Or g + CX* gy. 0 = O5* gq + CF q; » (44) 
aus denen durch Eliminieren von Fi 
O = — 6 CPyE1 q | (448) 


erhalten wird. Der Klammerausdruck stellt dabei, wie wir durch Vergleich mit der Definitions- 


gleichung (18) ersehen, die gesuchte, auf die n Kopplungspunkte b i eifi 
keitsmatrix des Gebildes dar: datetime ci nA oS 


© = OX — OX* (xX) Ox*, (44h) 


1 Die Steifigkeitsmatrix @* bedeutet hier die Matrix der i ifigkei i 
; statischen Steifigkeit ie i 
Weise aus dem System der Schwingungsgleichungen abgelesen werden (ees Scr aan agen ene 


2 Die an der Matrix (xx an. eschriebenen n 0 le die Spalte bzw. Zeile 
ig I dizes ns 1 n kennzeich en. bei de 1 
ni ° r d e€ 
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6. Experimentelle Bestimmung der Einflufzahlen. Es kann manchmal die Aufgabe vorliegen. 
daB die Kigenschwingungszahlen einer Anordnung zu bestimmen sind, von der einzelne Teilgebilde 
bereits baulich fertiggestellt vorliegen. In solchen Fallen kann es erwiinscht sein, die Ponaceiclisns 
den dynamischen Parameter der baufertig vorliegenden Glieder (besonders dann, wenn sie sehr 
komplexe Strukturen aufweisen), experimentell aus einer Schwingungsmessung zu ermitteln. Die 


durch 


dee 
jk Pp, 


definierten KinfluBzahlen ergeben sich beispielsweise fiir ein approximierendes Gebilde von f 
Freiheitsgraden gema8B (35) zu 


d, OW IO eV). SG) 
Bm)" ad | Yb Ag—2) 


Bringt man die rechte Seite dieser Gleichung auf einen gemeinsamen Nenner, so erhalt man im 
Zahler eine ganze rationale Funktion (f — 1)ter Ordnung in /, die als Produkt ihrer Linearfaktoren 
dargestellt werden kann 


Die GréBen Vaan ; haben dabei die physikalische Bedeutung von Antiresonanz-Kreisfrequenzen, 
wahrend die GréBen M;, als Ersatzmassen des Gebildes! aufgefaBt werden kénnen. Beide GréBen- 
gruppen stellen bezogene GréBen (bezogen auf die Punktepaare x;, x,) dar. Ihre Werte kénnen, 
unter Einschlu8 der Eigenwerte 7,, meBtechnisch verhaltnismaBig einfach ermittelt werden. 


7. Anwendungsheispiele. Die nachstehend angefihrten Beispiele, die verschiedenen technischen 
Bereichen entnommen sind, sollen die Anwendung des beschriebenen Verfahrens kurz erlautern. 


a) Einfach gekoppelte Gebilde. Die Resonanzbedingung fiir diese einfache Art der Kopp- 
lung (siehe Abb. 2) lautet nach (4) d,, + dy, = 0. Beispiele hierfiir sind das Zusammenbauen einer 
Verbrennungskraftmaschine mit einer Luftschraube, oder das Kuppeln zweier zweifach gelagerter 
Wellen zu einem durchlaufenden dreifach gelagerten Wellenstrang, u. a.m. Diese Falle sind in 
ihrer Anwendung allgemein bekannt, so daB es sich wohl eriibrigen diirfte, auf weitere Kinzelheiten 
einzugehen. 


b) Mehrfache Kopplung. a) Doppelrumpf-Flugzeug mit Leitwerk (Abb. 9a). Unter 
der vereinfachenden Annahme, daf der Doppelrumpf einseitig fest eingespannt sei, erhalten wir 
das in Abb. 9b bzw, 9c in zwei Projektionen dargestellte, in drei Teilgebilde aufgeléste System 
mit den an den Trennstellen eingetragenen Bindungskraften bzw. -Momenten. Die beiden Teil- 
gebilde J und 2 erscheinen hierbei als einseitig eingespannte, auf Biegung und Verdrehung bean- 
spruchte Kragtrager, wahrend das Gebilde 3 (Leitflache) ein frei-freies Biegungsgebilde darstellt. Bei 
geniigend grofer Steifigkeit kann dieser Teil unter Umstanden (zumindest fiir die Eigenschwingung 
1. Ordnung) als frei bewegliche starre Platte betrachtet werden. Das Kopplungsschema und die 
Verkniipfungsmatrix? @ sind in Abb. 9d bzw. 9e zu ersehen. 


B) Turbomaschinen-Fundament. Bei der Bestimmung der Eigenschwingungszahlen der 
in Abb. 10 in stark idealisierter Form dargestellten Turbomaschinenanlage mit Stahlfundament 
kann man die Zerlegung in verschiedener Weise vornehmen. 

Fiir eine orientierende naherungsweise Berechnung geniigt es unter Umstanden, die Gesamt- 
anordnung durch das in Abb. 11 gezeigte Ersatzgebilde zu approximieren, bei dem man nur die 
Elastizitat der Fundamentstiele und der Lagerkonstruktion beriicksichtigt, und den Fundament- 
tisch, einschlieBlich der darauf montierten Maschinengehause, als starres plattenférmiges Gebilde 
betrachtet. Dieses Ersatzgebilde gilt fiir die beiden Hauptschwingungsformen, Schwingungen in 
der Hochachse und Schwingungen in der Tischebene (allerdings mit unterschiedlichen Feder- 
konstanten) unverandert. 


1 Siehe etwa E. Hiibner, Forsch. Ing.-Wesen, 25 (1959) S. 148. 
2 Siehe FuRBnote 1, Seite 138. 
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GréBere Genauigkeitsanspriiche lassen sich durch eine starkere Beriicksichtigung des LY 


turellen Aufbaues der Anordnung erfiillen. ! 
Eine mégliche Form der Zeclegung des Gebildes ist in den Abbn. 12a und 12b in Aufri8 und 


Grundri®B dargestellt. Wenngleich man mit Riicksicht auf den hiermit verbundenen groBen Auf- 


Abb. 10, Idealisiertes Fundament einer Turboanlage. 


TTI WL 


7 


Abb, 11. Turboanlage: Stark vereinfachtes Ersatzsystem. 
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1 7 Abb. 9 (links). Doppelrumpf-Flugzeug. — a) Schematisierte Darstellung; 


Fas leal b) Aufgespaltenes System mit BelastungsgréBen im AufriB; c) desgi. im 
4 7 SeitenriB; d) Kopplungsschema; e) zugehérige Verkniipfungsmatrix. 


wand an Berechnungsarbeit in praktischen Fallen wohl kaum eine so weitgehende Aufspaltung 
vornehmen wird, soll dieses Beispiel doch zeigen, wie selbst auf erordentlich stark gegliederte 
Anordnungen rein formal der Berechnung zuganglich gemacht werden kénnen. 

Die als starr vorausgesetzten plattenférmigen Massen M , und My stellen dabei die auf den Funda- 
menttisch montierten Maschinengehause, einschlieBlich entsprechender Teile der Tischkonstruktion, 
dar. Durch die an diesen Massen angeordneten Federn wird die Biegungselastizitat der die Portale ver- 

_bindenden Langstrager, an denen die Gehausetragkonstruktion befestigt ist, beriicksichtigt!. 


1 Ebenso kann auch die Elastizitat der Rotorlagerung (einschlieBlich der Nachgiebigkeit des Olfilms in 
den Lagern) durch Anordnung zusatzlicher Federn im Ersatzgebilde leicht beriicksichtigt werden. 


XXIX. Band 1960 


ee Hiibner: Eigenschwingungszahlen zusammengesetzter Schwingungs-Systeme 149 


Alle ubrigen Teilsysteme, mit Ausnahme der starr gekuppelten Rotoren, sind homogene stab- 
formige Gebilde, deren Eigenwerte und Kigenformen leicht anzugeben sind. < 
; Die Langsschwingungs-Eigenwerte der Fundamentstiele werden im allgemeinen immer weit 
uber dem Bereich der interessierenden, im wesentlichen durch die Biegungselastizitat der Bauteile 
festgelegten Eigenschwingungszahlen der Gesamtanordnung liegen. Man lean dann die durch die 
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63-65 65 64 63 
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ANS 62 yh sf ONG? 

52 46 
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Abb. 12. Turboanlage: Aufgespaltenes System mit BelastungsgréBen fiir Schwingungen in der Hochachse. — a) AufriB; b) Seitenrif; 
c) Kopplungsschema (zur Vereinfachung sind mehrfach vorkommende Bauteile gleicher Art tibereinandergezeichnet dargestellt). 


Langskrafte beanspruchten einseitig eingespannten Stabe als masselose Federungsglieder betrach- 
ten, fiir die die Federungsmatrix nach (41c) wegen der festen Einspannung zu einem einzigen 


Element entartet. 
Der Vollstandigkeit halber ware noch zu erwahnen, daB die als starr gekoppelt vorausgesetzte 


Laufergruppe bei der Berechnung ihrer Kigenwerte und Higenschwingungsformen als frei-freies 
Biegeschwingungsgebilde anzunehmen ist. 


(Eingegangen am 14, Juli 1959.) 
Anschrift des Verfassers: Privatdozent Dr.-Ing. Erhard Hiibner, Aachen, Im Brockenfeld 17. 
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Bemerkung zur vollstandigen Auswertung des Spannungszustandes 
in Platten mit Hilfe der Spannungsoptik 


Von H. Bufler 


1. Grundlagen. Mittels der Spannungs-Dehnungs-Gleichungen (Hookesches Gesetz) bei Beriick- 
sichtigung der Geometrie der verformten Platte und der iiblichen Vereinfachungen (Kirchhoffsche 
Annahmen) lassen sich die Spannungen durch die Ableitungen der Durchsenkung W ausdriicken. 
Statt der Spannungen fiihrt man gewéhnlich die (durch Integration iiber die Plattendicke h ge- 
wonnenen) Schnittmomente und Schnittkrafte ein. Fiir die auf die Langeneinheit bezogenen Biege- 
momente M,, M, und des Torsionsmoment T ergeben sich dann in Bezug auf die kartesischen 
Koordinaten x, y die bekannten Beziehungen 


re ee " 
si aa 
B= ON" Seay @) 
wobei v die Querdehnzahl und 
Pe es “ 


die Biegesteifigkeit bedeuten. Es sei hier gleich vermerkt, daB diese Gleichungen auch bei verander- 
licher Plattendicke Giiltigkeit behalten. Der Zusammenhang mit den Spannungen o,, ¢, und tT am 
unteren bzw. oberen Plattenrand ist dabei gegeben durch 


Ox 6 M, 
of Ee (My (5) 


Aus den Gleichgewichtsbedingungen am Plattenelement h dx dy 1aBt sich eine Beziehung zwi- 
schen den Schnittmomenten herleiten, welche in Verbindung mit (1), (2), (3) die Differential- 
gleichung fiir die Plattendurchbiegung W liefert. Diese ist der Ausgangspunkt fiir eine theoretische 
Lésung von Plattenproblemen und wird fiir eine spannungsoptische Auswertung nicht bendtigt. 
Das Ziel der spannungsoptischen Auswertung besteht darin, die Schnittmomente M,, M, und T 
zu ermitteln, mit denen wegen (5) der Spannungszustand an der betrachteten Stelle vollstandig 
gegeben ist. 


Fihrt man die HilfsgréBen H, und H y gemaB 


H, ew 
ee i) yr? (6) 
H ew 

B -~—) sa (7) 


ein, so ergeben sich — wie man durch partielles Differentiieren von (6) nach x, (7) nach y und (3) 
nach y bzw. x leicht verifizieren kann — bei Beriicksichtigung von (4) folgende Beziehungen, welche 
das Torsionsmoment T und die GréBen H,, H y miteinander in Verbindung setzen: 


(i) + & (er) =°. 
> (i) a fa eae (9) 


(Bei konstanter Plattendicke kann man h wegkiirzen.) Diese beiden Gleichungen entsprechen 
iibrigens den Codazzischen Gleichungen in der Flachentheorie. 


XXIX. Band 1960 HH. Bufler: Auswertung des Spannungszustandes in Platten mit Hilfe der Spannungsoptik 51 


Die Biegemomente M, und M, lassen sich aus H, und.H, wie folgt berechnen: 


Mo — (HI, +H), (10) 


Bi Tecpeee v Hy). (11) 


2. Das Auswerteverfahren. Bemerkenswerterweise stimmen die Gleichungen (8) und (9) formal 
vollkommen mit den Gleichgewichtsbedingungen am Element dx dy beim ebenen Problem (Scheibe) 
tberein, wo fir H,/h°, H,/h® die Normalspannungen und fiir T/h? die Schubspannung stehen. 
Im ebenen Fall liefert der spannungsoptische Versuch die Linien konstanter Hauptschubspannung 
(Isochromaten, 6 = konst.) und die Linien konstanter Hauptspannungsrichtung (Isoklinen, Diz 
konst.); aus letzteren 14Bt sich das Netz der Hauptspannungstrajektorien konstruierem. In ent- 
sprechender Weise bedeuten bei Versuchen an Platten! (welche — damit sich der optische Effekt 
nicht aufhebt — entweder mit versilberter Mittelflache oder aus zwei verschiedenen zusammen- 
geklebten Platten (Zweischichtverfahren) hergestellt werden miissen) die Isochromaten Linien 
konstanten Haupttorsionsmoments (= Differenz der Hauptbiegemomente geteilt durch zwei) und 
die Isoklinen Linien konstanter Hauptbiegemomentenrichtung. Es ist daher naheliegend, fiir 
die vollstandige Auswertung des Spannungszustandes auf die Methoden der ebenen Spannungsoptik 
zuriickzugreifen. G. Haberland und H. Schwieger? verwenden fiir die Auswertung langs einer Ge- 
raden (sie sei die x-Achse) die dem Schubspannungsdifferenzverfahren analoge Methode, bei der 
neben der Auswertungslinie noch zwei dazu parallel verlaufende Hilfsgeraden, auf denen man 
im allgemeinen Fall der Platte mit veranderlicher Dicke den Ausdruck T/h? (und daraus % (Gs) 
_ bestimmen mu, benétigt werden. Das vom Verfasser zunachst fiir den ebenen Spannungszustand 


angegebene Verfahren®, das die direkte Bestimmung des in (8) erforderlichen Differentialquotienten 
4 ; 


ay 
wertungslinie auskommt. Damit werden die Hilfsgeraden iiberfliissig und infolgedessen auch die 
haufig mit groBen Ungenauigkeiten verbundenen Differenzenquotienten. 


IF : é : : 
(i) gestattet, besitzt den Vorteil, da man allein mit den Versuchsdaten lings der Aus- 


Im folgenden wird gezeigt, wie sich dieses Verfahren auch auf Platten (mit konstanter oder 
veranderlicher Dicke) iibertragen l4Bt. Wir gehen von der Hauptgleichung der Spannungsoptik, 
welche bei Platten 


‘ ge act = oh : (12) 


lautet*, aus. Im Gegensatz zum ebenen Fall steht in (12) die Dicke h im Zahler. Ferner gilt die aus 
dem Mohrschen Kreis folgende Beziehung 


deat sin 2 - (13) 


Aus (13) mit (12) erhalt man bei Beriicksichtigung der Veranderlichkeit von h durch partielle Dif- 
ferentiation 


a /T Sai ap Oh 26-4 
ait) = ana ay sin 2p + 5 2d cos2g— 55 sin 2). sew (1d) 
Wie ich gezeigt habe **, kann man schreiben 
00 oy 
eB a = Pa tw Oe (15) 


Dabei bedeutet a, bzw. a, den Winkel, den die Kurve 6 = konst. (Isochromate) baw. p = konst. 
_ (Isokline) an der Stelle x mit der x-Achse einschlieBt (Abb. la und 1b); die Werte fiir 0, p, 4, 


1 L, Féppl u. E. Ménch, Praktische Spannungsoptik, S, 116, 2. Aufl. Berlin/Gottingen/Heidelberg 1959. 
2 G, Haberland und H. Schwieger, Bauplanung u. Bautechnik 9 (1955), S. 350. 
3 H. Bufler, Z. f. angew. Phys. 8 (1956) S. 139. 
* Siehe G. Haberland, a. a. O. Gl. (42); dort wird die Isochromatenordnung 6 — 6 kann auch der Bruch- 
teil einer Ordnung sein — mit m, die spannungsoptische Konstante mit C (C = 1/S) bezeichnet. 
** Siehe FuBnote 3 
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kénnen direkt aus den spannungsoptischen Bildern entnommen und als Funktion von x aufgetragen 
werden (Abb. 2a und 2b). Hierbei sind f, bzw. B, die Winkel zwischen der x-Achse und der Kurve 
6 = 6(x) baw. y = g(x) und miissen auf graphischem Weg ermittelt werden. Die Ableitung dh/¢éy 


b 


Abb. 2. Die langs der Auswertungslinie benétigten Daten. 


a8 
S23 


1aBt sich aus dem Verlauf der Plattendicke bestimmen. Setzt man (15) in (14) ein, so folgt die # 
Formel 


Chwes i S ; CW.2'0—. 
ay (7) =~ gpa (ete ms te By sin 2g + ctg a, tg Bh, 20 cos 2m + me sin2 9). (16) 


Durch graphische Integration von (8) gelangt man dann zu 


Omaha: A, 
B= = fete | (7, (17) 


t= b 


wobei die mit dem Index 0 behaftete Integrationskonstante den Randbedingungen bei x = 0 ge- 
niigen muf*. Zur Berechnung von M, und M, ist neben der Kenntnis von H, gemaB (10) und (11) 


: ; (M,—» M,) (18) 


noch eine weitere Beziehung zwischen M, und M, erforderlich; diese liefert der Mohrsche Kreis: 


M,—M, = £2.T,,cos2@. (19) 


max 


H, =— 


z 


3. Zusammenfassung. Es wird gezeigt, wie man das vom Verfasser fiir den ebenen Spannungs- 
zustand angegebene Auswerteverfahren der Spannungsoptik auch auf Plattenprobleme iibertragen 


kann. _ 3 
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